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Karl Girkmann + 


Am 14. Juli dieses Jahres verschied 0. Prof. emer. Dipl.-Ing. Dr. techn. Dr. h. ec. 
Karl Girkmann nach langen, mit groBer Geduld ertragenen Leiden. 

Girkmann wurde am 22. Marz 1890 in Wien geboren. Schon als Student des 
Bauingenieurfaches an der Wiener Technischen Hochschule hatte er Gelegenheit, 
bei technisch bedeutenden Eisenbahn- und Briickenbauten mitzuwirken. Nach 
AbschluB des Stu- rechnung und Kon- 
diums arbeitete er se »  struktion — groBer 
als Bauingenieur in Stahlhoch- und 
der Privatindustrie. Briickenbauten des 
Er hatte Planungs- In- und Auslandes 
arbeiten fiir Hoch- verantwortlich. Ne- 
spannungsleitungen benher liefen aus- 
durchzufiihren und gedehnte Entwick- 
war an der Ent- lungsarbeiten werk- 
wicklung der elasti- stofftechnischer und 
zitatstheoretischen konstruktiver Art. 
Grundlagen des 1925 wurde er zum 
Leitungsbaues und Doktor der Tech- 


an der Ausarbeitung nischen Wissen- 
der _ einschlagigen schaften an der 
Vorschriften und Wiener Technischen 
Richtlinien mab- Hochschule promo- 
geblich __beteiligt. viert. 


Im Jahre 1931 
tratGirkmann, nach- 
dem eine Berufung 
an die Deutsche 
Technische Hoch- 


Als Vorstand des 
Technischen Biiros 
der Wanitsch-Hild- 
Werke und spater 
bei der Wiener 


Briickenbauanstalt schule in Brinn 
Waagner-Bir6é A.G. Kobé, Wien 20 der ablehnenden 
war er fiir die Be- Haltung der tsche- 


choslowakischen Regierung gescheitert war, zum Lehramt iiber und habilitierte sich 
1934 als Privatdozent fiir Stahlbau an der Technischen Hochschule Wien. Das Jahr 
1937 brachte ihm eine Berufung nach Graz, die er ablehnte, und eine zweite Be- 
rufung nach Wien, der er Folge leistete. Er wurde 1938 zum ordentlichen Professor 
und Vorstand des Instituts fiir Elastizitats- und Festigkeitslehre ernannt. Er leitete 
dieses Institut durch fast 20 Jahre, bis ihm ein schweres Leiden die weitere Austibung 
seiner Lehrtatigkeit unméglich machte. Im Jahre 1958 erfolgte seine Emeritierung. 

Im Studienjahre 1948/49 war Girkmann Dekan der Fakultat fiir Bauingenieur- 
wesen und im Studienjahre 1949/50 stand er als Rector magnificus an der Spitze 
der Hochschule. 1950 wurde er zum wirklichen Mitglied der Akademie der Wissen- 
schaften in Wien gewahlt. Die Technische Hochschule Graz verlieh ihm im Jahre 1955 
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in Anerkennung seiner hervorragenden wissenschaftlichen Leistungen das KEhren- 
doktorat. Er war Inhaber der Wilhelm-Exner-Medaille des Osterreichischen Gewerbe- 
vereins. Knapp vor seinem Tode erhielt er die goldene Ehrenmtinze des Oster- 
reichischen Ingenieur- und Architekten-Vereins. 

Girkmanns wissenschaftliches Werk ist auBerordentlich umfangreich. Im Jahre 
1938 erschien das gemeinsam mit dem Elektrotechniker Dr. E. Kénigshofer ver- 
faBte Buch ,,Die Hochspannungsfreileitungen“‘, dem bald nach dem Ende des zweiten 
Weltkrieges das Lehr- und Handbuch ,,Flachentragwerke“ folgte, von dem in kurzer 
Folge vier Auflagen erschienen. Fiir das englische ,,Civil Engineering Reference 
Book schrieb er ein Referat iiber die Elastostatik der Platten und Schalen. 
Daneben hat Girkmann noch iiber 40 Abhandlungen auf den Gebieten der Elasti- 
zitats- und Plastizitatstheorie, der Stabilitatstheorie und der Baustatik verfaBt. 
Uber die Ergebnisse seiner Forschungen hat er wiederholt auf internationalen 
Kongressen berichtet. Er war Mitherausgeber dieser Zeitschrift. 

Mit Professor Girkmann ist ein Ingenieur und Wissenschaftler von uns gegan- 
gen, der es wie kaum ein zweiter verstanden hat, Forschung und Lehre, Theorie 
und Praxis zu verbinden und gleichermafen zu pflegen. Seine Fachgenossen und 
seine Schiiler in aller Welt werden ihm ein treues Andenken bewahren. 


H. Parkus, Wien. 


Warmespannungen beim Abbinden von Massenbeton 
Von E. Tremmel, Salzburg 


Zusammenfassung: Die durch beliebige Verteilung von Warmequellen zeitlich abklingender 
Intensitét in Kreiszylindern und rechteckigen Prismen hervorgerufenen Temperatur- und 
Spannungsfelder werden allgemein abgeleitet. 

Im zweiten Teil wird die Anwendung der Theorie auf die Bestimmung der in einem Beton- 
zylinder zufolge der Warmeentwicklung beim Abbinden auftretenden Temperaturspannungen 
gezeigt. 


Erster Teil 


I. Einleitung 


1. In der vorliegenden Arbeit werden die von Warmequellenverteilungen zeitlich 
abklingender Intensitat in Kreiszylindern endlicher Lange oder unendlich langen 
rechteckigen Prismen hervorgerufenen Temperaturspannungszustinde untersucht. 

Instationire Temperaturfelder dieser Art treten z. B. beim Abbinden von Beton 
zufolge der Warmeentwicklung des Zements auf; in Anbetracht ihrer Bedeutung fiir 
den Massenbeton haben sich verschiedene Fachleute, wie F. Télke1, K. Hirschfeld2, 
A. Stucky und M. Derron® u.a., mit diesen Temperaturzustanden befaBt, ohne 
jedoch auf das mit ihnen verbundene Spannungsproblem einzugehen. 

Wir werden hier zunachst die rotationssymmetrischen bzw. ebenen Temperatur- 
felder fiir beliebige Warmequellenverteilungen unter Annahme homogener, zeit- 
unabhangiger Randbedingungen ableiten und uns hiebei einer allgemeinen, fiir rota- 
tionssymmetrische Zylinderkoordinaten und ebene kartesische Systeme giiltigen Dar- 
stellungsform bedienen. 


1 F. Tolke: Talsperren, Berlin. 1938. 
* K. Hirschfeld: Die Temperaturverteilung im Beton. Berlin. 1948. 
$ A. Stucky-M. Derron: Problémes thermiques... Lausanne. 1957. 
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Im weiteren werden die zugehérigen Spannungszustinde unter Anwendung der von 
E. Melan und H. Parkus‘ entwickelten Verfahren berechnet. Im rotationssymme- 
trischen Fall wird auf die umfassende Arbeit von R. Trostel zuriickgegriffen werden, 
wihrend das ebene Problem mit Hilfe einer vom Verfasser an anderer Stelle verdffent- 
lichten Analogiemethode* gelést werden kann. 

Die allgemeinen Berechnungsvorgange werden schlieBlich an Hand eines speziellen 
Beispiels erlautert, in dem der unter der Wirkung gleichmafig verteilter Warmequellen 
stehende Vollzylinder endlicher Linge, dessen Mantel- und Stirnflachen Warme in das 
umgebende Material abstrahlen, behandelt wird. 


2. Wir legen unseren Untersuchungen ein rechtwinkeliges kartesisches Koordinaten- 
system oder rotationssymmetrische Zylinderkoordinaten zugrunde und betrachten ein 
durch je zwei Koordinatenlinien oder -flachen begrenztes Rechteck oder zylindrisches 
Prisma, dessen Inneres durch eine Warmequellenverteilung exponentiell abklingender 


Intensitat Re 
W (a1, %2)e % (1.0) 


erfiillt sei. In diesem, fiir verschiedene Falle der Praxis in Betracht kommenden Gesetz 
bedeuten W (2, x.) die Quellverteilung zur Zeit ¢ = 0 und ¢, eine, den Grad des Ab- 
klingens kennzeichnende Zeitkonstante*. Die Differentialgleichung der Warmeleitung 
nimmt dann die Form 

1 08 Ww -- 


rae Oe el (1.1) 


an; mit a bzw. A sind hier wie tiblich Temperaturleitzahl bzw. Warmeleitfaihigkeit 
bezeichnet. 
Die homogenen Oberflachenbedingungen seien in der allgemeinen Form: 


eC, . 
ig in Xi = lip (1.2) 
gegeben, wobei der Index R (R = 1,2) auf die entsprechende Randflache x;r = konst. 


hinweist; unter mae ist der auf die auswarts gerichtete Oberflachennormale n;r bezo- 
gene Temperaturgradient zu verstehen. Die den einzelnen Oberflichen zugeordneten 
Konstanten h;r werden im allgemeinsten Fall fiir jeden Rand x; = /;r verschieden 
sein und kénnen auch die Werte 0 oder ~ annehmen. 


SchlieBlich gelte fiir die Anfangstemperaturverteilung @, mit ¢t = 0: 
Bo = Do (Ly, Xo). (1.3) 


II. Temperaturfeld: Homogene Liésung 


1. Die den Oberflaichen- und Anfangsbedingungen geniigende Lésung von (1.1) 
wird bekanntlich aus der Uberlagerung des mit passenden Koeffizienten multiplizierten 


allgemeinen Integrals 8 der homogenen mit einem Partikularintegral & der inhomogenen 
Differentialgleichung gewonnen. Setzen wir fiir die Lésungsfunktion der (1.1) zugeord- 
neten homogenen Differentialgleichung 


4 E. Melan-H. Parkus: Warmespannungen. Wien. 1953. 

5 R. Trostel: Instationare Wairmespannungen in Hohlzylindern mit Kreisringquerschnitt. 
Ing.-Arch. 24, 1ff. (1956). 

6 E. Tremmel: Zur Anwendung der Plattentheorie bei der Bestimmung von Warmespannungs- 
feldern. Osterr. Ing.-Archiv 11, 56ff. (1957). 

* Ein allgemeinerer zeitlicher Verlauf der Intensitaét kann durch Uberlagerung mehrerer 
analog zu (1.0) aufgebauter Funktionen erhalten werden. 
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Bake, (2.1) 
a 


Q 
~ 


die Méglichkeit einer Produktdarstellung voraus, dann gilt 


ae Oo) pS —a 9: ato m)é 
0 — Ss! Sy ee. es (P1,n X1) P.& es (P2,m X2) e€ ( oe 23 ) ’ (2.2) 


n=1m=1 


also 
o=0*e“ (hn tP2, me)! (2.3) 
worin die Funktionen X,, im allgemeinen Linearkombinationen unabhangiger Lésun- 


gen U;, (p; 2:), Viv (@iy %i) der mit dem Ansatz (2.2) nach Kiirzen der Zeitfunktion 
fiir 7 = 1,2 folgenden homogenen Differentialgleichungen 


A Xiv + giv Xiv = 0 (2.4) 
sind: 
X iy (Qiy 2%) = Ain Ui (Viv 2%i) ie biy Vi (Viv Xi)» (2.5) 


eee Ae 3 b; , 
Die Eigenwerte ¢;, sind, ebenso wie die Verhaltnisse ms der Integrationskonstan- 


UG 


ten, durch die Randbedingungen bestimmt, die mit (1.2), (2.2) und (2.5) fiir jede Rand- 
flache a;z = konst. = l;r, R = 1,2 und jedes Reihenglied y eine homogene Gleichung 
der Form: 


OX; 
ONiR 


Viv [Giv Uy (giv liz) + Biv Viv (iv lir)] + 


xi=lR 
+ hir [div Uiy (giv liz) + bv Viv (viv liz)] = 0 (2.6) 


liefern; die Ableitungen nach dem Gesamtargument (9;, 2;) sind hiebei mit ’ bezeichnet. 
Werden fiir jedes Reihenglied die mit den Konstanten a;, und b;, verbundenen 
Glieder zusammengefaBt, wobei wir annehmen, daB l;, > x; > 1;; und daher 


Ox; sail 4rd fiir ly» 
OR Me Ae firs Is 


gilt, da die Ableitung nach der auswirts gerichteten Normalen positiv ist, so erhalten 
wir 
diy [Piy Uiy (Pin U;,9) + hi,2 Vin (Piv U;,2)] + 
+ Os, [Piv Vin (Piv I; 2) == h;,» Vis (Viv l;2)] = 0 (2.7) 
diy [— Piv Vin (iv li,1) + Iy,1 Usy (Qin 4,1.) + 
+ 63, {= Viv Ve (Viv L; 1) ai h;,1 Wes (Viv l;1)] = 0. (2.8) 


Aus den verschwindenden Koeffizientendeterminanten dieser zweigliedrigen 
Gleichungssysteme folgen mit 
[piv Ui (Piv b:,2) + hi,2 iv (Pin Ui,2)) [— iv Viv (iw G1) + Ret Vi (Give. — 
— [= Pin Vin (Pin G1) + hi,1 iv (Viv Ui. )] [Qiv Viv (Piv b:,2) + hi,e View (Piv4i,2)] = 9 
(2.9) 
vy co Bestimmungsgleichungen fiir die Eigenwerte ;,; sind diese bekannt, dann 
b; 


kénnen schlieBlich auch die Verhaltniswerte ay berechnet werden. Der belanglose 


| 
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Proportionalitatsfaktor a;, wird in der Folge unterdriickt. In bestimmten Fallen, wie 
etwa bei symmetrischen oder antimetrischen Randbedingungen mit /;, = —J;.. und 
h;,1 = +h;,. an den geradlinigen oder ebenen Begrenzungen, oder wenn, wie beim 
Vollzylinder an Stelle einer Randbedingung die Endlichkeitsbedingung tritt, verschwin- 
det eine der beiden Konstanten a,,; die zweite Konstante 6,,, der dann der Wert 1 
erteilt werden kann, tritt damit gleichfalls nicht mehr in Erscheinung. 

Damit sind also die im Integral (2.2) der homogenen Differentialgleichung (2.1) 
auftretenden Funktionen X;, vollkommen bestimmt; iiber die Reihenkoeffizienten 
K,,, wird bei der Erfiillung der Anfangsbedingungen verfiigt werden. 


2. Koordinatensysteme: invariante Darstellung der mathematischen Beziehungen. 
Die Funktionen X,,,, (91,. 21) X2,m (P2,m 2) erfiillen die Orthogonalitatsbedingung: 


B 
as (MUR i== 0; Ha In 
fs. (P1,n %1) X1,v (G1,v %1) Xe, m (P2,m X2) Xeu (P2,n ave S : 
|= 0 fiir » +n, w+m; 
(2.10) 
das Integral ist hiebei iiber das gesamte Grundgebiet B zu erstrecken. Das in den Koor- 
dinaten x; ausgedriickte Inhaltselement ist durch 
_ 9 (X,Y, 2) _ 
dy ni) (%1, %2, X3) day diy dts (2.11) 
gegeben, wobei x, y, z kartesische Koordinaten bedeuten und die Funktionaldetermi- 
nante in tiblicher Weise mit dem symbolischen Bruch bezeichnet wird. Fiir kartesische 
Koordinaten hat sie demnach den Wert 1, wahrend sie bei dem auf Zylinderkoordinaten 
bezogenen Inhaltselement dy = r dg dr dz gleich r wird. 
Da wir hier nur ebene und drehsymmetrische zylindrische Gebiete betrachten, 
kann im ersten Fall die Koordinate z, im zweiten die Koordinate g auBer acht gelassen 
werden. Damit entspricht 


- { dx:dy dem Flachenelement in ebenen kartesischen Koordinaten 
’ ~\rdr-dz dem Raumelement in drehsymmetrischen Zylinderkoordinaten. 


Um die Untersuchungen in der Folge in einer fiir beide Koordinatensysteme giiltigen 
Darstellung vornehmen zu konnen, werden wir das Inhaltselement in der Form 


dy =da;-d x (2.12) 


einfiihren, womit die bei den Zylinderkoordinaten durch den _ ,,Gewichtsfaktor“ r 
bedingte Auszeichnung der Koordinate r zunachst nicht in Erscheinung tritt. 

Mit (2.12) kénnen wir nun die Orthogonalitatsbedingung nach zwei Veranderlichen 
aufspalten und gelangen so zu der fiir beide Systeme giiltigen Darstellung 


x x Ft ee ee ne 
/ iv (Qiv 08) in (Pix x) Seo tae oe ae ( ) 

Die gemeinsame Behandlung der auf kartesische und Zylinderkoordinaten bezogenen 
Temperaturprobleme erfordert ferner die Beibehaltung des Laplace-Operators auch bei 
den aus (1.1) entspringenden gewohnlichen Differentialgleichungen. Fiir eine in karte- 
sischen Koordinaten ausgedriickte Funktion F' (2;), die also nur von einer Veranderli- 
chen abhangt, gilt dann mit 


oe he 

A= sat ane (2.14) 
ar 

AE (aes) Jaf * (2.15) 
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Ist die Funktion F (2;) dagegen auf rotationssymmetrische Zylinderkoordinaten 
bezogen, dann wird mit 


a sae a 
A= or a r or ot 02” ay 
aon. 1. oe a. e 
or Seater fir Say 7 (2.17) 
AF (x)= oe 
for es 27 (2.18) 
oz? 


so dafS mit der Darstellung etwa einer homogenen Differentialgleichung in der Form 
AF+6F=0 (2.19) 


alle hier in Betracht kommenden Moglichkeiten erfaBt werden. 


III. Temperaturfeld: Partikulérlésung 
1. Zur Berechnung eines partikuliren Integrals der inhomogenen Differential- 


gleichung (1.1) entwickeln wir die Storungsfunktion ae in einenach den Eigenfunktionen 


X,.n X>,m fortschreitende Doppelreihe, deren Konvergenz wir voraussetzen wollen. 
Unter Einfiihrung einer der Warmeerzeugung entsprechenden Konstanten W,, mit der 
W = W,-w (2%, X2) gelten soll, wobei w eine dimensionslose Ortsfunktion bezeichnet, 
setzen wir demnach 


=< = Cum X1,n (Y1,n 21) X2, m (P2, m X2) (3.1) 


mit den Entwicklungskoeffizienten 
Cig = Aya Aaa (3.2) 


fiir die wir aus (3.1) zufolge der Orthogonalitat der Eigenfunktionen 


°B 
| Xn Xemdy 
Ca 


B (3.3) 
| xi. adr 


erhalten, wobei die Integrale — wie schon erwahnt — iiber das ganze Gebiet zu er- 
strecken sind. 


Wo @ te vy so 


an aad 


Das Storungsglied kann damit durch —— nmA1,nX2,m ausgedriickt werden. 

Aus (3.3) ist ferner zu ersehen, daf die in (3.2) formal angedeutete Aufspaltung der 
Koeffizienten C,,,, nach den Faktoren 4,,,, Ag, m nur dann moglich ist, wenn die gegebene 
Warmequellenverteilung Wy: w (x,, 22), also die Funktion w auf die Produktform 


W = Wy (Ly) * We (Xp) (3.4) 
gebracht werden kann. 
Wird nun das partikulare Integral mit den noch zu berechnenden Reihenkoeffi- 
zienten Gyn, gleichfalls nach den die Randbedingungen erfiillenden Eigenfunktionen 
entwickelt, dargestellt: 


ri) = un: ~ to ee >» Gh Xa Xan (3.5) 


oe peas 
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oder 
B= Fett (3.6) 
mit der Ortsfunktion 
Ge Cia Xa hs ee (3.7) 
n=1m=1 
so folgt aus (1.1) 
Age + ~ 9 + Hey 0 (3.8) 


und damit unter Beachtung von (3.5), (3.1) und (3.7) nach Kiirzen durch 22 lig 


1 
Gam E (Xa DET) te Fin Xa te On Ee Xie Gn a 0; (3.9) 


fiir die Laplace-Ableitung des Produktes X,,,: X,,m ergibt sich aus der homogenen 
Differentialgleichung (2.1) nach Kiirzen der Zeitfunktion e774 ¥':+m)# ; 


A DPE Kees == (p71, +- (2, m) D.Ga Kain (3.10) 


Kinfiihrung von (3.10) in (3.9) liefert schlieBlich die Bestimmungsgleichungen fiir 
die Entwicklungskoeffizienten G,, zu 


1 
Cam (see — Pun — Pan) + Cam = 0, (3.11) 
aus denen 
Gam = 1 Cum x (3.12) 
a P71,n — P22,m 
folgt, womit 8* durch 
= Wowew Onm 
o* = — eieek > a hee Wee (P1,n X1) X2,m(P2, m V2) (3.13) 
n=1m=1——_ — — , — y2,m 
a 


bestimmt ist. Wird (3.13) in (3.6) eingefiihrt, dann konnen wir das vollstandige Integral 
von (1.1), das sich aus # = 9 + #in der Form 


os Eee ey pe ee Xin (P1,n X1) X2, m(Y2,m 2) + 


iiccal G71,n — P22,m 


sie 2 DOES Xijn(P1,n 1) X2,m(P2,m V2) e Frnt F'.m)! (3.14) 
n=l1m=1 
ergibt, anschreiben und haben dann noch aus der Anfangsbedingung die Reihenkoeffi- 
zienten K,m zu berechnen. 


2. Bei einer allgemeinen Warmequellenverteilung W (x,, x.) erhalten wir demnach 
das partikulare Integral der Differentialgleichungen (3.8) bzw. (1.1) in der nach Doppel- 
reihen entwickelten Form. 

LaBt sich aber die gegebene Quellverteilungsfunktion W >: w (x, a) nach den 
beiden Variablen x, und 2, gema8 (3.4) aufspalten, dann kann die in (3.13) auftretende 
Doppelreihe, wie nachstehend gezeigt wird, durch einfache Reihen ersetzt werden. 
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Mit 


W = Wy (X1) We (X2) (3.15) 
also 
W = Wy, (X13) We (2) (3.16) 
und dem (3.6) analogen Ansatz 
B= I* tlle 


folgt aus (1.1) nach Kiirzen der Zeitfunktion e~“ 
AG* + J F* + 72 wy (21) we (#2) = 0. (3.17) 


Fiihren wir (3.15) in (3.3) ein, so ergibt sich 


B 
[v1 X12 Xe,m dy 
Ay pt Aa = Cy = B : (3.18) 


[Xn X22 m dy 


und wir kénnen nach Aufspaltung des Inhaltselementes die Integrale nach den beiden 
Veranderlichen x, und x, getrennt auswerten und erhalten fiir 7 = 1,2 die Koeffizienten 


lia 


{wi (ai) Xiv (Yiv a1) A xi 


Ay= 4 (3.19) 
[Xv (iv i) d x; 
ly 
der beiden Reihen 
WwW; = Om Aiy Xe. (Qiv i) 5 (3.20) 
vi 
Setzen wir nun 
= Huy, (x4) 3) Ain Xiv (Givi); (3.21) 
v—21 


entwickeln demnach die Quellverteilungsfunktion nur nach den Eigenfunktionen 


Xiv (viv x) und treffen fiir den ortsabhangigen Term o* des Partikularintegrals den 
die Randbedingungen fiir x; = l;r erfiillenden Ansatz 


oe = =. DDN Xiy (Div Xi) Fey (xx), (3.22) 
Ved 


in dem neben den noch zu bestimmenden Funktionen F’,, (a) und den zugehérigen 
Integrationskonstanten D;, gleichfalls nur die von xz; abhangigen Eigenfunktionen 
X;y (Py vi) auftreten. 


Mit (3.21) und (3.22) geht (3.17) iiber in die Differentialgleichungen 
v=1 


co if il [o-e) 
A ind Vee Xiv (Qiv Xi) Fey (Xx) + aa Diy Xiy (Div Xi) Fey (Xz) = 5 
= v=1 


+ wy (2%) DY” Aiv Xiv (Viv 2:1) = 0, (8.23) 
eth 


aus der die Funktionen D,, F’,, (a) zu berechnen sind. Die Hinfiihrung dieser Funktionen 
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in die Randbedingungen liefert schlieBlich die Bestimmungsgleichungen fiir die Inte- 
grationskonstanten D,,. 


Zufolge der Orthogonalitat der Eigenfunktionen X,, (g;, 2;) mu8 (3.23) fiir jedes 
Reihenglied erfiillt sein; es gilt daher, wenn noch jede der aus (3.23) folgenden Diffe- 
rentialgleichungen durch D;, dividiert wird: 


1 ; 
A [Xin (Viv Xi) Fev (t%)] + 7 X iv (Viv %;) Fey (t%) + wy (2) SX (Giva;) = 0: (3.24) 
Aufspaltung des ersten Gliedes liefert 
A (Xi Fyy) = HA Xie kV AP, (3.25) 


wobei hier auf das Anschreiben der Argumente verzichtet wird. Fiihren wir fiir die 
Laplace-Ableitung der Funktion X,, den aus (2.4) folgenden Wert — q?,, X;, in (3.25) 
ein, so ergibt sich weiter 


A (X;, Py) ae (= Priv Fy Ae A F,,) Xiy (3.26) 


und wir erhalten nach Einsetzen von (3.26) in (3.24) und gleichzeitigem Kiirzen von X,, 
die gewohnlichen Differentialgleichungen fiir die Funktionen F,,: 


A Fyy (x) + (= = piv) Fev (%) =++ a Wea) — 0), (3.27) 
Mit der abkiirzenden Bezeichnung 
aie = pig ty (3.28) 
— wir beachten, daB die Werte a;, mit g?,;, > A rein imaginaér werden — folgt aus 
(3.27) : 
A Frey (ae) + ary Fay (st) + 5 we (ae) = 0. | (3.29) 


Der homogene Anteil dieser Differentialgleichung ist vollkommen analog zu der 
Differentialgleichung (2.4) aufgebaut und fiihrt daher auch auf die von 2, abhangigen, 
gemaB (2.5) definierten Losungsfunktionen, wenn dort der Index 7 durch k ersetzt wird, 
wobei aber hier an Stelle der Eigenwerte ¢,,, die durch (3.28) erklarten, von den Eigen- 
werten g;, abhangigen Werte a,, treten. Wir werden daher von jetzt ab die Bezeichnung 
andern, indem wir X,, an Stelle von F,,, also 


Fy (Xk) = Xpy (aE) (3.30) 
schreiben und die Lésung des homogenen Anteils von (3.29) hier mit 
Ee (Giy rk) = Xi (Giy Xx) (3.31) 
bezeichnen. Das vollstandige Integral von (3.29) kann nun durch 
y Aiy 
Xv = Xev (Gv TH) + ae fe (x) (3.32) 


ausgedriickt werden, worin wee (a;, x) analog zu (2.5) durch 
Xi (Gy Lk) = Ay Oxy (Giy Xk) + Bey View (Gy x) (3.33) 


gegeben ist und as fx (a) ein partikulares Integral von (3.29) darstellt, das, wenn ein- 
fachere Verfahren nicht zum Ziel fiihren, jedenfalls durch Variation der Konstanten — 


hier also der Koeffizienten a,,, 6,, — gewonnen werden kann. 
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Die noch unbekannten Konstanten a,,, b,, und D,, werden nun aus den Oberflachen- 
bedingungen berechnet. 
Mit f, 


0% = D;, Ay (Viv %;) Xky (Giy Xk) (3.34) 


nehmen diese nach Kiirzen durch X,, die Form 


diy X' py (iv kr) an + hr Xkv(Givkkr) = 0 (3.35) 
xk =Tkr 
an. 

Wir fiithren hier X,, und X’;v gem4B (3.32) und (3.33) ein, wobei die Ableitung nach 
dem Gesamtargument @;, 2, aus (3.33) zu berechnen ist. Mit den beiden Randflachen 
a, = I,r zugeordneten Koordinaten a, = i,,, und 4% = ,, erhalten wir zur Bestim- 
mung der je drei Unbekannten a,,, b,,, D;, je zwei homogene Gleichungen, aus denen 
wir zunachst das Verhdltnis 6,,/a,, der beiden Koeffizienten berechnen; wieder wird 
der dabei auftretende belanglose Proportionalitatsfaktor unterdriickt. Damit sind die 
Funktionen X,, (a;, 7%) bestimmt, und man erhalt die Unbekannten D;, aus einer der 
beiden Gleichungen (3.35) unter Beriicksichtigung von (3.32): 


= OX’ 
De \ Gey X' ky (av kr) Gir 


+ ir Xiv (Gi ha) + Aiv Ps bl + her fe (xx)| =O (3.36) 


ONR 
xk= kr xk =lkR 
VAN 
o ¢ 
fe (me) + her fk (xs) | a 
Dy = —Aiy —— aa ans : (3.37) 
diy X’ky (Giv LR) Gap + her Xkvy (Giv kr) 
xk = Ikr 


oder unter Einfiithrung der abkiirzenden Bezeichnung 


ote (ax 
Ndi) + her fk (xk) | 
r; eee ier R xk=lkR 
iv a an = (3.38) 
diy X’ky (aiy ler) fae + her Xkvy (aiv kr) 
xk=lkr 
mit 
Diy 

dy — Ay . (3.39) 


Setzen wir nun (3.32) unter Beachtung von (3.39) in (3.22) ein, so folgt fiir das 
Partikularintegral der Differentialgleichung (3.17) und damit auch der Differential- 
gleichung (3.8) die Funktion 


- Wo = 
o* = D> Ain Xin (piv %:) Peder Xe (ey an) + fe (ar) ] , (3.40) 


vol 


die also durch eine einfache, nach den EKigenfunktionen X;, (viv x;) und den von den 
Kigenwerten Pir abhaingigen Funktionen X;,,, (@iy x) fortschreitende Reihe aus- 
gedriickt erscheint, deren Konvergenz im Einzelfall noch zu iiberpriifen ist. 


Da aber in den bisherigen Untersuchungen keine der Veranderlichen ausgezeichnet 
wurde, mu8 die aus (3.40) durch Vertauschung von i bzw. v mit k bzw. u hervorgehende 
Reihe, die also mit dem Summationsindex ~ nach den Eigenfunktionen Xiu (Pen Le) 
und den Funktionen X;,,, (Qi, Xi) mit need 
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Oey = — — Wey (3.41) 
entwickelt wird, gleichfalls eine Lésung von (3.17) und (3.8) sein. 


Wir sehen, da8 die Darstellung partikularer Integrale von (3.8) durch einfache 
‘Reihen, sofern sie tiberhaupt méglich ist, im allgemeinen auf zwei verschiedene Arten 
erfolgen kann, und zwar in Abhangigkeit von dem jeweils in Reihen entwickelten Faktor 
der in Produktform gegebenen Stérungsfunktion w = w; (2;) w; (a). Fiir das voll- 
stindige Integral der Differentialgleichung (1.1) gelten daher neben (3.14) die gleich- 
wertigen Darstellungen: 


ne We ent Sas) Sa, (91, 01) [din X2,n (rn 2) + fe (ae) ] + 
ts 2, DY Kam Xan (01.01) Xayml(pem ayes Se ease be 
und 
= 8 elt SI ee (2, m X2) [de,m X1,m (G2,m %1) +fr (a)] + 


m=1 


ae >, a Kor, Aath (P1,n X1) Xo, m(P2,m X2) E26 on + P'2+m)e (3.43) 


n=1m=1 


IV. Temperaturfeld: Gesamtlésung 


1. Zur endgiiltigen Festlegung des Temperaturfeldes sind die noch unbekannten 
Koeffizienten K,,,, den Anfangsbedingungen anzupassen. Legen wir der Bestimmung 
dieser Konstanten zuniachst die Darstellung (3.14) zugrunde, dann folgt gemaB (1.3) 
fir += 0 


Ww Sa <= Cam 
Do{ i, L2) = — 3 Dy 2 1 Xin (P1,n U1) Xe, m (P2,m Xa) + 
ni m=1—_ — gn — Y2,m 
ato 
ue ay > he DO (P1,n X1) Xo,m(P2,m X2) (4.1) 


n=1 m=1 


Wird die Moglichkeit einer Entwicklung der Anfangstemperaturverteilung nach den 
Eigenfunktionen vorausgeésetzt : 


Do (a1, Xz) = Pa pains Xin X2,m, (4.2) 
n=1m=1 
mit 
s 
[80X11 X2,mdy 
BAS =S5 B ’ (4.3) 


{ X27 X22,mdy 


so ergibt sich aus dem Koeffizientenvergleich 
ee 1 Cum + Nove: (4.4) 


2 ee yy) 
E> ey 1l,n 2,m 
ate Pl, Pp 


Di ne 
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und damit schlieBlich fiir die gesuchten Reihenkoeffizienten Kym 


On m 


Kam = Lass ie - 1 (4.5) 


ate = P*1, n — Y?2,m 
Mit (4.5) und (3.14) erhalten wir fiir die vom Warmequellenfeld unter den gegebenen 
Rand- und Anfangsbedingungen hervorgerufene Temperaturverteilung 


Wo —t/t al y Cum a 
@=— ee DS DAG Os) A AP cee 
=1m=1 = 2 2 
iL m ato Prin P*2,m 
. y Wo O 2 2 
a, Een te | 2 Ot) a Pz 8) a 
a= a ee. ae 


(4.6) 


Bei verschwindender Warmeentwicklung geht (4.6) in das bekannte Gesetz iiber, 
das den aus einer gegebenen Anfangstemperaturverteilung folgenden Temperaturzu- 
stand beschreibt. 


Liegt dem Partikularintegral der inhomogenen Differentialgleichung (3.8) aber eine 
der beiden Darstellungen (3.42) oder (3.43), allgemein also die Form (3.40) zugrunde, 
so haben wir, um den Koeffizientenvergleich gemaB (4.4) vornehmen zu konnen (3.40), 
in eine nach beiden Higenfunktionen Xj, (vy 7) > Xzu (Yu Vi) fortschreitende Doppel- 
reihe zu entwickeln, deren Koeffizienten M,,, in nachstehender Weise zu berechnen 
sind: 


Wir stellen das Patrikularintegral (3.40) unter Weglassung des konstanten Faktors 


“i nach den Eigenfunktionen entwickelt dar: 


Dy OMT Sertioe en Garten LD Xing iv Xi) | Kew (ain ax) + 3 fe (wx) (4. 7) 


v=l1p=1 


beiderseitige Multiplikation mit X,, (yy 71) Xu (Pen x) dy und Integration iiber das 
Grundgebiet liefert: 


B 


MM, [es (Piv Xi) rm (Pkw rr) dy a 


B 
ep. [X% (piv i) Xgu (Pay %e) | Xev (ass eack < i. (x) | dy (4.8) 


Zufolge der Verkniipfung der Funktionen (3.32) 


nae | Kev (ain m) + a fh (x)| (3.32a) 


mit den Kigenfunktionen Xiv (piy X;) verschwinden bei der gemaB (4.8) vorzunehmenden 
Integration simtliche Reihenglieder, deren Indizes von » verschieden sind, also auch 
die von x abhangigen. 


Um das verbleibende Integral (4.8) auswerten zu kénnen, beachten wir, daB die 
Funktion X,,, (vy vi) die (2.4) entsprechende Differentialgleichung: 


A Xi, igen =O (2.4a) 
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und X;, = F;,, die (3.29) entsprechende Differentialgleichung: 
A Xv + ay Xey + uk (7,) = 0 
befriedigt und sowohl X,,, wie X;, fiir a = ir den Randbedingungen 
st + Ir Xe=0 (1.2a) 
genugen. 


Wird (2.4a) mit X;,, (3.29) mit X zu Multipliziert, so ergibt sich nach Gleichsetzen 
der so entstandenen Ausdriicke und einfacher Umformung die Beziehung 


XivA Da = Xin A Xiv = (@ iy = Pen) Xiv X ky + anu O.G (4.9) 


Multiplizieren wir hier beide Seiten mit dem weiteren, in (4.8) auftretenden Faktor 
X*:, (~;y #;) und integrieren iiber das Grundgebiet, dann verschwindet die linke Seite 
und es verbleibt 

B 


B 
r= (a?:y — pen) | Xr eA Bg Ak dy + i | x4. Wr Xgudy, (4.10) 


da die auf der linken Seite vorgenommene Integration fiir den a,-abhangigen Anteil 
gleichbedeutend mit der Anwendung des Greenschen Satzes ist und damit auf 


(4.11) 


0X, OX ky |*k= lke 
Xev ky Xie a 


OnR OnMR |\xr=Ik.1 


fiihrt; dieser Ausdruck aber wird wegen (1.2a) an beiden Grenzen zu Null. 
Wird nun in (4.10) fiir d;, aus (3.39) eingesetzt, so erhalt man fiir die rechte Seite 
von (4.8) 
B z 
Di»{ Dh XG Xin dy = — 5 | Xv Wr Xin dy: (4.12) 


@iy— Pew, 


und damit schlieBlich fiir die gesuchten Entwicklungskoeffizienten M,,, aus (4.8) 


B 
ie [Xv weXkp dy 


M,,=— (4.13) 


iy — Pky B 
[Xv Xk dy 


Der rechtsstehende, die Integrale enthaltende Bruch aber ist identisch mit A,,, 
also den Koeffizienten des nach den Eigenfunktionen X;,,, (~j, %,) entsprechend (3.20) 
entwickelten «,-abhangigen Anteiles des Stérungsgliedes w, (%,). Setzen wir noch im 
Nenner fiir @?,, aus (3.28) ein, so folgt mit n = », m = uw 

Man = — =" _ =- 5 ies (4.14) 


SN pe Oe ae De 
ee Q*in QP*km ate Prin P°km 


und das gleiche Ergebnis erhalt man auch, wenn in (4.7) die Indizes 7 mit & und y mit 
jw vertauscht werden. 

Multiplikation von (4.14) mit as und Einsetzen in (4.4) fiihrt daher — wie es sein 
muB — wieder auf die durch (4.5) festgelegten Entwicklungskoeffizienten Kn. 

Aus dem Vergleich von (3.40) mit (3.13) folgt daher unter Beachtung von (3.2) und 
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Ak 


(3.39), daB die Summe der Reihe =——*—— X;,, durch 
—— — Piv— Pkw 
ato 
Sy aaa pc ee ee i, ae 4.15 
ees, 1 kw y vey kv ( A ) 
p=1 ae _—_ Priv — Prkw 


gegeben ist, wie mit (3.42) und (3.43) bereits vorweggenommen wurde. 


V. Spezielle Koordinaten 


1. Nachstehend werden nun die aus der Spezialisierung auf kartesische und Zylin- 
derkoordinaten folgenden Funktionen. angegeben. Wir werden uns dabei allerdings im 
wesentlichen auf die Darstellung der Eigenfunktionen, die zur Bestimmung der Kigen- 
werte dienenden Gleichungen sowie auf die Angabe einiger aus der Aufspaltung des 
Stérungsgliedes folgender Beziehungen beschranken. 


2. Kartesische Koordinaten (ebenes System) 
Wir ordnen der 
Veranderlichen x, = x die Kigenwerte 1,1 = Yn ZU 
Veranderlichen xz, = y die Eigenwerte 92,m = Ym ZU. 


Zufolge (2.15) nehmen die Differentialgleichungen (2.4) hier die Form 


d2 BXGsn 
d xj? 


+ piv Xiv = 0 (5.1) 
an, deren allgemeine Integrale (2.5) durch 


Ay = An = Gyn COS (Oy 2) 4 O;. 5 SIN (Gs 1) (5.2) 
und 
oar — de = Qo.m cos (Wm y) te Osx sin (Wm y) (5.3) 
gegeben sind. 
Mit den Koordinaten 2;r der achsenparallelen Randlinien 
Lae 
hye und hy, 
le >y >t 
h,, und h,, 


sowie den zugehérigen Konstanten 


beziehungsweise 


erhalten wir die gemaB (2.9) gebildeten Bestimmungsgleichungen fiir die Eigenwerte 
gn im allgemeinsten Fall zu: 


[— @n* Sin (Gn bee) + Hey > COS (Qn lx2)] [— Gn * COS (Pn ley) + fey Sin (Gp ly)] — 
— [+ n° Sin (Gn lx) + he1 + cos (Gn ler) 1 [+ Gn * COS (Pn Ixe) + hee > sin (Pn lxe)] = 0 
(5.4) 


und analog mit y und Ym an Stelle von x bzw. ¢, fiir die Begrenzungslinien y = konst. = 
= 1,,,, l,, die Bestimmungsgleichungen fiir die Eigenwerte yn. 


Fir (3.33) ergibt sich mit 2, = 2, » = m sowie den gemaB (3.28) berechneten Werten 


1 
Tae a pn = ato ae Te (5.5) 


X,m (Bn &) = Xn (Bn €) = G,m COS (Bn 2) + by, m SiN (Bm 2); (5.6) 
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bei wachsenden Eigenwerten y,, wird mit m > Mo 


Bm? = —— — Wn? < 0 
und damit 


Anes re on 


wobei die imaginare Kinheit Va 1 in tblicher Weise mit i bezeichnet ist. Fiihren wir 
(5.7) in (5.6) ein, so erhalten wir 


KAP te) = Xm (i Bp 2) = X_* (Bat), (5.8) 
wobei 
See (Bn x) — 7S Co} (Bn a) ta es Sin (Brn al (5.9) 
Analog gilt mit x, = y, wenn wir 
1 
we = Wee = ato 5 Pn2 (5.10) 
setzen 
Xa = Y,= cit cos (a, y) + ime sin (a, y) (5.11) 
und fiir n>n, 
Jee i 2— ja 5.12 
Dae ato Pr = UOn (5.12) 
wie oben 
Y, (an Y) = Yn (i dn y) = Yn* (Gn y) (5.13) 
Yn* (Gn Y) = O*5,n C0} (Gn Y) + b*o,n Sit (An Y). (5.14) 


3. Zylinderkoordinaten (drehsymmetrisches System) 
Hier ordnen wir 
der Veranderlichen x, =r die Eigenwerte 9,,, = qn, 
der Veranderlichen x, = z die Eigenwerte 92,n = yn ZU. 


Die Differentialgleichungen nehmen dann mit (2.16) bzw. (2.15) die Formen 


d2 Xiv 1 dXiy 


ie ee + Piiv Xin = 0 (5.15) 


bzw. (5.1) an; ihre (2.5) entsprechenden Lésungen lauten 
Xin = Rn = Ayn Fo (Gn) + Oy,n No (Gn) (5.16) 
X 5m = Lm = Ge_m COS. (Ym 2) Og, m8 (Yn 2), (5.17) 


wobei unter J, und N, die Besselschen Funktionen erster bzw. zweiter Art nullter 
Ordnung zu verstehen sind; die letzteren werden bekanntlich auch als Neumannsche 
Funktionen bezeichnet. 

Da die z-abhangigen Losungen bis auf die Bezeichnung mit den fiir den ebenen 
Fall angegebenen Funktionen iibereinstimmen, beschranken wir uns hier auf die An- 
gabe der aus (5.16) folgenden, fiir den drehsymmetrischen Zylinder charakteristischen 
Funktionen. 

Mit den Koordinaten z;r der beiden Mantelflachen und /, = r,, 1, = 7, demnach 
ta > 1 > 1; folgt aus (5.16) die Bestimmungsgleichung fiir die Kigenwerte unter Beriick- 
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sichtigung der bekannten fiir Zylinderfunktionen p-ter Ordnung giiltigen Differentia- 
tionsregeln 


Bo (x) = —By +1(2) + © Bp (x), (5.18) 


hier also mit p = 0 


Bo (2) = —8: (x) | (5.19) 
fiir den allgemeinsten Fall: 


[— @eJ1 (Gn %a) +h, Fo (Pn %a)] [+ Gn Ns (Qn 7i) + hy No (Gn 7i)] — 
— [+ gad (Gn 01) + hy To (Gn i) | [—Gn Ni (Yn Pa) + h,, No (Pn %a)] = 9. (5.20) 


Fiir die Lisung X;, der (3.29) zugeordneten homogenen Differentialgleichung er- 
halten wir, wobei x, =r gesetzt wird und damit 9; = yn, also nach (3.28) 


1 
Cae = Bar am ato aes Wn (5.21) 
in (3.33) einzufiihren ist 


Xm (Bn?) = Rn (Bm 1) = 41,m To (Bm) + 814m No (Bm 7). (5.22) 


Bei wachsenden Eigenwerten y,, gilt wieder fiir m >m, 


/ 


Sif il os 
Bm = 1) i — Pm = 1 Bm (5.23) 
Hinfithrung von (5.23) in (5.22) liefert 
Rn (Bn?) = Rm (6 Bt) = Bn* (Br ?)- (5.24) 


In diesen Fallen gehen wir wieder auf die fiir die praktische Berechnung zweck- 
maBigere Darstellung 


F* (Ba 1) a* a. by (But) Oh Kee eee) (5.25) 


iiber, wobei J, und K, die modifizierten Besselschen Funktionen nullter Ordnung be- 
deuten. Allgemein besteht zwischen den modifizierten Besselfunktionen p-ter Ord- 
nung und den gewohnlichen der Zusammenhang 


TA) te? doer) (5.26) 
und 
K, (x) = 3+? H,® (ix), (5.27) 
wobei unter H,) die durch 
Hy (x) = J, (x) + 7 Ny (2) (5.28) 


definierten Hankelschen oder Besselschen Funktionen dritter Art zu verstehen sind. 


4. Bevor wir diesen Abschnitt beschlieBen, wollen wir noch auf einen Umstand hin- 
weisen, der fiir die praktische Berechnung gegebenenfalls von Bedeutung sein kann. 


In (3.12) steht der Ausdruck 
1 
Ss ti (5.28) 


im Nenner, sein Verschwinden 148t also den mit ihm verbundenen Term ins Unendliche 
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anwachsen (Resonanz); ahnliches tritt auch beim Verschwinden eines der durch (3.28) 
definierten Werte a;, ein. 

Da die GréBen a und ¢, empirisch bestimmt sind, hat es kaum einen Sinn, in diesem 
Falle umstandliche Grenziibergange vorzunehmen. Die EKigenwerte g,,, und g.,, sind 
unabhangig von fy. Es ist daher immer méglich, ohne die ,, Genauigkeit“‘ der Ergebnisse 
zu beeintrachtigen, den Wert ¢, innerhalb der in Betracht kommenden Grenzen so zu 


es : ée 1 1 ; 
wahlen, daB keiner der Ausdriicke ale Y1.n — P2,m und ah ~ Viv verschwindet. 


5. Wie schon eingangs betont, gelten die vorliegenden Untersuchungen vor allem 
den Aufheizungsvorgangen im abbindenden Beton und den von ihnen hervorgerufenen 
Spannungszustanden. 

Die nachstehend angegebenen Grenzwerte?® fiir die in die Rechnung eingehenden 
Konstanten beziehen sich daher nur auf Beton. 


Warmeleitfahigkeit “2 .............00. 1,0 bis 2,2 Cal/mh °C 
Temperaturleitzahl —@ _. 26.6.5 ee nee es 0,003 bis 0,004 m2/h (a = 4) 
Spezitische Warmer ...:.:,... 5.0.0.0. 0,21 bis 0,27 Cal/kg °C 
Raumgewicht Visita. t wanatisales te 2300 bis 2500 kg/m? 
Warmeiibergangszahla ................ 12 bis 40 Cal/m?h °C 
Warmeentwicklung We ............04.. 100 bis 300 Cal/m*h 
Warmeausdehnungs- 

koeffizient re ae en ore Pere 0,9: 10-5 bis 1,2- 10-5 
Abklingungskoeffizient ¢ aaare te artesfo rues 0,005 bis 0,020 h-} 


VI. Spannungsfeld: Partikularlésung 


1. Nach Festlegung des Temperaturfeldes sollen im folgenden die Spannungen be- 
rechnet werden, wobei wir uns der von E. Melan und H. Parkus entwickelten all- 
gemeinen Verfahren’ bedienen, die vor allem durch die Verwendung des thermisch- 
elastischen Verschiebungspotentials gekennzeichnet sind. 


2. Zur Erklarung dieses Begriffes schreiben wir zunachst die auf thermische Wir- 
kungen erweiterten Grundgleichungen fiir den Verschiebungsvektor u; an: 
Tb Pde 9 (211 ty) b8 


o 
a aS ] — 2» Cra a OF (6.1) 


hierin bedeuten e = >’ du;|dx; die kubische Dilatation, » die Poissonsche Zahl — 
Verwechslungen mit dem oben verwendeten Summationsindex » sind kaum zu befiirch- 
ten — und a7 den Temperaturausdehnungskoeffizienten. Partikularlésungen der Grund- 
gleichungen sind nun durch den Gradienten 

ee (6.2) 


+ OR; 


einer Funktion ® gegeben, die daher als thermisch-elastisches Verschiebungspotential 


bezeichnet wird. 
Die Einfiihrung von (6.2) in (6.1) liefert die Differentialgleichung 


Ab =** "ard, (6.3) 


Caan O eae ING CXL 
Aras ane. 


Ingenicur-Archiy XIIT/2 6 
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aus der ein partikulares Integral von ® zu bestimmen ist. Mit ® gilt ferner fiir den 
dem Verschiebungsfeld wu; zugeordneten Spannungszustand 


Sie = 26 (Se Adit), (6.4) 


OX; OX, 
wobei G den Schubmodul und 


eee: 
b= 44 fiir i + & 988 


Mit (6.3) kann (6.4) auch in der Form 


cay we a2 D 1+ ¥y ne ‘ 
Gin = 24 (SS — 7 ar 0 bi) (6.44) 


Kronecker-Symbol bezeichnen. 


geschrieben werden. Die aus dem Verschiebungspotential gemaf (5.4) berechneten 
Spannungen erfiillen aber die Randbedingungen — wir werden hier spannungsfreie 
Rander oder Oberflachen verlangen — i. a. nicht. Es ist daher ein zweiter Spannungs- 


zustand o;, zu ermitteln, mit dem die aus (6.4) folgenden Oberflachenspannungen getilgt 


werden kénnen. Die endgiiltigen Spannungen sind dann durch o;, = o;, + oi% gegeben. 
Wir haben also vor allem das thermisch-elastische Verschiebungspotential ® aus (6.3) zu 
bestimmen. Eine Partikularlésung kann unmittelbar aus der Differentialgleichung (1.1) 
der Warmeleitung gewonnen werden, die wir hiezu auf die Form 


o=a (fa Odt + Whe d | (6.5) 


bringen; da es sich hier bloB um eine andere Darstellung des Temperaturfeldes handelt 
und nicht etwa um ein Partikularintegral der Differentialgleichung (1.1), muB (6.5) den 
Randbedingungen (1.2) und (1.3) geniigen. Schreiben wir auf der rechten Seite fiir 


d=94+4, 
wobei # wie oben das allgemeine Integral der homogenen Differentialgleichung (2.1) 


und #@ ein Partikularintegral der inhomogenen Differentialgleichung (1.1) bedeuten, so 
k6onnen wir (6.5) unter Verwendung von (3.6) nach den einzelnen Zeitfunktionen auf- 
spalten: 


t 
ga al (ade 4 )emars (sda). (6.6) 
Wird die im Storungsglied a = row (%1,%2) auftretende Funktion w als Laplace- 


Ableitung einer noch zu bestimmenden Funktion Y aufgefaBt: 
Ae eee (6.7) 


dann folgt aus (6.6) bei Vertauschung der Reihenfolge von Integration und Laplace- 
Ableitung und unter Beriicksichtigung der Zeitunabhangigkeit von 0* und W: 


t t 
= Ww ~— 
eel (o* +2 v)fema+ af oa), (6.8) 
Ausfiihrung der Quadratur im ersten Glied liefert weiter 
$s gr (0 +2 v) en tlt 4 aA|ddt, (6.9) 


und wir erhalten nach Hinfiihrung von (6.9) in (6.3) und Durchfihrung der zur Laplace- 
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Ableitung inversen Operation, d. h. Weglassung des auf beiden Seiten stehenden 
Laplace-Operators ein Partikularintegral von (6.3) 


t 
] oa af 
ace (a eos v) o£ af Dat), (6.10) 


in dem die Funktion Y noch unbekannt ist. 


Bevor wir aber Y berechnen, soll zunichst das zweite, in (6.10) innerhalb der eckigen 
Klammer stehende Glied bestimmt werden, das den EinfluB des Integrals der homogenen 
Differentialgleichung (2.1) auf das Verschiebungspotential angibt. 


Aus (2.2) erhalt man in entwickelter Form 


— co co K 
ee Se wy X. PTE a Ses, © 
| = Pint Pam Th M 1) X2,m(P2, m t2) € , (6.11) 


woraus nach Hinfiihrung von K,,, gemaB (4.5) schlieBlich 


t 


1 P71,n + P2,m 
pe 2 > ; 
fe P*1.n ee te} (6.12) 


n=1m=1 


folgt. 


Zur Berechnung der Funktion Y haben wir nun ein partikulares Integral der Diffe- 
rentialgleichung (6.7) aufzusuchen. 


Fiir den allgemeinsten Fall einer nicht in Produktform darstellbaren Stérungsfunk- 
tion w setzen wir w gema&8- (3.1) in einer nach den Eigenfunktionen X,,,,-X,_, ent- 
wickelten Doppelreihe an und erhalten damit die Differentialgleichungen 


A= > pif We mame g Sy bp. Cab. Pier (6.13) 


— — 
n=1m=1 aS mn 


die fiir jedes Reihenglied erfiillt sein miissen. Setzen wir auch VY in einer nach X,,,° Xo, m 
entwickelten Doppelreihe an, so folgt unter Beachtung von (3.10) aus (6.13) 


gees — t Ope Mase D. = (6.14) 
mit r 
Cum 
Bam = — P71,n + P22,m 
also 
rae = past aide ny Xo,m- (6.15) 


P71,n + Y?2,m 


Zusammenfassung von ie Ym mit dem durch oe gemaB (3.7) mit (3.12) gegebenen 


Partikularintegral von (1.1) und Multiplikation mit — at, liefert schlieBlich fiir den 
ersten in der eckigen Klammer stehenden Term von (6.10): 


= su (Be4 Bepem = SSM = Susuiey 
A fae 4 ae tT 3, 1 — Pam) (p1 n+ 92, m) 6.16 
ato 7 : ( . ) 


6* 
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Setzen wir (6.16) und (6.12) in (6.10) ein, so folgt fiir das thermisch-elastische Ver- 
schiebungspotential der Ausdruck: 


oo C 
Dp — rer laws v 1 = tie Xe ment — 
f einai (Ge Pn — Pm) (21,5 + 9%,m) 
an Sy y" 2 ou ats a 1 es Xin Xe, me aPrn + F'aom) t 4 
iors Sess P*1,n Pr2,m (Ga — phn — am) (e212 + 9%2.m) (6.17) 


Kann aber der mit der Zeitfunktion e~/” verbundene Teil des Temperaturfeldes 
durch einfache Reihen gema8 (3.20) durch (3.42) oder (3.43) dargestellt werden, dann 
wird w in (6.7) zweckmaBig auch in der nach einer einfachen Reihe entwickelten Form 
(3.20) einzufiihren sein. 

Mit 


bas 


AY = w, (2%) Aiv Xiv (Qiv 23) (6.18) 


v=1 


I 


und dem (3.22) entsprechenden Ansatz 


—— reXee (Viv Xi) Wy (Xx) (6.19) 


v=1 


erhalten wir aus (6.7) die Differentialgleichung 


A (Xi Vey) == 40% Aiy D, a = OF (6.20) 
Fiir die Laplace-Ableitung des Produktes X;, Y%, folgt unter Beachtung von (2.4): 
A (Xiy Ves) oe? (— Py ees + A bares y. Cee (6.21) 


Wird (6.21) in (6.20) eingefiihrt, so ergeben sich nach Kiirzen von X; und einfacher 
Umformung die gewohnlichen Differentialgleichungen: 


A es = Priv ee = As WwW = 0. (6.22) 


Da fiir die Funktionen ¥%, keine Randbedingungen vorgeschrieben sind, geniigt es, 
irgendein die Differentialgleichung (6.22) erfiillendes Partikularintegral aufzusuchen. 
Wir schreiben diese Lésungen in der Form 


Py = Aiy Pry (6.23) 
an und erhalten mit (6.19) 
Y= a Aiv Xiv (Piv Li) Dev - (6.24) 


Diese Funktion wird nun i. a. mit der aus (6.15) berechneten nicht iiberein- 
stimmen. Da aber sowohl (6.23) wie auch (6.15) Partikularintegrale von (6.7) sind, 
kénnen sie sich nur um eine Funktion mit verschwindender Laplace-Ableitung unter- 
scheiden, die wir als Lésung der (6.20) zugeordneten homogenen Differentialgleichung 
erhalten; sie beeinflu8t den Spannungszustand nicht und bleibt daher im weiteren 
auBer acht. 

Kinfiihrung von (6.24) in den ersten Klammerterm von (6.10) ergibt unter Ver- 
wendung von (3.40) 


ar ato e— [to (9* — a v) = ~Aty Y e-t Bs A* Es (dix Xe — fe (xx) a Pk (2x)) . (6.25) 
A A 1 


\ y= 
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Da keine der Veranderlichen ausgezeichnet ist, kénnen wir in (6.25) die Indizes i 
mit k bzw. v mit 4 vertauschen. Die Zusammenfassung von (6.25) mit (6.11) ergibt 
fiir das thermisch-elastische Verschiebungspotential die beiden Ausdriicke: 


lip Wo = oo 
P= 5;— sar | Oty 7 tite DAY kate (diy Nek (2a) = 
n=-z 
Fe - See X4 OG: x ea (P'1sn + PP2sm)t \ (6 26) 
Fal male be + eae z at i A 
und 
l+y» { Wo cee 
@ = ke | ato 7] e—t/to Ds es P. Ce (de,mX1,m SE WO cy (21)) £23 
m=1 
Aue y Si elem xX Xe —a(@*isn + P2:m) t ¢ 
a 2 etme 1,n 02, mé 3 (6.27) 


in denen die Funktionen X,,,,, X,,, und p,,,, P;,m durch (3.32) und (6.23) gegeben sind. 

Bei der praktischen Berechnung wird man von Fall zu Fall entscheiden miissen, 
ob dem thermisch-elastischen Verschiebungspotential ® der Ausdruck (6.17) oder eine 
der beiden Darstellungen (6.26) bzw. (6.27) zugrunde gelegt werden soll. Den Nachteilen 
der Doppelreihen steht im ersten Fall der Vorteil gegeniiber, da8 im thermisch-elasti- 
schen Verschiebungspotential nur die Ortsfunktionen X, ,»(Q1,n 2%) X2,m (P2,m L2) auf- 
treten, die Spannungsberechnung daher nur fiir diese durchgefiihrt werden mu. Im 
zweiten Fall dagegen kommen, wie aus (6.25) zu ersehen ist, zu diesen noch die 
Funktionen 


Ny Oig We) vbw (Civ Lb )on 4 aiv Gin Vi RALE) 
und Xiy (Viv Li) Pe (Xk) 


hinzu, aus denen die Spannungen o:, nach (6.4) berechnet werden miissen. Jeder dieser 
Funktionen entspricht ferner ein besonderer Randspannungszustand 6;4, dessen Uber- 
lagerung mit o;, die endgiiltigen Spannungen liefert. 


VII. Spannungsfeld: Homogene Lésung 


1. Mit der Bestimmung des thermisch-elastischen Verschiebungspotentials und der 
aus diesem entspringenden Spannungen o;, soll die allgemeine Behandlung des vor- 
liegenden Problems fiir rechteckige und zylindrische Kérper abgeschlossen werden. Im 
rotationssymmetrischen Faille wird der zur Erfiillung der Randbedingungen erforder- 
liche Spannungszustand o;, bekanntlich aus einer Loveschen Verschiebungsfunktion 
erhalten; in analoger Weise konnte man fiir den ebenen Fall die Marguerresche Ver- 
schiebungsfunktion verwenden und die entsprechenden Randspannungszustande wie 
oben in einer fiir beide Systeme giiltigen Form ableiten. Wir werden hier von dieser 
Méglichkeit Abstand nehmen und im ebenen Fall, dem iiblichen Gebrauch folgend, auf 
die Airysche Spannungsfunktion zuriickgreifen. Die sinngemaBe Ubertragung der 
in der Scheibentheorie zur Lésung des Randwertproblems angewandten Verfahren auf 
rotationssymmetrische Falle kann sich dagegen in manchen Fallen als zweckmaBbig 
erweisen. 

2. Die hier fiir das ebene Problem verwendete Airysche Spannungsfunktion F ist 
durch die Bipotentialgleichung 

AAF=0 (7.1) 
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und durch ihre Randwerte, also 


Gi ae 03, = —Gii (7.2) 


und 


+k Oi = —Ojk (7.3) 


bestimmt, wobei die Spannungen mit F allgemein durch 


cin = — ( 5 — 5 AF) (7.4) 


OX; OX’ 


verkniipft sind. ; 
Dieses Problem hat K. Girkmann fiir den hier in Betracht kommenden Fall einer 


beliebigen Randnormal- oder Schubkraften ausgesetzten Rechteckscheibe gelést’. 
Temperaturspannungen in einfach zusammenhangenden Bereichen kann man 
jedoch auch, wie der Verfasser in einer friiheren Arbeit gezeigt hat, auf Grund der 
Plattenanalogie berechnen. Die auf Temperaturwirkungen erweiterte Differential- 
gleichung der Airyschen Spannungsfunktion (7.1) lautet bekanntlich: 


AAQ=—7*"arAd (7.5) 
und stimmt in dieser Form mit der Differentialgleichung der Biegeflache einer pro- 
portional zu A # belasteten diinnen Platte tiberein. Ist diese Platte langs ihrer Beran- 
dung s unverdrehbar gelagert, dann mu8 die nunmehr als Biegeflache der querbela- 
steten Platte gedeutete Funktion 2 bei einfach zusammenhangenden Bereichen in s 
den Bedingungen 


C= o (7.6) 
und 

AQ 

free 0 (7.7) 


geniigen ; die gleichen Bedingungen gelten aber auch fiir die einem geometrisch ahnlichen 
Querschnitt zugeordnete Airysche Spannungsfunktion bei verschwindenden Rand- 
spannungen o,, und o,,. Wir kénnen daher auf die Berechnung des thermisch-elastischen 
Verschiebungspotentials im ebenen Falle iiberhaupt verzichten und die Spannungen 
unmittelbar durch Umrechnung aus den Momenten einer proportional zu 

1 0d Wo 


belasteten eingespannten Platte gewinnen. Da fiir derartige Platten die EinfluBfelder 
der maBgebenden Momente bereits berechnet vorliegen, bedeutet dieses Verfahren 
eine erhebliche Verringerung des Rechenaufwandes. Beziiglich weiterer Einzelheiten 
iiber Anwendung der Plattenanalogie auf Temperaturspannungsprobleme sei auf die 
erwahnte Arbeit des Verfassers® verwiesen. 


3. Uber die von quellfreien instationaren rotationssymmetrischen Temperatur- 
feldern hervorgerufenen Spannungszustande hat R. Trostel in einer ausfiihrlichen 
Arbeit® berichtet. Wir werden in der Folge an die Ergebnisse dieser Untersuchungen, 
soweit sie fiir unsere Zwecke in Betracht kommen, ankniipfen. 


7 K. Girkmann: Flachentragwerke, 4. Aufl., S. 89ff. Wien:. 1956. 
8 S. FuBnote 6 auf S. 59. 
® §. FuBnote 5 auf S. 59. 
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Wie schon oben erwahnt, werden die Randspannungszustande og; o;x hier aus Loveschen 
Verschiebungsfunktionen L abgeleitet, die der Bipotentialgleichung 


ADEA (7.9) 


geniigen miissen. Der einer Loveschen Verschiebungsfunktion zugeordnete Spannungs- 
zustand o;, ist mit ihr durch die nachstehenden Beziehungen verkniipft: 


a =i (r4t- SS) (7.10) 
Foe = jects Oe (> Agha . fo set 
a... le NAL — FT| (7.12) 
oe la AL—T), (7.13) 


Orr + Sgy 1 Oz2= 7-9, Ta Ne +y)AL. 


Den belanglosen Faktor ;— z, werden wir hiebei im weiteren weglassen. 


Der in der zitierten Arbeit von R. Trostel eingeschlagene Weg besteht nun darin, 
zuerst die aus dem thermisch-elastischen Verschiebungspotential abgeleiteten, an den 
beiden Mantelflachen auftretenden Spannungen o,, und o,, durch passend gewihlte 
Linearkombinationen von Loveschen Spannungszusténden zu kompensieren. Wir 
schreiben zu diesem Zweck die zwischen dem thermisch-elastischen Verschiebungspo- 
tential und den Spannungen bestehenden Beziehungen (6.4) in Zylinderkoordinaten 
an: 


es 2G(22 cenitdl ©) (7.14) 

i 26(— eee Ao) (7.15) 

nye 2G(S2 naif, °) (7.16) 
a2 D 

ip yieS 2G. (ie) 


Unter Beachtung von (6.3), (7.14), (7.17) sowie der Differentialgleichung (2.4) folgt, 
da8 die z-abhangigen Terme in den Ausdriicken fiir o,, als nullte und zweite Ableitun- 
gen mit den gema8 (5.17) definierten Funktionen Z(z) tibereinstimmen, wahrend in 
o,, die ersten Ableitungen Z'(z) auftreten. Die aus der Loveschen Verschiebungsfunk- 
tion folgenden Spannungen 0,, bzw. o,, sind dagegen durch ungerade bzw. gerade Ab- 
leitungen mit den in L auftretenden z-abhangigen Funktionen verkniipft. Zur Erzwin- 
gung des identischen Verschwindens der Spannungen Or + o;, und o,, + o,, langs der 
beiden Mantelflachen miissen die z-abhangigen Terme in der Loveschen Verschiebungs- 
funktion also gleich Z’(z) gesetzt werden. Um einen allgemeinen Uberblick iiber die 
hier in Betracht kommenden Lésungen von (7.9) zu erhalten, treffen wir den Kreis- 
und Hyperbelfunktionen erfassenden Ansatz 


L = L* (r) e®, (7.18) 
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aus dem unter Beachtung von (2.16) vorerst 


AAL=Al(SF Beatie wt L*)ee=| = 0 (7.19) 


or? y or 


folgt. 
Der in (7.19) innerhalb der runden Klammer stehende Term entspricht der Differen- 


tialgleichung (5.15) und verschwindet daher fiir die durch (5.16) definierte Funktion R. 
Fiihren wir ferner in (7.19) die Funktion 


oR 
ee (7.20) 


ein, so ergibt sich — wie man leicht nachrechnen kann — 


a2 (r. BR: 1 d(r.R = ae 
eS ee OCR) ee & (7.21) 
und damit aus (7.19) 
AAL=A (Re) =0, (722)aa 


so daB auch r R: eine Lésung der aus (7.18) folgenden, gewohnlichen Differentialglei- 
chung 4. Ordnung ist. Setzen wir noch fiir R: ein, indem wir (5.16) gemaB (5.19) nach r 
differenzieren, dann erhalten wir die allgemeine Losung von (7.19): 


Ly* = P*y,y J q (yt) + P*o,4 No (ovr) + Q*1,u 0 Jy (yt) + Q* 2,00 Ny (@yr) (7:28) 


mit den Konstanten Po, Baa und OF ew: QF as 
_ Fir die Z(z) gemab (5.17) und damit auch Z’(z) aufbauenden Kreisfunktionen mu 


w rein imaginar, also w = 7 sein. Ersetzen wir w durch die in (5.17) auftretenden 
EKigenwerte y,,, dann folgt entsprechend (5.25) sowie der unter Beriicksichtigung von 
(5.26), (5.27) und (5.28) angewandten Differentiationsregel (5.19): 


L*® = P¥ 1m Lg (mt) + P*o,m Ko (Wnt) + Q* 1m Ly (Ym) + Q*o,m% Ky (ymr). (7.24) 


Unter Beachtung der zwischen Exponentialfunktionen und Hyperbel- bzw. Kreis- 
funktionen bestehenden Zusammenhange kénnen wir die Funktionen L,* bzw. Ln* mit 


beliebigen Linearkombinationen der Argumente w, z bzw. yz verbinden und erhalten 
so aus (7.23) die Lovesche Verschiebungsfunktion : 


LL, = eres. (wy 7) ae jeepville (wy 1) Qt Jy (wy 1) +O. ut Na (wy r)] ° 
- (A, Gof (wy z) + By Sin (wy 7)] (7.25) 


und aus (7.24) speziell die fiir die Kompensation der Mantelspannungen o,, und a,, 
in Betracht kommende Verschiebungsfunktion: 


Ln = [Pim Lo (Ym?) + Pam Ko (mt) + Qimt I, (Ym 1) + Qomt Ky (pm 1)] Z (pm). 
(7.26) 
Bildet man nun aus (7.26) die Spannungen g,, und 3g, z, 80 erhalt man mit den an 
der 4uBeren und inneren Mantelfliche zu erfiillenden Randbedingungen a,, + a,, = 0 


und o,, +4,,=9 vier Gleichungen, aus denen die Konstanten Pals ae oune 
Q2,m Zu bestimmen sind. 
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In einem zweiten__ Rechnungsgang werden dann die an den Stirnflachen verbliebenen 
Spannungen o,, + o,, und o,, + o,, zum Verschwinden gebracht. 


Zufolge der Drehsymmetrie des Spannungsfeldes stehen die Schubspannungen in 
sich im Gleichgewicht. Aber auch die Spannungen a,,, liefern — wie sich leicht zeigen 
1aBt — keine Resultierende. Nach dem St. Venantschen Prinzip kann also das gemak 
dem friiheren berechnete Spannungsfeld, das den Randbedingungen im strengen Sinne 
nur langs des Zylindermantels geniigt, fiir Zylinder mit groBem Héhen-Durchmesser- 
Verhaltnis als endgiiltig angesehen werden. Ist dieses Verhaltnis klein, dann ist auch 
der Kinflu8 der an den Stirnflachen verbliebenen Restspannungen zu beriicksichtigen. 
Fir diesen Fall hat R. Trostel ein Verfahren entwickelt, auf das hier nur verwiesen 
werden soll. SchlieBlich kénnte aber auch das oben erwihnte von K. Girkmann zur 
Lésung des Randwertproblems fiir die Rechteckscheibe entwickelte Verfahren sinn- 
gemaB auf den rotationssymmetrischen Zylinder iibertragen werden. 


4. AbschlieBend sei noch betont, daB zur wirklichkeitstreuen Erfassung der beim 
Abbinden des Betons auftretenden Spannungszustinde auch die Zeitabhangigkeit der 
Elastizitatskonstanten zu beriicksichtigen ware. 


Wahrend die Querzahl des nunmehr fiir Massenbetonbauwerke ausschlieBlich 
verwendeten Riittelbetons bereits nach etwa einer Woche als konstant anzusehen ist 
(die Unterschiede zwischen den Werten, die » bei Zug und Druck annimmt, bleiben hier 
unberiicksichtigt), erreicht der Elastizitats- bzw. Schubmodul seinen Endwert erst 
nach einem erheblich groBerem Zeitraum™'. Die Zeitabhangigkeit kann in erster An- 
naherung etwa durch ein Exponentialgesetz der Form 


E (t) = E, (1 — et) (7.27) 


beschrieben werden, in dem die empirisch zu bestimmenden Konstanten L, den End- 
wert, dem sich der H-Modul asymptotisch nahert, und tz eine den Grad des zeitlichen 
Anstieges von H kennzeichnende Zeitkonstante bedeuten. 


Das gleiche gilt aber auch vom Temperaturausdehnungskoeffizienten a7, dessen 
Veranderlichkeit von wesentlichem Einflu8 auf die hier untersuchten Spannungs- 
zustande ist. 


Da die elastischen Grundgleichungen (6.1), auf die die vorstehend verwendeten 
Beziehungen aufgebaut wurden, bei den in Betracht kommenden langsamen Temperatur- 
anderungen ihre Giiltigkeit auch fiir zeitabhangige Elastizitatskonstanten behalten, 
kann in den hier entwickelten Spannungsformeln das Produkt a7G analog zu (7.27) 
durch : 


ap Eo 
(1 + ») 


arG = 5) (1 e~tIta B) (7.28) 
ersetzt werden, wobei die Werte ar EZ, und t,, von Fall zu Fall aus Versuchen zu be- 
stimmen sind. 

Fiir iiberschlagige Rechnungen geniigt es, arG mit einem dem Alter des Betons 
entsprechenden Mittelwert einzufihren. 


10 Hirohiko Yoshida: Uber das elastische Verhalten von Beton. Berlin. 1930. 
11 Sixth International Congress on Large Dams, 1958, Communication 24. 
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Drehschwingungen von Zahnradgetrieben 
Von F. Bohm, Stuttgart 
Mit 13 Textabbildungen 


Zusammenfassung: In der grundlegenden Arbeit von J. Zeman?’ wird das Zahnradgetriebe 
als linearer Zweimassenschwinger behandelt. Die lineare Schwingungsgleichung mit konstanten 
Koeffizienten beschreibt jedoch das Drehschwingungsverhalten eines stationar laufenden Zahnrad- 
getriebes nicht ausreichend genau. Praktisch hat sich ein Einflu8 des ibertragenen Drehmomentes 
auf die dynamischen Zahnkrafte gezeigt. Als zusatzliche Eigenschaft wird daher die periodisch 
verinderliche Federzahl, welche sich durch den Zahneingriff ergibt,in die Rechnung hineingenommen. 
Die nichtlinearen Eigenschaften des Getriebes, das Zahnspiel sowie die statische Unbestimmtheit 
der Kraftibertragung, wenn zwei Zahnpaare gleichzeitig kimmen, werden vorerst unterdriickt. 
Dies erscheint fiir gut bearbeitete Zahnrider durchaus plausibel. Die sich ergebende Bewegungs- 
gleichung ist eine Hillsche Differentialgleichung. Es gelingt dann mittels der zugeordneten 
Volterraschen Integralgleichung, eine beliebig genaue Naherungslosung fir die Schwingungs- 
beanspruchung zu finden. AnschlieBend wird der Einflu8 der Nichtlinearitaten bei kleinen Werten 
des Flankenspieles und gleichmaBigem Tragbild auf diese Naherungslosung naher untersucht und 
als im allgemeinen nicht erheblich festgestellt, wodurch sich ein relativ einfaches Rechenverfahren 


fiir die dynamische Beanspruchung ergibt. 


I. Einleitung 


In neuerer Zeit hat sich die Ansicht durchgesetzt, da die dynamischen Zusatz- 
krafte fiir die Lebensdauer eines Zahnradgetriebes, besonders bei hoheren Drehzahlen, 
maBgebend sind. Das schwingungsfahige System besteht aus den Drehmassen der 
beiden Zahnrader und den deformationsfahigen Zahnen, die gerade im Eingriff sind. 
Das zu iibertragende Leistungsmoment erzeugt eine gewisse Vorspannung dieser Zahne, 
die bei Schwingungen im allgemeinen positiv bleiben wird, aber nicht mehr konstant 
sein kann. Sinkt die Vorspannung durch die dynamischen Krafte auf Null, so wird 
das Zahnflankenspiel durch- 
laufen und der Zahn schlagt 
an der gegeniiberliegenden 
Flanke an (s. Abb. 1, welche 
die Verhaltnisse beziiglich 
der Riickstellkraft zeigt). 
Abb. 1. Verlauf der Rickstellkraft R als Funktionder Fir negative Werte von x 

Zeit und der Zahnverformung x hat man eine nichtlineare 


S ... Zahnspiel, «= 1+ " a ... Uberdeckungsgrad Schwingung zu _ erwarten. 
uy 2 
C ... Federzahl 


II. Aufstellung der Bewegungsgleichung 


Die Federzahl des Systems ist davon abhangig, wieviel Zahne gerade im Eingriff 
sind. Nach + ? ist sie praktisch als linear abhangig von der Zusammendriickung 
anzunehmen, da der nichtlineare Teil, die Hertzsche Pressung, sehr rasch ansteigt 
und dadurch nur geringen Einflu8 nimmt. Bei Geradzahnradern ergibt sich eine 
wesentlich groBere Schwankung der Federzahl als bei Schragzihnen. Diese Schwankung 
erfolgt ziemlich sprunghaft? und ist von der GroBe + 30 bis 35%. Ist der Uberdeckungs- 
grad des Getriebes z. B. ¢ = 1,5, so sind wahrend einem Drittel der Hingriffsdauer 
ein Zahnpaar und wihrend zwei Drittel zwei Zahnpaare im Eingriff. Die Federzahl 
1aBt sich gut durch eine Fourier-Reihe, abhangig vom Drehwinkel y, darstellen. 


1 J. Zeman: Dynamische Zusatzkrifte in Zahnradgetrieben. VDI-Z 99 (1957). 


> F. Karas: Elastische Forminderung und Lastverteilung beim D leingri i 
ae Hoa sat g beim Doppeleingriff gerader Stirn- 


| 
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Sind zwei Zahnpaare im Eingriff, so liegt ein statisch unbestimmtes System vor. 
Ks soll aber vorerst angenommen werden, da8 beide Zahnpaare gleichmaBig zum Tragen 
kommen. Der Ubergang von C,,,x auf Cyin ist auch bei stationarer mittlerer Drehzahl 
des Getriebes von der tiberlagerten Schwingung abhangig. Die Schwingungsamplituden 
sind jedoch sehr klein gegeniiber der Zahnteilung, so daB der Einflu8 auf C unwesentlich 
sein diirfte. 

Fir die Bestimmung der Erregung des Schwingers ist nur der Sammelfehler eines 
Zahnradpaares von Interesse. Er ist die Differenz zwischen dem wirklich zuriick- 
gelegten Drehwinkel des Rades 2 gegentiber dem Drehwinkel des Rades 2 bei fehler- 
freien Verzahnungsverhaltnissen und verschwindend kleinem Zahndruck. Bei starkem 
Zahndruck verandert sich dieser Fehler infolge Einebnung des Oberflichengebirges 
der Zahnflanken geringfiigig, auBerdem entsteht ein stiickweise konstanter Fehler- 
anteil durch elastische Deformation der Zahne. Auch der Uberdeckungsgrad wird 
dadurch, wie schon oben erwahnt, in gewissen Grenzen schwanken. Wenn man aber 
von extremen Fehlern® absieht, kann man den Uberdeckungsgrad naherungsweise 
konstant annehmen. Wir werden uns also darauf beschrinken, nur den WaAlzfehler 
bei kleinem Zahndruck zu beriicksichtigen, und voraussetzen, daB dieser Fehler sich 
bei grofem Zahndruck nur insofern andert, als ein stiickweise konstanter Anteil 
infolge der rechteckfo6rmig veranderlichen Elastizitaét hinzukommt. 

Wie in der Arbeit von J. Zeman soll auch hier der WAlzfehler mit der Periode der 
Zahnteilung angesetzt werden, wodurch sich eine periodische Stérkraft ergibt; als 
weitere Storkraft sei das eingepragte konstante Drehmoment M, angesetzt. 


Fir die in Abb. 2 gezeichnete Anordnung gelten folgende Bezeichnungen: 


eel AN Oe ae Massentragheitsmomente 
ee Mgt ck ges eingepragte Drehmomente 
err at Sn ok Ubersetzung 
iy AN EE RR OR Zusammendriickung der Zahnflanken 
in Zentimeter 
Nyt a ae Flankenfehler gemessen im Bogen- 
ma, positiv wenn Rad 2 voreilt FEES 
ie tik « Beenie toss Federzahl "aS 
aa A Se a Zahnezahl des Rades 1, 2 
028 AE OR lee eM Winkelgeschwindigkeit des Rades 1, 2 
Ra Ra eee Uberdeckungsgrad 
Deo. 4 isitay Zahneingriffsdauer 
1 CTS oe Se Periode der Federzahlschwankung 
Se a Ae A A Kingriffsfrequenz Abb. 2. Schema eines Zahnrad- 
| ao ees a Zeit fiir eine Umdrehung des Rades 1,2 getriebes 
NS GA ae Dampfungszahl 


Die Federzahl des ideellen, statisch bestimmten Schwingers ohne Spiel ist abhangig 
von gy, und x. Da x aber sehr klein ist gegentiber der Eingriffsteilung, soll von dieser 
Abhangigkeit abgesehen werden, so daf8 man ansetzen kann: 


C=C (gy, ) = C (pr) 
Wegen 9, = — t+ Aq,, wobei 4g, von gleicher GréBenordnung wie x ist, 
1 


kann schlieBlich geschrieben werden: 


Opldat) a= Go dae ne dort) (1) 


3 —. Niemann und H. Rettig: Dynamische Zahnkrafte. VDI-Z 99 (1957). 
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Mit w, wurde die Eingriffsfrequenz bezeichnet, welche mit der Kreisfrequenz der Zahn- 
raddrehung iiber w, = 2, ©; = 2, @, zusammenhangt. Aus Abb. 3 sind Zahneingriffs- 
dauer 7'z, Uberdeckung «, Federzahl Cy) und Walzfehler fy) in Abhangigkeit von 


der Zeit t ersichtlich. Die Zahneingriffsdauer ist 7'— = ¢ T, = es. Dadurch wird 


WO, = 20, = 01% = W222 (2) 
if z i T. 


: woe ie 
at oe re AAA eae pet Ved omar ke Fiir die Drehwinkel ¢,,, gilt nach 
a le ea Abb. 2, wenn Drehungen entgegen 


ft) a B. dem Uhrzeigersinn positiv ange- 
nommen werden 


Lei t 
Abb. 3. Federzahlschwankung und periodischer Walz- P1-t + fo — ivatmrase P2 
fehler in Abhangigkeit von der Zeit , 
r=, (%t+ F2t+fo) (3) 


Die kinetische Energie des Systems ist 
T=1,9;+ 1, 92 (4) 
Fiir die potentielle Energie gilt 


2 


2 
U = —0w ys = —Cw FS (G2? + G2? + fo? + pr g2tt Wifwot+ 2wafo) (5) 


Mit der Lagrangeschen Funktion ZL = 7’ — U lauten die Bewegungsgleichungen 


d/ aL ol 
5 (Sg) — Spc (6) 
Da wir eine reine Rechtecksschwankung fiir C annehmen, wird oy = 0. Die Ab- 
* 


weichung von dem wirklichen Verlauf (s. Abb. 4) ist sehr gering, d. h. a ist klein 
gegen C. Setzt man 7’ und U ein, so ergibt sich 


LT, Gi + Cw (G1? + Got + fy ) 1? = M,|\ 

I G2 + Cw (Qrt + 2 + fo) 72? = MJ 

Der nichtlineare Term ~ oS wurde wegen der oben er- 

wahnten Kleinheit von a sowie wegen x < 1 unterdriickt. 

Die Zahndeformation 2 ist sehr klein gegeniiber der Ein- 
griffsteilung, das Quadrat daher vernachlassigbar. 


Abb. 4. Gesamtverformung 
“(statisch) und Federzahl CO 


eae te Mit Gleichung (1) ergibt sich 
i+ soe x ry = hi 
b+ oar, = 
lo=fo 
wtOn (G+ 7) ons ite + hi r 


Praktisch wird immer eine gewisse Dampfung vorhanden sein, entstanden durch 


baie ae des Olfilms, Schallabstrahlung und durch Werkstoffdampfung in den 
ahnen. 


Ea a 
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Ks ist schwierig, die Dampfung exakt vorauszuberechnen, so daS man in der 
Rechnung auf einen geschwindigkeitsproportionalen Ansatz angewiesen ist. Aus Ver- 
suchen* sind Dampfungskonstanten bestimmt worden, mit denen wir uns vorlaufig 
begniigen. 

Vervollstandigt lautet also Gleichung (8): 
1 I, TI: *. ; M,. M: 
ler) Sees sO OMe (7% aie bas fis) nie (9) 
Von F. Karas wird die Abhangigkeit der Feder- 
zahl vom Drehwinkel bzw. vom Eingriffspunkt = c/,t/ 
bestimmt und eine derartige Federzahlfunktion 
auch praktisch berechnet (Abb. 4). Man sieht, daB 
selbst die abschnittweise Konstanz nicht ganz 
zutrifft. 

Fir die rechteckige Federzahlschwankung nach 


Abb. 5 kann man nun leicht eine Fourierentwick- Abb. 5. Idealisierte Federzahl- 
lung angeben. Mit ¢ als Uberdeckungsgrad ergibt schwankung C(«,1) 
sich aus Abb. 5: 
C (wt) = a +>) a, cos nwt = 2e—3 + 2 (1) sin nx (¢ — 2)-cosna,t 
1 1 
(10) 
AbschlieBend fiihren wir noch folgende Abkiirzungen ein: 
ES OT ie Soot RL a in airs Ae a a Kigenzeit 
: Ii Te 1 ; < 
Mra = The Trcic reduzierte trage Masse 
i /eMinns. M. ITs : Shs : 
P= (7% aia 7) 5 oy Umfangskraft. (Bei stationirem Lauf ist 
Op. = VC miter /Mrea Merce DRE, Cohay. keer mittlere Eigenfrequenz 
ENO Ne CT ee AS Frequenzverhaltnis 
D= [les Ge cus mak 2; edie eet ae Dampfungsma8 nach Lehr 
teh Cmax — Cmin 1 “ . 
Semllaa ai aren ee sles oe Schwankungsziffer 
Cmitter ergibt sich mit Gleichung (1) zu 
Cmax — Cmi 
Oraitter = Omin +5 Oo 
Die Federzahl ergibt sich somit zu 
C (t) — C mittel + eee ae ey An COS N WW t (11) 
” 1 


Bezeichnet man die Ableitungen nach t mit ’, so folgt aus Gleichung (7) bzw. (9): 


P = wr 
sae tech fila 378 (12) 


co 
” 2 , 2 \’ , Sn r= 
ve’ +2Dyvea'+y P+ oo An COS NT iS 


Gleichung (12) ist eine inhomogene Hillsche Differentialgleichung. Die mathematische 
Theorie dieser Gleichung wurde im Zusammenhang mit dem Problem der Mondbahn- 


4H. Glaubitz und K. Gésele: Entstehung und Frequenzzusammensetzung der Gerdausche 
von Kraftwagengetrieben. Dtsch. Kraftfahrtforsch. Nr. 64, Berlin: 1942. 
5 §. FuBnote 2. 
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berechnung und der Schwingungen von elliptischen Membranen weitgehend ausgebaut, 
jedoch ist es nicht gelungen, exakte Lésungen zu finden. 


Ill. Stabilitét der Lésung der homogenen Bewegungsgleichung 


Welche minimale Diampfung mu8 vorhanden sein, damit instabile Eigenschwin- 
gungen vermieden werden? Im Falle der abschnittweisen Konstanz der Federzahl 
laBt sich die Bestimmung des sogenannten charakteristischen Exponenten oder eines 
gleichwertigen Faktors H relativ einfach durchfiihren. Wir gehen in derselben Weise 
vor wie Meissner*. Die Lésung ist hier abschnittsweise eine gedampfte harmonische 


Schwingung. Mit 


a” +262'+ [A+ yop(t)la=0 (13) 

wobei g(t) = +1, OS tS 2(e—1)a, und p(t) =—1, 2(e—l)asxts2a 
Es gilt my = a9" = (Cysin Nit + C2008 Ny 7) e“** (14) 
Yq = & ste = (C3sin Net + Cy cos Net) e~ ®* (15) 


N,, N;, die Schwingungszahlen in den einzelnen Zeitabschnitten, sind gegeben durch: 
M=VAty—2& Ne= Va—y — & (16) 
Die Lésungen x, und x, miissen im Zeitpunkt t, = 2 (e« — 1) a stetig meinander iiber- 
gehen und auch die Tangente mu8 dort stetig sein, d. h. 
eet ene end. ee (t= ee (eh (17) 
Weiters muB nach dem oben Gesagten nach einer Periode 22 eine Amplituden- 
vergroBerung um den Faktor H eintreten. Es gilt also 
ts (20) — Ae, (0) und 2. (2-2c) = A a (0) (18) 
H ist hiebei ein Faktor, der das Anwachsen oder Abklingen der Schwingung beschreibt*. 


Die vier Bedingungen (17) und (18) ergeben vier homogene algebraische Gleichungen 
in den Unbekannten C, ... C4. Dieses System hat nur dann eine nichttriviale Losung, 
wenn die Koeffizientendeterminante verschwindet. Dies ergibt eine quadratische 
Gleichung fiir den Faktor H: 


H?—22AH+¢2=0 (19) 
wobei 
NY+N2? . : 
A =cos[2(e—1)a Nj] cos[2(2—«)2 2 — yay sin [2 (e — 1) Ni] sin [2 (2—«) 2 No] 
und — e-52(e—1)x 
Es wird HH, , SAA ae 
2 gg 


Wenn Stabilitit vorliegt, also | H,, 2 | < 1 sein soll, mu8 gelten 
| H, + H,| < 2. 


¢ E. Meissner: Schiittelschwingungen in Systemen mit periodisch veradnderlicher Elastizitit. 
Schweizer Bauzeitung 72 (1918). 

* Bei der allgemeinen Behandlung der homogenen Hillschen Differentialgleichung wird besser 
das Floquetsche Theorem verwendet, wobei der charakteristische Exponent fiir die Stabilitat 
der Lésung mafSgebend ist. 
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Aus (19) folgt H, + H, = 24 x2, 


' daher Drs lone? 422 
. oder Aa eaten” (20) 
| Die GréBe| A |, aufgetragen als Funktion von 2, y, 6, ist daher ein MaB fir Stabilitat 


; oder Instabilitaét der homogenen Gleichung. Schneidet 1/x? die Kurve A, so gibt es 
' instabile Bereiche, anderenfalls nicht. 

Die Kurve A wird bei Abwesenheit von Dampfung auch als Eigenwertspektrum 
bezeichnet. Ihr Schnitt mit der Geraden + 1 ergibt die Eigenwerte, die dann in der 
Meissnerschen Karte als Grenzkurven erscheinen. 

In unserem Falle sind die Parameter A, y, 6 alle vom Frequenzverhaltnis abhangig, 
daher werden in dem in Abschn. IV folgenden Beispiel die Kurven A und 1/z? als 


|. Funktionen von y bestimmt. 


IV. Bestimmung der Resonanzbeanspruchung 


Falls der Zustandspunkt des Schwingungssystems in einem stabilen Intervall 
von » liegt, sind nicht anwachsende Eigenschwingungen mdglich. Da in wirklichen 
Systemen immer Daimpfung in irgendeiner Form vorhanden ist, klingen diese Schwin- 
gungen ab, auBerdem werden die Labilitatsintervalle eingeengt. Bei groBerer Dampfung 
kann Labilitét ttberhaupt vermieden werden. 

Daher handelt es sich nun im wesentlichen darum, die Partikularlosung der Glei- 
chung (12), umgeformt mittels der Transformation x = e-Y?*- z (um das Dampfungs- 
glied links wegzuschaffen) 


v,Dt 


” ~~ — D? Pe = ” ‘ 
z +ni(1 +6 S'aqc08n1}2— 0! ge (er +f «P| (21) 
1 m 


zu finden. In Gleichung (21) ist v2 (1 — D?) = »,? und & = &, (1 — D?) gesetzt worden. 
Es wird angenommen, da’ die Eigenschwingung fiir t > oo gegen 0 geht. Aus der 
Volterraschen Integralgleichung, die der Gleichung (21) zugeordnet werden kann, 
wollen wir mit Hilfe der Neumannschen Reihe eine Naherungslosung nach Potenzen 
von &, konstruieren. Da die Funktion ® (rt) im ganzen Bereich endlich bleibt, wei’ 
man, da die Neumannsche Reihe sicher konvergiert. Da &, <1, ist die Konvergenz 
sogar sehr gut, so da man die Potenzentwicklung nach wenigen Gliedern abbrechen 
_ kann. Die obiger Differentialgleichung (21) zugeordnete Volterrasche Integralgleichung 
lautet: 

z(t) =cosmt + 5 sinnt — £,[@ (uw) sin vy, (t — u)z(u)du 
<u, Dt j ( 22) 

pie ag | VI-D* A (u)siny, (tr — u) du 


De 
0 


Die Riicktransformation der Lésung von z nach x geschieht wieder mittels 
“Dt 


ae te i ee er 


In unserem speziellen Fall ist die Erregung und die Federzahlschwankung als Fourier- 
Entwicklung gegeben: 


= oo 
E(t) = >) asin(nt — yn), O(t)= Sa, cosnt 
0 i 
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Die Neumannsche Reihe, die als Naherung eingefiihrt werden soll, lautet 
z(t) = 29 (t) +E. % Ct) 4 Fy? 22 (7) + E523 (tT)... (23) 
Setzt man (23) in (22) ein, so ergibt sich 
Zo 4 by 2 + 12 22... 
= cosv,t + ~ sin W,t — §, ry [(2 cos nu) (Zo (w) + & 2 (uw) + &,? 22 (w) -. Py 
«\4 


: rt 


+ wuDu BS 
-siny, (t —u)du+ asle-? be sin (nur) sin vy, (tr — u)du 
1 0 

Koeffizientenvergleich bei gleichen Potenzen von &, liefert folgenden Satz von 
Gleichungen 


_c vy Du 

: 1 sl : cs 

zo = connst + tain e+ fe [elt 3S ce aes 
0 


a= [ (225 cos) 2 (n) sin (x —u)du 

7 \ 1 

: (24) 
2= — 9; [(2 emcosmu) (wu) sin vy, (t¥ — u)du 

ry 1 


Tv 


err ea= —n{ (3a cosnn)s (uw) sin vy, (t —u)du 
1 


0 


Diese Gleichungen lassen sich iterativ lésen. Fiir z, ergibt sich natiirlich die Losung 
der linearen Schwingungsgleichung mit konstanter Federzahl Cyittei: 


v,Dt 
Alert = Gn Sin (NT — Yn — Yn) — 
Rehr nek oars = : Spm" D 
0 1 1 ¢ 
| ( es ma-™) i pe | (25) 
ohare 
sae v2 — n? (1 — D2) 


Die ersten zwei Terme stellen die Eigenschwingung dar. Infolge der Dampfung klingt 
diese in x rasch ab-und es bleibt nur der dritte Term iibrig. Geht man nun mit diesem 
in die Gleichung fiir z, ein, so ergibt sich die erste Naherung* 


Tv 
antes 7 Qn Sin (NU — Yn — Yn) 
a4,=— — Vv Am COS M U —— ¢ 
1 S(3 m i / v2 5 2 4y,2 D2 n2 (26) 
[apes Be 1 =e 


Wir wollen vorerst nur den konstanten Term a, der Fourier-Entwicklung E.,) beachten 
und a, = 0, fiir » + 0 setzen. Mit dem Ausdruck 
7 sin 
: (ntm+yv)tr—A 
erwin nem ses) = 21 | 
V+ (mimo 
ntmty 
0) 


A = arctan 


* Die untere Grenze des Integrals leistet ebenfalls keinen bleibenden Beitrag in a. 
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vy Dt 
ergibt sich Zool an Dd: on (1 — D2) 
“Dt = 
eee V1 SD ad sin nt — Ain sin nt — Aan 
Siam : Vy oa Db?) ae / D2 = D2 
a a |» fee 
(27) 
wobei 
n—yv 4 n+y 
Ai, = arctan SEY I= pDy Ao, = arctan aD (1 — D2, 


Fiir n = », tritt Resonanz ein, wobei 4,, = 0 wird. 
Durch den zweiten Term verschiebt sich allerdings die Resonanzstelle ein wenig. 
Betrachten wir die Funktion 


at 


[2 1 == sie (n as V4)? | 


Z Vv 
foi? a-v,hr04) 
70 7 
9 \ 
8 
Ya 
2*07 
6 
g 05 
4 OY 
é G3 
2 G2 Y, =z 
5 
& 
7 Qi 3 
Zz 
: 7 
7 7 2 3 |n-yl g 7 2 3 |n-y| 
a 1) 


Die spezifische Dampfung in a) und b) wurde zu D = 0,1 gewahlt. 


Abb. 6. a) Funktion [»,? D? + (n — »,)? (1 — Dy) 2, vy, = l, 2, 3, 4, 5 sind die ersten finf Resonanz- 
stellen und m ist der Summationsindex. 


b) Resonanzfunktion V(,—y, = % D [»,? D? + (n — »,)? (1 — Dy}? 


Waene in Abb. 6a dargestellt ist. Fiir n = + », hat sie ein Maximum von der Grofe 


—p -/1— D?. Fir n = », 4 1 ist ihr Wert nahezu 1 und fiir groBere Entfernung 


yon der Resonanzstelle, d. h. fiir | + »,| >1 wird der Funktionswert kleiner als 1 


Ingenieur-Archiy XTII/2 7 
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Man kann daher den zweiten Term in (27) vernachlassigen, da negative »,-Werte 
unmoglich sind, der Term also nie zur Resonanz kommen kann: 


“Dt n=oco ’ A 
it WG yt eee SALONS Se oD 2 
2 = — € = > ym Pi (1) (28) 
1 oe Tp + (n — 9)? 


Um die Konvergenz der Reihenentwicklung besser ersichtlich zu machen, soll als 
Resonanzfunktion V,, siehe Abb. 6b, 
v1 D vy D 


Va—w) Vox? D? + (n — 1)? (1 — D2) it Vv D? + (n — » V1 — D2)? 2a 
eingefiihrt werden. Die Gleichung (28) schreibt sich dann 
uD 
= eH we Vin — vy) Gn Sin (1 — Ain) (30) 


Die zweite Naherung ergibt sich nun mit Hilfe der zweiten Gleichung (24) und der 
Gleichung (30) zu 


MDT n=co m= if 
a ° ie 2 n— 
23 = — V1 el1—D? 1 = 5a : Y ae Qn Vighee i sin = Mt Ag Ay | 
ja Pe T= De == (n aia m + 1)? 


sin n+ m) 1 = An — As sin (n — m)t + Ain — As 


D | DB 
Ye i- pe + (n+ m — 01)? oe ip? + (n — m-=91)* 


sin (n — m)t — Ain — Aa We i) 
D? 
oe ID? + (n — m — 1)? 
7 n+m-+y Wi = 
A, = arctan ay =(1 D2 , As = arctan ie (1 — D2)? 
es = r= aL — panes et 
A, = arctan ~ ah "1 (1 — D®)? , A, = arctan — a “(1 Dy? | 
Die Resonanzfunktionen V haben die Form 
D 
Vien +m+v.) = = (32) 


Vv12 D? + (n + m+ 11) (1 — D?) 
Wieder kann der Ausdruck mit » + m+ », weggelassen werden, da sein Beitrag 
gering ist. Die GréBe der Resonanzspitze wird dadurch kaum beeinfluBt. Es ergibt sich 
somit naherungsweise 


v, Dt n= co m=oo 
1— D? =) 
De aaa tm 2 On Om Vin —v,) {Vintm— vy) Sin [(n + ie 2 Ais} 
(33) 


+ Van-—m—y,) sin [(n ae Ain — Ag] — Vign—m4 vy Sin [(m — m) t + Ain — As]} 


Die stationare zweite Naherungslésung in x lautet daher 


a 25 QIN=ID) x 


or n=co m=oco 

m= So — (3 5) Be 5m Vv — dag + (8) ceo 2 ton Vrs 

(Vin+m—v, Sin [(m + m) t — Ay, — Ao] + Vaq—m—vy Sin [(0 Ae Ain Aa] 
— Vit n—m+m) sin [(m — m) t + Ain — As]} 


(34) 
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Will man noch eine dritte Iteration anschlieBen, so ergibt sich diese zu 


f, 


a "x ; is 
m=3—(plen 2 On Vin —y, Sin (n t — Ayn) + +(% by 2, 2 Oa Vie a 3 
J 
) 


{Vint m—v,) Sin [(n + m) t — Ain — As] + Ven m—vy Sin [n — m)t — An — 


EV meas sin — in) 2 =A Ce. ei oanl wit; 
1 


1 
{Vin mw) {Vineme1—vy 8in [(n + m +l) t — Ay — A, —Ai] 
fe) eee einai lr Ay — As — Ag] 
— Ving m—14+ vy Sin [(m + m —1)t + Ain + Ae — As]} 
+ Vin—m—vy {(Vin—m-41—vy Sin [(m — 1m + 1) t — Ain — Aa — A] 
So See pe sila — 1) Ai it, = Ao] 
— Vig ae oon Bin G—-m— 1) x + Ai, +Ay —Ao]} 
— Vin—me my {Von—m41—vy Sin [(m — m + 1) t — An — As — Ay] 
+ Vin—m—1—vy Sin [(m — m — 1) t — Ap — As — Ae] 
WW ens ieyg sin (2 — mm — 1) 7 A, Ag Ag]}* 

mit den Phasenwinkeln 


= arctan 


A,= arctan Pee 6 = py? A, We Ur l— (1 — DP 


v1 D v, D 
= nt+m—l— 3 = —m—1— 
sind —pyt A= tee ate 


arctan 
VI D 


n—m—I1lt+n l 


A, = arctan OTTER Pe 


(35 
Fir die Aufstellung einer Naherungslésung ist es wichtig, ein Kriterium zu ge- 


winnen, welches besagt, bis zu welcher Ordnung (3) entwickelt werden muB. 


Wir betrachten den Ausdruck 
K= (3) On Om A+ - Oh Vin—v) Vine maw) Vinimeltevy + Vingmel..thim) (36) 
imal imal 
Ist dieser von der GroBenordnung 1, dann hat der ganze Term wegen des 
Faktors = ; (3) ~" nur mehr einen kleinen Wert gegentiber a, und kann daher, je nach 
der geforderten Genauigkeit, vernachlassigt werden. 


Nun ist es nicht schwierig, die Naherungslésung fiir die periodische Erregung 
E (t) 5 ce sin (nT — Yn) 
1 


anzuschreiben. 
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Es gilt fiir x mit Hilfe von (25) 


(Paes An Sin (NT — Yn — Yn) 


a can » ee 


1— D? 


Die erste Naherung ergibt sich dann sogleich zu 


ave \ An Am 
i (3) i= a v2 3\?, 4%? D2 n? 
(ae) + Foe 


{Vea+tm—w) sin [(n m) T 


Vn 


Wn 


(37) 


As| 


+Vin—m—v,) Sin[(m — mM) tT — Yn — Wn Aa] Vin—m+y) 8in [(m — m) t + yn + Yn — As]} 


Die zweite Naherung wird 


Xe = %1 + 


+ (Bf we 3 2 ee 


V2 une ao 49,2 D2 n2 
1— D2 o = be 


Fee | 4 Viewny. vain ee a ee 
— Ving m—t4yy sin [(m + m —1) 0 + Yn + Yn + 42 — As]} 


st Ga m= va em ee ay SID [1 ete b) Ty, Az SAA 


) 
( 
+ Va—m—1—y) 8in[(n —m —1)t — Yn — Yn ey ee 
— Vin—m—t4my) sin [(m — 1m — 1) + yn + yn + Ay — Ae]} 
( 
) 
) 


— Va-men) (Vaimti-vy sin[(an — m + It ay, — @, — As — AG] 


+ Vie—m—i1—) sin [(n — m — 1) t — yn — Yn As As] 
— Va—m-i+y) Sin[(n — m — 1 t+. + yn + As — Ay]}} 


(38) 


(Vntm—m) { Vintm+i—v,) Sin [(n-+-m+1)t | 


(39) 


Auch hier ist die Konvergenz mit obigem Kriterium (36) priifbar und man kann damit 


unwesentliche Teilschwingungen aussieben. 


V. Beispiel zum Zahnradgetriebe ohne Zahnspiel und statisch bestimmten 


Zahneingriff 
Es werden dieselben Annahmen getroffen wie in 1, 8S. 252: 
Wibersetzarniee sy 0 tte. eee ane eee (eee ieee | 
VAN ahaley AN Ups ortsaub cusicharae in eS Gibran ts 5 Jn ey 0 
Modul es. sct? Suet ard oe 8 i. 6 OE cee ania m = 10mm 


Federzahl (Zahnbreite 6 = 1 cm) 
Trage Masse bezogen auf WaAlzkreis 
Spezitische Damipfings 77, jvasen. eee ee D = Of 


gies 


ha tdi tee ce a vt er d= d,.— 30cm 
WaAlzsprung (am Teilkreis gemessen) ...... fi = 0,005 em 

sie Ses amet C = 116000 kp/cm 
Lees m, = 0,00285 kp sec?/em 
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Der periodische WAlzfehler sei ebenfalls die von Zeman angegebene Funktion f(t) = 


= 0,41 cos wt + 0,235 cos 2, ¢ + 0,088 cos 3 w, t + 0,0195 cos 4a, t + 0,00195 
cos 5@,¢... 


Mit dem Uberdeckungsgrad «¢ = 1,65, nach Dubbel, Taschenbuch fiir den Maschinen- 
bau, ergibt sich aus Gleichung (10): 


Cx) = 0,30 + 1,135 cos t — 0,515 cos 27 — 0,0674 cos 37 + 0,300 cos 47 — 
— 0,180 cos5t... 


Nach tise Coie 1.4556 
Die mittlere Federzahl wird dann nach Gleichung (11): 


Cmittel = 116 000 4 


(1,45 — a 116000 | 0,30 = 150 450 kp/em. 


Daher wird die Federzahl 


On) = 150 450 1 + 0,19 (cos t — 0,455 cos 21+ — 0,0647 cos 37 + 0,265 cos 47 — 
— 0,151 cos 57...) 


Darin ist die Schwankungsziffer € = 0,19. Die in der Klammer stehende Fourierreihe 
stellt die Funktion ®,, = >/a,cosnz dar. Die Koeffizienten a, sind sinngemaB 
: 


abzulesen. 


Die mittlere Eigenfrequenz ergibt sich zu: 


1 /150450 24 
OO, = eee — (210 S0Ge-. 


Das entspricht einer Drehzahl des Rades 1 von n, = 2308 U/min. 


Stabilitatsbetrachtung 
In unserem Fall ist « = 1,65 und mit Gleichung (16) wird 


N+ N,2 


N, = 7650 oe ) N2 = 6335 sec™’, 2N,N, 


= 102! 


A ergibt sich daraus zu 


45.040 17.200 
A= 101 cos 2 =, 01.cos ce 
Cc 


Oc 


5 


10 : 
Setzt man noch ~- =u, so wird 
Cc 


A = 1,01 cos 0,4504 wu — 0,01 cos 0,1720 u. 


Diese Funktion ist in Abb. 7 dargestellt. 
k 2(e-1)x 
Die Funktion — — ee °  schneidet hier nirgends die Kurve|A|, daher sind 


die Lésungen fiir alle w, stabil. Es ist hier auch jene Dampfungskurve eingezeichnet, 


fiir die im niedersten Mitnahmeintervall gerade noch Instabilitat eintritt. Ist der 
Getriebes entsprechend schlecht, so ist das Getriebe an der 


Dampfungsgrad des 
ie eee Stelle w, = 2 @ und deren Umgebung 


2A kaye nicht verwendbar. 
WO)” %/mtb-a77 Wie man aus Abb. 7 sieht, ist die 
6 Gefahr instabiler Eigenschwingungen 
103 mit O- 0007 nicht groB, denn D = 0,001 wird wohl 


fr 


immer vorhanden sein, selbst wenn 

stabi be/D-o keine Olfilmdampfung zugegen ist. 
( Ware die Dampfung so gering, dah 
% instabile Eigenschwingungen mdglich 
sind, so k6nnen sich diese bei w, = const. 
unter der Bedingung M, >71M, unbe- 
bu grenzt aufschaukeln. Ist kein auBeres 
Moment vorhanden, dann kann diese 
Aufschaukelung nur soweit gehen, bis 
ein stabiles Gebiet der Parameterebene erreicht wird, wobei Drehzahlabfall eintritt. 


IAI 


7 
yapaes 
Pzan 


Abb. 7. Stabilitatsdiagramm (EKigenwertspektrum) 


Resonanzbeanspruchung 


Es handelt sich also im folgenden darum, die Resonanzen mit der Storung und 
damit auch die dynamische Beanspruchung zu finden. Bei einem Walzsprung von 
fi = 0,005 cm lautet das Stérglied der Bewegungsgleichung (12) hier: 


+ y2 
He ee + (0,00205 cos r + 0,00470 cos 2t + 0,00396 cos 37 + 
+ 0,00156 cos 417 + 0,000244 cos 5r...). 
Mittels Gleichung (35) kann nun die Partikularlésung von (12) fiir den Stérungsterm 


err ; ‘ 4 
@o = 750459 gefunden werden. Resonanz tritt auf, wenn n = y, ist. Um zu entscheiden, 


welche Glieder der Mehrfachsummen von Bedeutung sind, geht man mit 


&1 

= = 1,804 

und der Resonanzfunktion V nach Bild 6b in das Kriterium (36) und bestimmt die 
GroSenordnung der Zusatzschwingung. Wir nehmen etwa i = 1, », =n =1, 2,3... 
usw., V (0) ergibt sich gleich eins, es wird daher 


Die erzwungene Schwingung sin rt hat also in der Resonanz die 0,902fache Amplitude 
von a, der Verformung durch die Umfangskraft! 


5 : ] 1é& 
Mit n = 2 wird 3 Koja 5 p 22 = 1,804 - (—0,455) = —0,41. 


ea) 
oy 
Die erzwungene Schwingung cos 2 t hat daher bei », = 2 die 0,41fache Amplitude 
VON Go. 
Durch die konstante Umfangskraft werden demnach beachtliche Schwingungen 
angeregt, welche sich den Schwingungen durch ungenaue Flankenform der Zahne iiber- 
lagern. 
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; Betrachten wir die K-Werte fiir i = 2, so ergibt sich mit Gleichung (36) und 
Bild 6b an der Resonanzstelle », = 1 folgende GréBenverteilung in % von ap: 


g Kann’ 100 =F (3) an Vin—1) Vinm—s) 100 
jak ae os i. 
m=1]| 4,07 | —o,094 | 0,065 S 
m=2!| —0,094 0,0298 = — 
m=3| —0,09 = a = 
m= 4 == = i = 
FEY ayae Volt Vin -1)-100 
Aad | 4,07 | << 0n Siew 0-0TSTA lo. = 
m = 2 | —0,094 0,085 = — 
m= 3 | —0,090 zs is 3 
m=4 = = cs am 
(By An Om Vin—1) Vin—m+1)* 100 
m1 | 4,07 | =1),08. « |c-0:00043. |. « ua 
m= 2 | 718.4 | 0,0841 0,000275 — 
31). 0,265) —0,12 = te 
m= 4 0,52 as =: a 
Diese K-Werte sind mt a zu multiplizieren und liefern dann die Amplituden der 


Teilschwingungen. Beschrankt man sich bei der zweiten Naherung auf Glieder gréBer 


als 1%, so sind nur die umrandeten Zahlenwerte zu beachten. In der dritten Naherung 


; : US NG 
sind nur die Terme von der Form (ap) a3 Va—1) Va—1-1) Va—141~-1 von Interesse. 


Die zugehérigen Glieder haben die Form (s. Gleichung 39): 


32 en ee 0,1 
+1 6 +1 —_— 


bee Ng, ‘ n+m—yv 
(3) 4,3 ees Vo tust= yy Sin (x + m—l)r—O arctan” "10 
—— ——S en ee 
0 


: n—m—y 
Bee ie Vi ines avjaellt [(~ —m + 1)7r —0— arctan” 5"! — 0| 
——$S ~~ —_e_—_so ’ 
cer 0 i 
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py he AN I [(~ ayo arctan” —"—"!_9| 
‘ LE aa MBO SETAE! — oct ese 
=I 0 —1 —84,3° 
ee Vin —m-+%;) Vin —m-+-1— 4), sin [i Sell =i 1) tra Uae arctan” "1 0] 
+1 0 +1 84,30 
SOV Gh eats vy) ieee ey) ILS [(@ —m—l)r+0+ arctan” 7" J) 
— ——— es ee ene pacts 
+1 0 —1 843° 
Ay 
=(+) oo ay? Vo) Vay Vo) COS T. 
Die Lésung 2, bei der Resonanzstelle », = 1 wird nun gebildet durch 
ee 
; ; em 
my = — t- 1804 E - 1 sin (tr — 0) + 0,1 - (—0,455) sin (21 — arctan 5) 
+ 0,048 (—0,0647) sin(3 1 — arctan ia} + 0,032 (0,265) sin (4: — arctan 5): ; ; 
=— a [0,902 sint + 0,041 cos 27 + 0,0027 cos 37 — 0,0076 cos 4T...] 
a + |0,0407 sin (« — 0— arctan =) 


+. 0,0407 sin ( —r —0—arctan =) 


— 0,0407 sin (« — 0— arctan i) — 0,184 sin (0 -7— 0—arctan . Se ] ) 


= — © (0407 cos t 
v 


Ogu, ; ee hea 
39 2 O:: cos — a 0,0185 cost. 


Xz — X_ = —1,8048 


%, =—% (1 — 0,902 sin r — 0,0592 cos tr —0,041 cos 2 r — 0,0027 cos 31 
+ 0,0076 cos 47...) 


Zn, Soekirumn ver 
Federzallsthwalkung 


50 


x, 


S 7 
nn ee ee 5 
© 2310 1155 710 57 432 —7,[U/min] 


Abb. 8. Schwingungsausschlag in Prozent der konstanten Vorspannung x) = P/Omittel. Die strich- 
lierte Kurve ergibt sich durch Iteration der Gleichung (12) 
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Von dieser Funktion ist das Maximum zu bestimmen, indem man die Teilschwingun- 
gen graphisch addiert. Dies wird wesentlich erleichtert, wenn man ein Gerat zur 
harmonischen Synthese zur Verfiigung hat. Auch eine Gezeitenrechenmaschine laBt 
dafiir verwenden. Man kann dann auch die Phasenwinkel A genauer beriicksich- 
igen. 

Die sich ergebende maximale Amplitude ist in % von P/Cyitte: in Abb. 8 iiber », 
bzw. der Drehzahl des Zahnradgetriebes aufgetragen. Begniigt man sich mit einem 


vin 


Skolaz Bow! EPUB bBA 1 nbn hs 2 
chains / 


— tn 
By UES BB RZ 1-7 ek Bel VIE D638 Yh Y 


Abb. a: VergroBerungsfunktion V* des linearen Schwingers mit konstanten Koeffizienten bei 
einer spezifischen Dampfung D = 0,1. Die Skala fiir V* mu8 entsprechend der zu berechnenden 
Resonanzstelle », gewahlt werden 


Naherungsverfahren, welches keine unmittelbare Abschatzung der Konvergenz zulaBt, 
kann man in bekannter Weise vorgehen, indem man Gleichung (12) iteriert: 


Lp +2Dvaxy +? a = 


>XLy = 


Omittel Cmittel 


al 


n 
xy" + 2Drxy' + v2 x = y?— Ex v? > an COSNT 
1 


Omittel 


ae IP ei £ I Vo An COS (NT — Yn) 
1 Omitted = T° V(v? — n2)? +4 D2? n2 
2Dyn 
tan yn = oliere usw. 


Die Kurve (x, — %,) ist in Abb. 8 strichliert eingetragen. 


Auf analoge Weise untersucht man den dynamischen EinfluB des periodischen 
Walzfehlers. Hier beniitzt man die VergréBerungsfunktion V* zum Aufstellen der 
Teilamplituden, siehe Bild 9. Das Kriterium zur Aussiebung der unwesentlichen Teil- 
schwingungen wird handlich, wenn man » ~ », setzt. Der Fehler betragt 5°/o). 

‘ Die dynamische Zusatzkraft wird in guter Naherung bei Resonanz, v = 1: 


Payn a Ce) =X (¢)) = 


Cmitter ‘| 1 + & "a, cost] - (0,0108 sin t + 0,00048 sin 27 — 0,00156 cos 27...) 
1 
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In Abb. 10 sind die drei GréBen Payn, C, 2, tiber t aufgetragen. Abb. 11 zeigt den 
Verlauf der Resonanzkurve von « iiber »,. Der Unterschied gegeniiber dem Schwinger 
mit konstanter Federzahl ist bei der hier angenommenen Dampfung nicht sehr groB. 


769.000 


776.000 


OX 
- 000156 C08 20 \ \ 


! 
ae ie tarts 


Se 


Abb. 10. Verlauf der Schwingung z, (t), der Federzahl O (rt) und der dynamischen Zusatzkraft Payn 
beta ee 


rhealineater Schwinger 


Uineater Schwinger lit E(C)= a, cost 


0 7 z 3 ¥ me 

_——— ee ee eee ee -7, 
oo 7240 3620 2416 7810 HUT ~—w,[sec’/ 
oo 2310 155 770 577 432 <7, LYmmn] 


Abb. 11. Resonanzkurve des rheolinearen Schwingers, welcher durch den Wiilzfehler erregt ist 
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Er zeigt sich hauptsachlich in der Krregung durch das konstante Leistungsmoment M,. 
In Abb. 12 sind nun die Diagramme Abb. 10 und 11 zusammen gezeichnet. Man sieht, 
dafs die Vorspannung P praktisch nie ausreicht, um ein Abheben der Zahne bei pian. 
zu vermeiden, wenn ein Zahnspiel vorhanden ist. 


Zleim 
4-420 kp 
4-10 kp 
P-0 
GOO? 
AbhEbLEN Wer 
LIANTIANKED ' 
P-0 
AChECEN 
GOO. eget 


VOrsp. f2* 420 kp 


1 
1 
! 
1 
1 
{ 
1 
' 
‘ 
1 
1 
' 
' 
1 
1 


Vorsp. F; = 160 kp 


7 2 3 4 Y, 


Abb. 12. Resonanzkurve des rheolinearen Schwingers, welcher durch den Walzfehler und die 

Vorspannung P erregt ist. Um ein Abheben der Zahnflanken auch bei v, = 1 zu vermeiden, mite 

die Vorspannung P = 400000 kp betragen. P = 1000 kp wiirde geniigen, um bis zun = 6500 U/min 
spielfreies Laufen zu erzielen 


VI. Einflu8 der statisch unbestimmten Kraftiibertragung 


Die bisherigen Ergebnisse zusammenfassend, kann gesagt werden: Das Zahnrad- 
getriebe mit konstanter Drehzahl nach Abb. 2 ist bei Vernachlassigung des Zahnspieles S 
und der statisch unbestimmten Kraftiibertragung ein Schwinger mit periodisch ver- 
anderlicher Federzahl. Fiir D = 0,001 ergeben sich instabile Schwingungen bei einer 
Kreisfrequenz von w, = 2 w,. Dieser Fall ist allerdings recht unwahrscheinlich, da aus 
Messungen Dampfungszahlen um D = 0,1 bekannt sind. Der Einflu8 der periodischen 
Federzahlschwankung wird daher hauptsachlich in einer schwachen Frequenz- und 
Amplitudenmodulation der abklingenden Eigenschwingung des Systems bestehen. 

Die erzwungene Schwingung, welche durch die konstante Vorspannung und den 
WaAlzfehler entsteht, bleibt dann beschrankt, ist aber groBer als die des linearen Schwin- 
gers mit konstanter Federzahl. Wahrend die Vorspannung dort keinerlei Schwingung 
anregt, geht hier die Federzahlschwankung gewissermaBen als zusatzliche periodische 
Erregung ein. Wenn wir nun den statisch unbestimmten Fall mit C = C (9, x) betrach- 
ten, so ist klar, daB diese Funktion innerhalb des Streifens (Cyax — Cmin) liegen muB 
und die Ubergange von Cy,i, auf Cyax oder umgekehrt nun nicht mehr so regelmabig 
sein werden. Der Eingriffsteilungsfehler f, spielt hier eine groBe Rolle, siehe *. Man 
kann aber jetzt trennen zwischen der periodischen Rechteckschwankung und einzelnen 
Abweichungen (Sté8en) von dieser. Es kommt ganz auf die Herstellungsgenauigkeit 
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an, ob die dem StoB jeweils folgende gedimpfte Eigenschwingung oder die periodisch 
erregte erzwungene Schwingung die gréBere Beanspruchung bringt. Wenn C (9, 2) 
gemessen wird, was durch Verspannen zweier Zahnrider und Messung der Federzahl 
in allen Stellungen leicht geschehen kann, so laBt sich die StoBfunktion experimentell 
ermitteln und ihr Einflu8 auf die Schwingung feststellen. Da jedoch die der Federzahl- 
schwankung folgende Eigenschwingung stark gedimpft ist und Stabilitat bis zu 
D = 0,001 herrscht, kann hiedurch keine wesentliche Veranderung der dynamischen 
Beanspruchung durch die Umfangskraft eintreten. 

Schnellaufende Zahnradgetriebe, die merkliche Kingriffsteilungsfehler aufweisen, 
sind von vornherein untragbar und miissen bei der Erzeugung ausgeschieden werden. 
SchlieBlich ist noch zu beachten, daB die Abnutzung unverhaltnismaBig groBes Zahn- 
spiel hervorruft. 


VII. Beriicksichtigung des Zahnspieles S 


Im weiteren Verlauf der Rechnung wird die Bewegung des Getriebes in zwei Phasen 
zerlegt. Die erste Phase ist durch das Anliegen der kimmenden Zahnflanken gegeben. 
Die Bewegungsgleichung ist bekannt. Die zweite Phase stellt das Durchlaufen des Zahn- 
spieles dar und kann an die erste leicht angeschlossen werden. Zu diesem Zweck bedient 
man sich der Darstellung in der Phasenebene x und z. Fiir die Darstellung der ersten 
Phase geniigt vorerst eine grobe Naherungslosung der Gleichung (12). Diese kann mit 


eager e . 
i Prt iga tiberfiihrt werden in 


aa G U 2 \7 Ep IC! ane. P. v7 - . 
y' +0,2vy' +» c + _ Qn COS ney pes Tas ~ eT Ti, Vo. 
Bei vy = 1 ergibt die Naherungsrechnung, wenn wir die Erregung auf die Terme a, und 
a, beschranken: 


ei a a! ea oN, eee - 
Lea lé Cnerrn x5 cos (t — 90°) = Asint (40) 


und die homogene Gleichung (12) liefert 
Yn = ae-P* cos (t — A) (41) 
worin @ und A von den Anfangsbedingungen abhangen. 


Nach Abb. 10 nehmen wir » = 1, A = 0,0160 cm sowie eine Umfangskraft P = 
= 420 kp an. Um die Bewegungsvorginge wenigstens genahert weiter verfolgen zu 
kénnen, wenn ein Zahnspiel S = 0,005 cm vorhanden ist, wollen wir folgende Ver- 
einfachungen treffen: 


1. Gleichung (40) sei giiltig, unabhingig davon, ob die Zahnflanken des treibenden 
Rades links oder rechts anliegen (nachdem sie das Zahnspiel durchlaufen haben). Der 
WaAlzfehler ist dann in beiden Drehrichtungen der gleiche. 


2. Das Zahnspiel ist immer gleich gro8 fiir alle Zahnpaare. Diese Forderungen 
schranken die Giiltigkeit der Rechnung nicht wesentlich ein, wie noch gezeigt wird. 
Wenn das Zahnspiel durchlaufen wird, gilt die Bewegungsgleichung 
P 
Ye (42) 


~~ ‘Mred 


Man kann nun den Bewegungsvorgang in der x, x-Ebene sehr anschaulich darstellen 
(Abb. 13). 
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Die Kreise K, und K, stellen die Partikularlésung Gleichung (40) dar. Die homogene 
Lésung (41) ergibt eine logarithmische Spirale. Die Gesamtlésung, welche von einem 
beliebigen Anfangspunkt 
der x, x-Ebene ausgehen 
kann, schmiegt sich daher 
immer an den Kreis K,,, an, 
d. h. dieser ist die stationare 
Lésung. Gehen wir etwa vom 
Punkt P, aus, so bewegt 
sich der Systempunkt ent- 
lang K,, bis die Ordinate 
x =O erreicht wird. Von 
hier ab gilt nun die Glei- 


{WM sec 


chung (42). Es ist ay AE [ay 
: Je 
hi aay “ a 
ae (43) 
nas +C,t+C, 


Mit den Anfangsbedin- 
gungen zur Zeit ¢t = 0 


iP, 


Yo = —- >— Yo =-—h Abb. 13. Bewegung des Zahnradgetriebes mit Spiel (Beispiel 

Cmittel Abschnitt V) und konstanter Federzahl Cmittel bei der Reso- 

5, SI nanzstelle », = 1 nach Abb. 11 
wird 
Ie 
ca C26 tn ‘Omittel 
Fiihrt man U=yt+- aaa 
: ° ; : 2a i d : : 

und wu = y+A ein, so wird aus Gl. (43): «=w? oP u~ > , also die Gleichung 


einer Parabel. 


Ist das Spiel S durchlaufen, wird w= S 


S = ted (uy? — Zhu) 


Te 
Daraus 
- 2h (Qh? PSs 
ee Soy 
Die Tangente bei t = 0 also fiir uw = u = 0 ergibt sich zu 
di 
page ae 


An der Ordinate x = 0 ist h = 111 cm/sec abgeschnitten. 


Nach Durchlaufen von S ist « nach Gl. (44) 


215,6 
gtieoe— at 6.4 cm/sec. 


Es ist natirlich der zweite Wert zu nehmen. 
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Nun gilt wieder die Bewegungsgleichung (40), der Kreis K, ist um S verschoben. 
Die Losungskurve nahert sich dem Kreis K, und durchlauft dabei den Winkel «. Dann 
setzt wieder ein Parabelstiick an, das Spiel wird durchlaufen und der Zustandspunkt 
kommt auBerhalb K, zu liegen. Die Lésung ist dadurch nicht exakt periodisch, sondern 
wachst an, kann aber maximal nur den Punkt P, erreichen. Verfolgt man die Lésung 


von P, aus weiter, so ergibt sich ein ae = 127,8 cm/sec und ein A pax = 90,0166 cm, 
welches bestimmt nicht iiberschritten wird, da die Parabelscheitel links von S liegen. 

Die oben aufgestellte Forderung 1 mu8 nicht unbedingt eingehalten werden; es 
ergibt sich dann ein anderer Kreis K, + K,. Indem man die Parabel von P2 aus zeich- 
net, findet man eine obere oder eine untere Schranke fiir Geschwindigkeit und Aus- 
lenkung; im zweiten Fall ist aber die Maximalbeanspruchung durch die Stationar- 
losung (40) gegeben. 

Ist die zweite Forderung nicht erfiillt, also S nicht konstant, sondern sind wegen 
der periodischen Wiederkehr des Zahneingriffes S; Zahnspiele vorhanden, so miBte 
man eine 7-blattrige Phasenebene zur Darstellung der Vorginge verwenden. Aber wie 
kompliziert der Verlauf der Schwingung auch sein wird, sie kann nie zur Instabilitat 
fiihren. Man kann sogar sagen, das die Losungskurve iiber die Kreise K,, , nicht wesent- 
lich hinausgelangen wird, wenn die Zahnspiele klein und nur wenig verschieden sind 
und das eingeprigte Moment konstant ist. Die Dampfung fiihrt dann immer zur Ver- 
nichtung der hineingeleiteten Energie. 


VIII. Dynamische Beanspruchung 


Man mu8 hier dreierlei Beanspruchungen unterscheiden: 

1. Reine Druckbeanspruchung durch die Schwingung. 

2. StoBformige Beanspruchung durch Aufeinanderprallen zweier Zahnflanken. Dies 
tritt auf, wenn das Zahnspiel durchlaufen wird, aber auch, wenn der zweite Zahn 


dazukommt. Denn da die Radkorper die Relativgeschwindigkeit x (am Teilkreis ge- 
messen) haben, gilt fiir den dazukommenden Zahn das Verzahnungsgesetz nicht. 

3. SchlieBlich wirkt an den Zahnflanken selbst noch die Reibkraft Pay, «uw; sie ist 
jedoch bei giinstigem Schmierzustand (~ = 0,01) sehr klein. Der entstehende Schub 
ist gering. Die Griibchenbildung entsteht namlich nicht durch kombinierte Oberflachen- 
beanspruchung Druck plus Schub, sondern infolge feiner Haarrisse, welche durch 
Walzpressung zustande kommen, so da Ol eindringen kann, die Risse von innen her 
erweitert und schlieBlich ganze Teile der Flanke heraushebt. 

Diese StoBbeanspruchung beim Ubergang vom Einzel- zum Doppeleingriff soll 


fs ; : 6 Sere ; 3 2 : 
naherungsweise ermittelt werden. Wir setzen die kinetische Energie ~ am = (wobei 


m, die Zahnmasse und P,, die maximale StoBkraft darstellen) des stoBenden Zahnes 
gleich der bei der Pressung auftretenden potentiellen Energie 
3, cages 
4 oem 
und vernachlassigen die elastische Verbindung der Zahnmasse mit dem Radkérper, 
da sie sehr weich ist gegeniiber der Oberflichenfederung nach Hertz. Die Konstanten 
k, und k, sind gegeben durch 


k 2k 
evs asl fe) 


leg 8 (1 a a Vidi: Sesh ae Poissonsche Konstante 
a 
: , Roc Elastizitatsmodul 
ee ra i 1.3 ekg) an RL Zahnbreite 


& oo a LEV PRR RTS 5, Kriimmungsradien der Zahnflanken. 
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Wird — dem Berechnungsbeispiel entsprechend — der Wert a+ ~ 100 cm/sec ein- 
gesetzt, so liefert dies eine maximale StoBkraft P,, — 505 kp, die also nicht viel gréBer 
ist als die Umfangskraft. . 

Die Hauptbeanspruchung ergibt sich daher aus der Schwingungsrechnung bzw. 
aus Bild 12. Die auf einen Zahn wirkende Kraft ist gegeben durch 


ava = Aas (r;) : Crain (46) 


Geht man mit Gl. (46) in die Hertzschen Gleichungen ein, so lat sich die Werk- 
stoffbeanspruchung bestimmen. VersuchsmaBig zu klaren ist noch der Zusammenhang 
zwischen Schmierzustand und Dampfung, womit dann eine den mechanischen Vor- 
gangen gemaBe Getrieberechnung médglich ist. 

In einer weiteren Arbeit soll der Einflu8 von zwei an das Getriebe angekuppelten 
Rotoren untersucht werden. 

(Hingegangen am 28. Februar 1959) 


Der Tragheitspol des Kreuzschiebers 
Von Q. Bottema, Delft 
Mit 2 Textabbildungen 


Zusammenfassung: Es wird ein einfacher Zusammenhang zwischen der festen Trigheitspol- 
kurve und der Bahnkurve des Schwerpunktes der bewegten Scheibe gezeigt. 


Als Beispiel zu allgemeinen Betrachtungen iiber die Ermittlung des Tragheits- 
poles fiir ein ebenes Getriebe hat Weirich? die Lage dieses Punktes beim Kreuzschieber 
bestimmt. Fiir die Tragheitspolkurve wird in der bewegten Ebene ein Kreis, im festen 
System eine Ellipse gefunden. Die letztere ist zwar durch die Neigung der Symmetrie- 
achsen und die Groen der Halbachsen vollig bestimmt, doch ergeben sich diese als 
ziemlich komplizierte Wurzelausdriicke, so daB die einfache Sachlage nicht klar her- 
vortritt. In dieser kurzen Note zeigen wir, wie die feste Tragheitspolkurve mit der Bahn- 
kurve des Schwerpunktes der bewegten Scheibe zusammenhangt. Es wird dabei keine 
von den verschiedenen Methoden zur Konstruktion des Tragheitspoles benutzt, sondern 
lediglich auf die Definition dieses Punktes zuriickgegriffen. 


Die Bewegung des Kreuzschiebers wird erzeugt, wenn ein Kreis K, (Mittelpunkt //, 
Radius /) innerhalb eines Kreises K, (Mittelpunkt O, Radius 2 /) rollt (Abb. 1). Jeder 
Punkt von K, beschreibt eine Gerade durch O. Ist S der Schwerpunkt des bewegten 
Systems und sind A, und A, die Schnittpunkte von SM mit K,, dann wahlen wir OA, 
und OA, als rechtwinkelige Achsen OX und OY im festen System. Die Lage der beweg- 
ten Scheibe wird durch den Winkel OA,A, = » gegeben; es sei A, S = b, A, S =a, 
wobei diese Zahlen in Richtung A, A, bzw. A, A, positiv gedeutet werden. Wir setzen 
voraus, daB S nicht auf K, liegt, somit ab + 0 ist. 


Fiir die Koordinaten von S hat man 


t,—=12 sin ©, ys = b cos (1) 
Die Bahnkurve von S ist eine Ellipse mit den Halbachsen |a| und |b|; die Sym- 
metrieachsen sind OX und OY. 
Fiir die Komponenten der Beschleunigung von S geht aus (1) hervor 
z= —asing:g?+acosp:¢, y= —bcosy: g— bsing: @ (2) 


1H. Weirich: Zur Ermittlung des Tragheitspoles und der Tragheitspolkurve, Osterr. Ing.- 
Arch. 9 (1955), 230—238. 
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Sind K, und K, die Komponenten der auf die Scheibe wirkenden Kraft K und M ihr 
Moment um O, m die Masse der Scheibe, mp? das Tragheitsmoment um den Schwerpunkt, 
dann ist 

K,=mz,, K,=my,, M=mp' (3) 


Die Wirkungslinie von K hat die Gleichung 
K,x«—K,yi+iM=0 (4) 


oder 
(—bcosy:x+asing: y)g?+ (—bsing: —acosy: y+ p*)p = 0 (5) 


Fiir jeden festen Wert von @ und bei veranderlicher g und g wird durch (5) ein Ge- 
radenbiischel dargestellt. Das heiBt: in jeder Lage des Systems geht die Wirkungslinie 
der Kraft durch einen Punkt 7’, der von der Geschwindigkeit und der Beschleunigung 
unabhangig ist. 7 ist der Tragheitspol. Fiir jedes ebene zwangslaufige Getriebe kann 
die Existenz eines Tragheitspoles in derselben Weise nachgewiesen werden. Man erhalt 
namlich fiir die Wirkungslinie der Kraft immer eine Gleichung, welche in g? und @ 
homogen und linear ist. 

T ist der Schnittpunkt der Geraden 


—bcosg:x+asng:y=0 
und 
—bsing:x—acosp:'y+p?=0 


welche die Wirkungslinien sind fiir g = 0 und g = 0. Die Koordinaten von 7’ sind 
2 


Ci - BURD, ogg P* cos yp (6) 


he 
Wir haben also 
p? p? 

C7 Sleepy I BeAr Ak (7) 


womit ein sehr einfacher Zusam- 
menhang zwischen der Lage des 
Tragheitspoles 7’ und der des 
Schwerpunktes S gewonnen ist: 
T liegt auf OS, so daB OT = k- OS, 

== z. Die feste Trdghettspol- 
kurve geht aus der Bahnkurve des 
Schwerpunktes hervor, durch eine 
Homothetie mit dem Zentrum O und 
dem Faktor k; sie ist eine Ellipse 
mit den Symmetrieachsen OX und 
OY und den Halbachsen | ka | 
und |kb|. 

Wir bestimmen in der beweg- 
ten Ebene ein rechtwinkeliges 
Koordinatensystem (é, 7) mit dem 
Ursprung S; die é-Achse fallt mit 
Abb. 1 SA, zusammen. Man hat dann 


x= sing + 7 cosg + asing 
y = —fcosp + nsing + bcos 


‘ 


Der Trigheitspol des Kreuzschiebers 105 
oder 
§ = xsing — y cosp — asin? y+ b cos? p (8) 
7 = «cosy + ysiny — (a + b) sing cosp 


i Aus (7) und (8) findet man fiir die Koordinaten des Tragheitspoles in der bewegten 
ene 


a—b 


Er = (lb — 1) (asin? g — b cos?) = (kt n(S “F cos 29) (9) 


nr = (k — 1) 


a+b. 
"9 Sin 29. 


Der geometrische Ort von 7' ist 
der Kreis mit der Gleichung 


(SF =p)P shin? == BA. > (10) 
a—b 

9 ? 

ee eee 


mit p = (k — 1): 


Der Mittelpunkt P dieses Krei- 
ses liegt also auf A,A,, und zwar 
so, daB SP = (k — 1) — = 
= (k — 1) MS,d.h. MP=k MS. <7 
Der Radius des Kreises ist Ky~7------ Zn 
|k—1|-1. Abb. 2 


Man kann aus diesen Resultaten leicht ableiten, wie die Tragheitspolkurven erschei- 
nen, wenn man die Bewegung auf ein willkiirliches rechtwinkeliges Koordinatensystem 
bezieht. Die Scheibe bewege sich so, daf die Punkte B, und B, bzw. auf OX und OY 
bleiben; B,B, = 21. Der Schwerpunkt S sei durch die Polkoordinaten r und # bestimmt 
(Abb. 2). Die Symmetrieachsen der Bahnkurve von S sind OA, und OA,; offenbar ist 
der Winkel B,OA, gleich > . Weiterhin ist a =1 + 7, b=1—r. Wir haben also: die 
feste Tragheitspolkurve ist eine Ellipse mit den Halbachsen | k (/ + r)| und |k (1 — r)|, wo 
is a die Neigung der einen Achse ist 2 Die bewegliche Tragheitspolkurve ist 
ein Kreis mit dem Mittelpunkt P; P liegt auf der Geraden MS und MP = k MS; der 


nec ey es 
Radius des Kreises ist a bd: 


(Hingegangen am 13. April 1959) 
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Elastisch-plastische Biegung des nichthomogenen orthotropen 
Bogenstreifens 


Von W. Olszak und S. Zahorski, Warschau 
Mit 2 Textabbildungen 


Zusammenfassung: Die vorliegende Arbeit nimmt auf die Veroffentlichungen * * 7» ® ° Bezug, 
von denen die drei ersten die Grundlagen der Theorie der plastisch orthotropen nichthomogenen Kor- 
per behandeln, wahrend die zwei letzteren Arbeiten die Losung von zwei konkreten Problemen 
betreffen, welche als Anwendung der oben erwihnten Theorie anzusehen sind. 

Gleichzeitig wird in der vorliegenden Arbeit der Versuch unternommen, das durch B.W. Sch af fer 
und R.N. House Jr. ” diskutierte Problem zu verallgemeinern, und zwar auf den allgemeinen 
Fall eines orthotropen nichthomogenen Materials, bei Beriicksichtigung gewisser Voraussetzungen 
beziiglich Belastung und Charakter der Orthotropie sowie der Nichthomogenitat. 


I. 


Betrachten wir ein breites Bogensegment, das durch ein stetig und monoton mit der 
Zeit t anwachsendes Momentenpaar M(t) einer reinen Biegungsbeanspruchung ausge- 
setzt ist. 

Um das so definierte Problem hinreichend zu kennzeichnen, nehmen wir an, dab 
der Bogen als im zweidimensionalen Formanderungszustand sich befindend angesehen 
werden kann; bei Einfiihrung eines zylindrischen 
Koordinatensystems (vgl. Abb. 1), dessen Ur- 
sprung mit dem Kriimmungsmittelpunkt des 
Bogens zusammenfallt, ist dieser Zustand gleich- 
bedeutend mit der Bedingung 


éz = 0; (1.1) 


es gelte auch die Annahme, daf sich die Ab- 
messungen des Bogens wahrend des Belastungs- 
prozesses nicht merklich andern. 
Abb. 1. Der kreisférmige Bogenstreifen Wir nehmen weiterhin an, daf das Material 
durch eine zylindrische Orthotropie gekenn- 
zeichnet ist und da alle mit den elastischen Eigenschaften zusammenhangenden Kenn- 
groBen von den mit den plastischen Eigenschaften zusammenhangenden GréBen 
unabhangig sind, daB hingegen die einen sowie die anderen stetige und monotone 
Funktionen des Radiusvektors r darstellen*. 

Die obenerwahnten Annahmen haben zur Folge, da der allgemeine Fall eines zwei- 
dimensionalen Formanderungszustandes sich auf den einfacheren Fall eines ebenen 
Formanderungszustandes zuriickfiihren laBt**. 

Wir nehmen ferner an, da8 das zylindrische Koordinatensystem r, #, z mit den Haupt- 
richtungen der zylindrischen Orthotropie des untersuchten Systems zusammenfallt. 

Im Hinblick auf den axialsymmetrischen Charakter des untersuchten Problems 
fallen die Hauptrichtungen der Orthotropie mit den Hauptspannungsrichtungen 


* Selbstverstindlich sind nachstehende Kombinationen méglich: elastische Homogenitait — 
plastische Nichthomogenitat; elastische Nichthomogenitat — plastische Homogenitat usf. Auf 
eine dieser Moglichkeiten kommen wir bei den weiteren Betrachtungen in Abschnitt IX zuriick. 

** Beziiglich der zweidimensionalen Probleme in der Theorie der nichthomogenen Kdérper 
siehe Arbeit 4; in dieser werden auch die Kriterien formuliert, welche den Ubergang von dem 
(allgemeineren) Fall des zweidimensionalen Zustandes auf den (einfacheren) Fall des ebenen Zu- 


standes, als einen seiner méglichen Sonderfille, gestatten. In der Theorie der homogenen Korper 
tritt dieser Unterschied nicht auf. 
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zusammen; die Hauptspannungen, somit hier die Radialspannung o,, die Umfangs- 
spannung 0; und die Langsspannung o,, sowie die entsprechenden Formanderungen, 
€, und &, sind unter diesen Umstainden Funktionen des Halbmessers r und der Zeit t. 


I. 


Unter Benutzung der bekannten Ausdriicke fiir die Hauptkomponenten des Form- 
anderungstensors 


F) 
mere ia! rag ae ee (2.1) 
worin die Radialverschiebung wu = u (r, t), sowie nach EKinfithrung der GréBen (vgl. !) 
1 1 /ou U 1 /ou 
e=9 (& + &+ &) = x(3 | =| e=3(F ae 2) 


und bei gleichzeitiger Benutzung der Komponenten des Deviationstensors der Form- 
anderungen in der Form 


e=e—e ( =7;1,2), (2.3) 
erhalten wir die nachstehenden Zusammenhange 

C= & —€ =F(e + 39), 

= —e=>(e— 39), C) 

é;= &— E=— €. 


Die allgemeinen Kompatibilitatsbedingungen fiir die Formanderungen, ausgedriickt 
in Zylinderkoordinaten (vgl. z. B. 4), reduzieren sich in unserem Fall auf eine einzige 
Gleichung, welche die nachstehende Form annimmt 


2 0K Oe, ios; 
l 
r or | or r or 0. (2.5) 


Werden die Formanderungen, die in der obigen Gleichung auftreten, durch ¢ und g, 

unter Benutzung der Beziehungen (2.4), ausgedriickt, so erhalt man 
er — gr +e —3q' = 0; (2.6) 

die ersten und zweiten Ableitungen beziiglich r sind hierin mit einem bzw. zwei Strichen 
angedeutet. 

Die Gleichung (2.6) enthalt zwei unbekannte Funktionen « und q, zu deren Be- 
stimmung noch die Differentialgleichung des inneren Gleichgewichtes in der bei 
fehlenden Massenkraften giiltigen Form 


00; Oy -— Of 


FP erie ne oe 2a) 


herangezogen wird. 
Ii. 


Betrachten wir nunmehr die Beziehungen zwischen dem Formanderungs- und dem 
Spannungszustand im elastischen Bereich. Im Falle einer zylindrischen Orthotropie, 
wenn die Hauptrichtungen des Spannungszustandes sich mit den Hauptrichtungen der 
Orthotropiestruktur decken, kénnen diese Beziehungen in der nachstehenden Form 
(vgl..z.B.. *) 


E, = Ay Or — Aro Oo se Ais oz; 
Et = Aye Or oa Boo Ot =P Ao3 Fz, (3.1) 


Ez = Aye O, + Agg Ot + gg Oz, 
8* 
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oder in der invertierten Form 
0, = Ay & + Aye & + Ais &; 


0; = Ajyy & + Age & + Ags &; | (3.2) 


o, = Aj & + Ags & + Ass &, 
geschrieben werden, wobei die Elastizitatskoeffizienten a;; und die Elastizitatsmoduli 
A;; nunmehr als mit der Lage des betrachteten Elements verinderliche GroBen auf- 
zufassen sind, d. h. A;; = Aj; (r), ai; = ai; (r), und dies mit Riicksicht auf die voraus- 
gesetzte Nichthomogenitat des betrachteten Systems. 
Es kann gezeigt werden, da im allgemeinen Fall der Orthotropie die GréBen A;; 
und a;; durch die Beziehungen 


1 i 
Ay = FT (422 33 — A" 93) « Ay = FT (413 G23 — 12 3); 

1 1 
A, = A (412 G23 — Ag Ao), Ay; = Se (412 413 — 1 M3), 
sed Sis pee a 

ant (211 G33 — 743), sage ay (Q11 G22 — 742), (3.3) 
| G1, Gra, As 
A =| Gy, Age, Aas = det | ai; |, 
413, U3, Ags 


verkniipft sind, unter der Voraussetzung, daB ¢,, e, und e, voneinander unabhangig sind. 

Bei Vorhandensein einer beliebigen (auch zylindrischen) orthotropen Struktur ist 
es i. a. nicht méglich, die physikalischen Beziehungen in ein Volumsanderungsgesetz 
und ein Gestaltsinderungsgesetz aufzuspalten, im Gegensatz zu den bekannten Gesetz- 
maBigkeiten, die aus der Theorie der isotropen Koérper her bekannt sind. Die Arbeit 1° 
behandelt, unter anderem, diese Frage recht ausfiihrlich. Im untersuchten Fall liegen 
jedoch die Verhaltnisse anders. 

Versuchen wir die Beziehungen (3.2) formal auf eine fiir die weiteren Uberlegungen 
besser geeignete Form zu bringen, indem die mittlere Normalspannung 


1 
C= (o, + o; + @;) (3.4) 
sowie die Komponenten des Tensors der Spannungsdeviation 


$8; = 0; — oO (Cc (3.5) 
eingefiihrt werden. 
Bedienen wir uns auberdem der nachstehenden Bezeichnungen 


1 
G, (r) = 5 (2 Ay — 3 Aga Ay Als Aaa) 
1 
G, (r) = 9 (2 As, — 3 Ay, + Ay, + Aig — Ags), 
1 
Gs (r) = 9 (2 Aj3 — 2 Ass — Ay, + Aoy), 
; (3.6) 
A, (r) =_ (2 Ay, + 3 Ay, — 2 Az, — Ayg — 2 Ags), 
1 
Hf, (r) = 3 (2 A, + 3 A, — 2 Ay, — 2 Aj, — Ass), 
1 
A, (r) Sale Aj; + 4A,, — 2 A,,— Aj, — 3 Aj,), 


a a Sf 
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so kénnen unter Beriicksichtigung von (1.1) die Beziehungen (3.2) auf die nachstehende 
Form gebracht werden 


G, (r) e, + Hi, (r) «, 
G(r) & + H, (r) «, (3.7) 


oe eae 3 (") e, + Hs (r) e; 


wenn zusatzlich noch die Bezeichnungen 


u 
Fi(r= 9 (Au Age | Aig Aas); 
4 (3.8) 
Virb= oe (Ay, + 2A,. + 2 Ag, + 3 Ay, + Ais), 
eingefiihrt werden, driicken wir die mittlere Spannung aus als 
2 
Qasr (1), er. tos (i),.e- (3.9) 
Bei Beriicksichtigung von (3.5) kann die Gleichgewichtsbedingung (2.7) in der Form 
00 St — Sr OS, 
OT at OAT or (3.10) 


angeschrieben werden. 

Werden in die obige Gleichung die Beziehungen (3.7) und (3.9) eingefiihrt, so erhal- 
ten wir, unter Beriicksichtigung von (2.4) sowie der Kompatibilitatsbedingung (2.6), 
ein Gleichungssystem, das zur Bestimmung von « und @ dienen kann; dieses lautet 


é'r— yr + e' —3q' = 0, 
Fa VE je IP 
pt Gt V+ Ale| + UF +O) 9 = HG — 4) — (3.11) 
£(G, + G,) +—(H, — H)), 


wobei die Striche Ableitungen beziiglich r bedeuten. 

Das Gleichungssystem (3.11) wurde zwecks allgemeiner Fassung des Problems an- 
gefiihrt. Fiir einen beliebigen Fall der Nichthomogenitat kénnen die entsprechenden 
allgemeinen Integrale in geschlossener Form, ausgedriickt durch elementare Funktionen, 
nicht angegeben werden; aus diesem Grunde fiihren wir im Laufe der weiteren Uber- 
legungen gewisse zusitzliche vereinfachende Annahmen ein. 


IV. 


Die Méglichkeit der Erzielung einer einfachen Losung des untersuchten Problems 
erhalten wir durch Einfiihrung der Annahme der Inkompressibilitat des Materials, 
sowohl im elastischen als auch im plastischen Bereich, d. h. fiir 


(es (4.1) 


Diese Annahme erfordert eine eingehendere Besprechung. In der Arbeit * wurde 
gezeigt (vgl. ebenfalls sowie * und °), daB im Falle eines homogenen Materials bei 
ebenen Formanderungszustanden die Annahme der Materialzusammendriickbarkeit 
im elastischen sowie im plastischen Zustande (was einem Wert des Poissonschen 
Koeffizienten von » = 0,5 entspricht) bei Behandlungen von elastisch-plastischen 
Problemen zu Ergebnissen fiihrt, die den unter Voraussetzung elastischer Kompressi- 
bilitat des Materials erhaltenen Ergebnissen sehr nahekommen. Dies betrifft insbeson- 
dere die Werte der Spannungskomponenten o, und o;. 
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GréBere Unterschiede treten bei der Bestimmung der Spannungskomponente a, sowie 
der radialen Verschiebung u auf. Es besteht jedoch die Méglichkeit der Erzielung genauerer 
Ergebnisse, welche — durch Einfiihrung der Annahme elastischer Kompressibilitat — 
aus der entsprechend einfacheren (inkompressiblen) Losung erhalten werden kénnen, 
indem geeignete korrigierende Multiplikatoren 2 » beziehungsweise 2 (1— y) eingefiihrt 
werden, wobei » dem tatsichlichen Wert des Poissonschen Koeffizienten fiir das gege- 
bene Material entspricht. 

Im Falle eines nichthomogenen und orthotropen Materials wiirde die Richtigkeit 
eines ahnlichen Gedankenganges einer genaueren Priifung bediirfen. Die Lage ware 
aber analog, wenn man annehmen diirfte, daB fiir ein reales nichthomogenes Material 
alle Elastizitatsmoduli A;; durch ,,ahnliche‘‘ Funktionen der Lage gegeben sind, d. h. 
also, daB diese Funktionen sich untereinander lediglich durch unveranderliche Multipli- 
katoren unterscheiden (Fall der sogenannten ,,gleichformigen“‘ oder ,,monofunktionellen“ 
Nichthomogenitat®) ; in einem derartigen Fall kénnte jedes entsprechend kleine Element 
als ein in analogen Bedingungen arbeitendes Element angesehen werden, fiir welches 
die in der Arbeit ? enthaltenen Bemerkungen als zutreffend angesehen werden kénnen. 
Wenn dann also in jedem von diesen Elementen die Werte der Spannungskomponenten 
o, und o praktisch von der Annahme der Kompressibilitét oder Inkompressibilitat 
des Materials im wesentlichen unabhangig sind, so kann auf Grund der Annahme der 
Stetigkeit und ,,Gleichformigkeit’’ der die Nichthomogenitét charakterisierenden 
Funktionen angenommen werden, da o, und a; fiir « =0 sich nicht wesentlich von 
den entsprechenden Werten des kompressiblen Materials (« + 0) unterscheiden. 

Es mége gleichfalls erwihnt werden, da8 die unter Annahme der Inkompressibilitat 
des nichthomogenen Materials erhaltenen Ergebnisse beim Ubergang zum Fall des 
homogenen Materials zu Resultaten fiithren, die aus der Literatur bekannt sind (ver- 
gleiche Abschnitt IX) und die in unseren allgemeineren Ergebnissen als Sonderfialle 
mitenthalten sind. 

Die Annahme der Inkompressibilitat des Materials fiir die orthotrope und nicht- 
homogene Struktur zwingt eine zusdtzliche Bedingung fiir die Elastizitatskoeffizienten 
(bzw. Elastizitatsmoduli) auf, und zwar 


Qy1 + Geo + Ag + 2 (Aig + og + g1) = 0, 
l (4.2) 
Ay, + 2 Ayy + G39 + Pee (G43 + Gog) = 0; 


nach Kinfiihrung entsprechender ,,technischer‘‘ Moduli H,, H, £3, 12, . . ., ¥31 (vgl.z. B.3) 
— welche als Lagefunktionen behandelt werden — kommt man erst zu den nachstehen- 
den Beziehungen 
3 
M12 + 93 + Yay Se 


‘ (4.3) 


a si al Gs as a) Sane, (%31 + 39)?. 


Da jedoch fir tatsaichlich in der Praxis verwendete Materialien die Inkompressi- 
bilitétsbedingungen (4.2) bzw. (4.3) im allgemeinen nicht zutreffen, entsteht die Frage, 
in welcher Weise die Koeffizienten a,;; eingeschatzt werden sollen, um méglichst an 
die Wirklichkeit angenaherte Materialeigenschaften zu charakterisieren. 

Sind die Beiwerte a;; dieser Art, daB die linke Seite der Beziehung (4.2) — zu ver- 
schiedenen Werten (in Abhangigkeit von der Lage) fiihrend — sich nicht viel von Null 
unterscheidet, so kann die Inkompressibilitatsbedingung als angenihert erfiillt betrach- 
tet werden; bei den weiteren Uberlegungen kénnen dann die tatsachlichen Beiwerte 
a;; in die Rechnung eingefiihrt werden. 
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Kine andere Behandlungsart kénnte auf die Einfiihrung einer gewissen Zahl der 
Elastizitatskoeffizienten mit deren Ausgangswerten beruhen, wobei die Wahl der rest- 
lichen derart erfolgen kénnte, da8 die Inkompressibilitatsbedingungen fiir das Material, 
(4.2) baw. (4.3), erfiillt seien. Man kann z. B. die GréBen E,, E,, Es, 1. als gegeben an- 
nehmen, wohingegen die GréfSen v3, und ¥, (welche einen geringeren Einflu8 auf den 
Spannungszustand in der Ebene 1, 2 ausiiben) derart im Einklang mit der Bedingung (4.3) 
zu wahlen waren, dafi im Endeffekt die elastischen Eigenschaften des Materials 
gewissermafen auf dessen Inkompressibilitat abgestimmt werden. 

Im besonderen Fall der monofunktionellen Nichthomogenitaét kann angenommen 
werden, daB mit der Lage des betrachteten Punktes nur die Moduli £; eine Anderung 
aufweisen, da ja gerade diese bekanntlich einen entscheidenden Einflu8 auf den Charak- 
ter der untersuchten Erscheinungen ausiiben; hierbei kénnten die GréBen »;; als unver- 
anderliche Werte gelten; die Erfiillung der Bedingung (4.3) auf Grund entsprechender 
Wahl von »,, und 3, ware dann leicht erreichbar. 


V. 


Fir inkompressibles Material nimmt die Kompatibilitatsbedingung (2.6) die folgende 
Form an 
rp’ +39' =0; (5.1) 
deren allgemeines Integral ist 
B 
P= ga (5.2) 


darin bedeuten A und B von den Parametern des Problems abhangige Konstanten. 
_ Die mittlere Normalspannung o wird unmittelbar aus der Gleichgewichtsbedingung 
(3.10) erhalten, wobei die Beziehungen (3.7), (2.4) und (4.1) zu beriicksichtigen sind: 


PG +4) -— 2&9). (5.3) 


Im Hinblick auf (5.2) erhalten wir weiterhin 


C) B A 0d (BG 
= (G+ @) —-2G+@) +3(G2-44G), (5.4) 
oa (OS gy —A(St ay FO Aa +e, (5.5) 


worin C die Integrationskonstante bedeutet. 
Unter Benutzung von (3.5) schreiben wir die Komponenten des Spannungszustandes 
im elastischen Bereich als: 


oo = f(r) + Agi) +0, 


=F Uf) +e(At@)1— Alo) + (G+ GN) +C, (5.6) 
B 


wobei nachstehende Bezeichnungen eingefiihrt wurden 


G,+4 G,+G 
ts dp = f(r), [= tear ygi(ry: (5.7) 
Befindet sich der Bogen im elastischen Zustande, so werden die Konstanten A,B 
und C aus den folgenden Bedingungen bestimmt: 
b 
[o1_, =9, [o,],_, = 0, | rar See (5.8) 
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wobei M(t) das an den Stirnflachen des Bogens angreifende Moment bedeutet (vgl. 
Abb. 1) 
2 


M [f (0) —f@)I, 
M (g (b) —g(@)I, (5.9) 
M 


[9 (a) f (6) —- 9 (6) F(@)], 


AG 
ne 
Aaaeais 
worin 
L = 2a (f(b) — f (a)] — 28 [g (6) — g (a)] — (6? — a*) g (a) [f (6) — f (@)] + 
+ (6? — a*) f (a) [g (6) —g (@)]- 
Die Ausdriicke f (a), f(b), g (a), g (6) sind als Werte der Funktionen f (r) und g (r) auf- 


zufassen, die durch die Gleichungen (5.7) in den Punkten a und 6 bestimmt werden, 
hingegen sind die GréBen a und f gegeben durch 
b 


5 
: : l 
a=flgr)+(G4+G@)dr, p= {[f) +r (GitGs) ar. (6.10) 
Es sei bemerkt, daB die Komponenten o, und o; auf bekannte Art von einer Span- 
nungsfunktion @ (r) abgeleitet werden kénnen, so da 


1 d@ d? @ 


Oo; = ae Pak: Of ae 5 (5.11) 


dann ist die zusitzliche Bedingung, die das Fehlen von Normalspannungen (herriihrend 
von Zug- oder Druckeffekten im Bogenquerschnitt) ausdriickt 


. * BO d@ : 
Jodr= he dr = °F 


a a 


=i; (5.12) 


den zwei ersten Bedingungen (5.8) gleichwertig. 
Die Lage der neutralen Schicht c wird aus der im allgemeinen transzendenten 
Gleichung 
[o:],_,= 9 


erhalten. 

Erreicht das stetig und monoton anwachsende Moment M (t) einen gewissen Wert, 
M, = M (t,), welchen wir als das erste kritische Moment bezeichnen, so weist das Bogen- 
material erstmalig plastische Deformationen auf*. 


VI. 


In diesem Abschnitt befassen wir uns mit dem véllig plastischen Zustand, vorlaufig 
unter Ubergehung der Analyse der Zwischenstadien, welche beim Ubergang des Bogen- 
materials vom elastischen in den plastischen Zustand auftreten. 

Wird das Material durch die sogenannte ,,ideale Plastizitiat‘‘ charakterisiert, so kann 
die verallgemeinerte Huber-Misessche Plastizitatsbedingung fiir den Fall eines ebenen 
Formanderungszustandes und eines inkompressiblen Materials [vgl. Arbeit 7, Glei- 
chung (4.7)], unter Beriicksichtigung der Axialsymmetrie des Problems, wie folgt dar- 
gestellt werden 


(o, — oy)? = 4 [K (r)}*; (6.1) 


* Der Wert des Moments M,, sowie die Art wie auch die Stelle des Auftretens der plastischen 
Zone werden (gemaf} den Annahmen in Abschnitt 1) weiter unten besprochen. 
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es gilt dann 


O, 0; =) 2K () "wenn o; > 'o;; | (6.2) 
Op 0, aa (1). Wenn, = ea i 
Diese Ausdriicke, welche formal mit der Plastizitatsbedingung fiir das isotrope und 
nichthomogene Material (vgl. * °) im Einklang stehen, haben jedoch eine ab- 
weichende physikalische Bedeutung. Der Ausdruck K (r) ist zu verstehen als 


[K (NP = 4 Eos (r) Ku (r) = att Bia (1) Koy (7)? ca 
_worin (vgl. 7) 
Bis = COLON: ONCE * Coa 
2 Ky (1) = ar IAT a TGA ae pe 
2 Ky (1) = aor eee TaaT: | 


Q, (r) Q: (r), Q- (7) bedeuten hierin die Plastizitatsgrenzen, die mit den Hauptrichtungen 
der Orthotropie 7, ¢, z im Zusammenhang stehen. 

Da unterhalb der neutralen Schicht, die durch den Radius c gekennzeichnet ist 
(also in der Nachbarschaft von r = a), o, > o,, hingegen oberhalb der neutralen Schicht 
(also in der Nachbarschaft vonr = b) o, > o; gilt, ergibt die Gleichgewichtsbedingung 
(2.7) gemeinsam mit der Plastizitaétsbedingung (6.2) folzende Lésungen: 

Seo <7 "0: 


6, = h (r) — h(a), (6.5) 
o,=h(r)+2K (r) —h (a); 
fOr O< 770% 
OP — hr) +h (d), | (6.6) 
a= —h(r)—-2K(r) +h), J 
wobei die nachstehende Abkiirzung eingefiihrt wurde: 
h(r) = 2 (= ar, (6.7) 
Die Integrationskonstanten wurden aus den Randbedingungen 
en = 0, bel aN 0 (6.8) 


bestimmt. 

Die Lage der neutralen Schicht wird aus der Stetigkeitsbedingung fiir die Span- 
nungskomponente o, beim Durchgang durch die neutrale Schicht bestimmt. 

Wir erhalten fiir c im allgemeinen Fall die nachstehende transzendente Gleichung 


h (0) = [h (a) +h (0). (6.9) 


Das Moment M, = M (ts), welches wir das dritte kritische Moment oder kurz das 
Grenzmoment nennen, entspricht der volligen Plastifizierung des Bogenmaterials (vgl. 
z. B. 12); dieses wird durch den Ausdruck bestimmt 

4 : 
Mo Me = | [2 (r) 42K (r)—h(b)]rdr — [ (hr) + 2K (r)—h(a@]rdr, (6.10) 


wo an Stelle von c der aus der Gleichung (6.9) erhaltene Wert einzusetzen ist. 
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Das Moment M,= M,, hangt lediglich von einer einzigen Materialfunktion K (r), 
welche die Nichthomogenitit des Materials darstellt, sowie von den durch die Radien 
a und b bestimmten Abmessungen des Bogens ab. Erreicht die auBere Belastung in 
Form des Moments M(t) einen Wert, der dem Wert M,, gleichkommt, so mu8 angenom- 
men werden, da8 das Grenztragvermégen des Bogens erreicht wurde. 


Vil. 


Betrachten wir nunmehr die Zustinde, welche dem Ubergang des Bogens in den 
vollig plastischen Zustand vorausgehen. 

Nehmen wir an, daB die Veranderlichkeit der elastischen Eigenschaften des Materials 
G,(r) und G(r) sowie die Veranderlichkeit des Plastizitatsmoduls K (r) so beschaffen 
sind, daB die Plastifizierung des Bogens 
von ,,unten“ her beginnt, d.h. also von 
jener Zylinderflache, die durch den 
Halbmesser r = a bestimmt wird. In 
diesem Fall erreicht der Unterschied 
der Spannungen og, und o, (a, > o,) — 
mit dem Anwachsen des Moments M (t) 
: auf M, — einen kritischen Wert, der 

I Plastische Zone aus der Plastizitatsbedingung (6.2) her- 

Abb. 2. Die elastische und die plastischen Zonen vorgeht, und bewirkt in der ,,unteren“‘ 
: Partie des Bogens das Auftreten einer 

plastischen Zone, die zwischen den Flachen r = a und r = 9, gelegen ist (vgl. Abb. 2)*. 

Wir erhalten auf diese Weise die Erscheinung des ,,gehemmten“ plastischen Flie- 
Bens, eingeschrankt durch das Vorhandensein einer elastischen Zone, welche an die 
plastifizierte Zone angrenzt; diesen Zustand nennen wir den ersten elasto-plastischen 
Zustand. . 

Den Wert des Moments M, bestimmen wir aus den Gleichungen (5.6) und (6.2) 


K (a) 4a? 


II Plastische Zone 


Elastische Zone 


Ms = Ga) +G,(a) (BF 2 AF al) ° (72 
worin unter Beriicksichtigung der Zusammenhange (5.9) die Bezeichnungen 
Z| ; B 
At =, BY (7.2) 


eingeftihrt wurden. 

Die Spannungen in der plastischen Zone werden durch die Ausdriicke (6.5), in der 
elastischen Zone durch die Ausdriicke (5.6) dargestellt; die in den letztgenannten Glei- 
chungen auftretenden Konstanten A, B und C haben jedoch nunmehr eine andere 
Bedeutung ; deren jetzige Werte A,, B, und C, werden namlich aus der Randbedingung 
und den Stetigkeitsbedingungen fiir die diesbeziiglichen Spannungen sowie fiir die 
Verschiebungen an der Grenze r = 0, berechnet: 


olan = 0, Loe <= ees ) ho = [ofl : (7.3) 


Die oberen Indizes e und p nehmen Bezug: e — auf die Spannungen in der elastischen 
Zone, p — auf die Spannungen in der plastischen Zone. 
Die Ausdriicke fiir A,, B, und C, sind mit Riicksicht auf deren komplizierten Bau 


im allgemeinen Fall der Nichthomogenitaét in Abschnitt X (Anhang) zusammen- 
gestellt. 


* Die Bedingungen, welche im oben genannten Fall von den Funktionen G, (r), G(r) und K (7) 
erfillt werden miissen, besprechen wir in Abschnitt VIII. . 
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Die das elastische Gebiet vom plastischen Gebiet trennende Grenzschichte r = 01 
wird bei gegebenem Moment M aus der Bedingung 


Q1 b 
— M=)ofrdr+ |orrdr. (7.4) 
ermittelt. 

Mit Zunahme der Zeit ¢, sowie des mit dieser verbundenen Moments M , erreicht — 
fiir den Wert ¢ = t, und den Wert des Moments M, = M (t,), der per analogiam als 
zweites kritisches Moment bezeichnet wird — der Unterschied der Spannungen o, und o; 
(o, >o,) den zweiten kritischen Wert, der durch die Plastizitatsbedingung (6.2) bestimmt 
wird; dies bewirkt das Erscheinen einer plastischen Zone im ,,oberen‘‘ Teil des Bogens; 
diese liegt zwischen den Zylinderflachen r = @, und r = 6 (vgl. Abb. 2). 

Die Erscheinung des FlieBens des Materials in den beiden plastischen Gebieten, die 
sich in das Innere des Bogens hinein fortpflanzen, wird durch die Existenz des angren- 
zenden elastischen Gebietes gehemmt, welches dazwischen gelegen ist, 0,< 1< 0; diesen 
Zustand nennen wir den zweiten elasto-plastischen Zustand. 

Der Wert des Moments M, betragt 


K (b) 4b? 


M2 =) +0) AFB’ (7.5) 
wobei die Bezeichnung 
A 
A,* = a : Poa — (7.6) 


eingefiihrt wurde. 

Die Spannungen in der elastischen Zone werden durch die Ausdriicke (5.6), in den 
plastischen Zonen hingegen durch die Ausdriicke (6.5) und (6.6) bestimmt. 

Die neuen Konstanten A,, B, und C, berechnen wir aus den Stetigkeitsbedingungen 
fiir die diesbeziiglichen Spannungskomponenten und die Verschiebungsfelder auf 
Grund der nachfolgenden Bedingungen 

Lora, mb lor | atyi8 Lox Waite ts Lot leap; 
; 3 (7.7) 
CA ee ae ee kee weston ia 


F 2 


Deren Werte sind gesondert in Abschnitt X (Anhang) zusammengestellt. 

Die Bedingungen (7.7) (vier bei nur drei Konstanten A,, B, und C;,) liefern eine 
zusatzliche Beziehung, welche die Grenzen zwischen den plastischen Zonen und der 
elastischen Zone, 0, und 0,, verkniipft, welche gemeinsam mit der Beziehung 

Qi Qe b 
= tf.) ofrdr+ | orrdr+ | of rar (7.8) 
a ou Qe 
die gleichfalls bendtigten GréBen e, und g, fiir einen gegebenen Wert des Moments M 
zu ermitteln gestatten. 


Vu. 


Wir gehen nun dazu iiber, diejenigen Bedingungen fiir die Funktionen, die die ela- 
stischen und plastischen Eigenschaften des Materials charakterisieren, festzulegen, 
damit die Forderung nach gleichzeitiger Plastifizierung des. ganzen Bogenquerschnittes 
erfiillt ist, und zwar in jenem Augenblick, wenn iiberhaupt die ersten plastischen De- 
formationen auftreten, d. h. also beim Anwachsen der auf eren Belastung M (t) auf den 
Wert des ersten kritischen Moments /,. 
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Diese Forderung ist durch technische und wirtschaftliche Vorteile begriindet (vgl. 
89) im Hinblick auf die gleichformige Anstrengung des Materials der ganzen Kon- 


struktion. 
Wir stellen nunmehr die Forderung, da die fiir den elastischen Zustand ermittelten 


Spannungen (5.6) fiir jedes 7, (a<r<b), die Plastizitatsbedingung (6.2) identisch 
befriedigen : 


lor — 0, |= sar 1G ) +@2 (r)] |B — 2 Art|=2 K (r), (8.1) 
wobei A und B durch die Ausdriicke (5.9) festgelegt sind. 
Die Identitat (8.1) fordert, daB 
G, (7) 2G. (r) 4r? 
Shan ae LE ae Ne (8.2) 


Den GréBen A und B kommt die Rolle von Materialkonstanten zu, welche den 
Wert des kritischen Moments @, = M,, beeinflussen, wobei dieses den ganzen Bogen- 
querschnitt vollig in den plastischen Zustand tiberfiihrt; dieses betragt 


A= A*M,, B= Bt M,. (8.3) 


Im Hinblick auf die Singularita&t der durch die Identitat (8.2) bestimmten Funktion 
y (r), welche im Punkt r = 1p, 


ee eee en GB ge) 8.4 

4 2 Veoater’ oe 
auftritt, ist das Erfiillen der Bedingung (8.2) praktisch schwer zu verwirklichen, wenn 
die Annahme der Stetigkeit der Funktionen K (r) sowie G, (r) + G@, (r) aufrecht- 
erhalten wird. 

Ks kann leicht festgestellt werden, daB fiir r <r, die Funktion y (7) eine wachsende 
Funktion, hingegen fiir r > 7) eine abnehmende Funktion ist. Die Feststellung, ob das 
Intervall [a,b] den singulaéren Punkt 7 einschlieBt, ist nur bei gegebener Funktion 
G, + G, moglich und fiihrt daher das Problem — im Falle einer negativen Antwort auf 
die so gestellte Frage — auf eine entsprechende Wahl von K (r) zuriick*. 

Das Bestehen der Singularitat der Funktion y (r) la8t sich dadurch erklaren, daB 
im elastischen Gebiet die Umfangsspannungen o, — beim Ubergang durch die neutrale 
Schicht — ihr Vorzeichen zwar andern, dabei aber stetig bleiben, wohingegen im pla- 
stischen Bereich dieselben Spannungen — beim Durchgang durch die neutrale Schicht 
— durch eine ausgepragte Unstetigkeit gekennzeichnet sind (Vorzeichenanderung 
sowie sprunghafte Anderung des Wertes der Spannungen). 

Legen wir nunmehr die in Abschnitt VII angesagten Kriterien fest, welche die Funk- 
tionen G,, G, und G; erfiillen miissen, damit die Bildung der plastischen Zonen in der 
oben beschriebenen Weise vor sich gehe. 

Man kann auf Grund von (8.2) feststellen, daB die hiefiir notwendige und hinreichende 
Bedingung in der Form ausgedriickt werden kann, daB die Funktion 


EK (r) 73 y (r) K (7) 85 

(4, +4) (B24, — BiG, @) +4, M7) ee 
eine wachsende Funktion im Intervall [a, 6] sein muB. 

Ist die Funktion R (7) eine in diesem Intervall abnehmende Funktion, so tritt das 


plastische Gebiet zuerst an der ,,Oberseite‘‘ des Bogens auf, d. h. auf der durch den 
Radius r = b bestimmten Zylinderflache. 


EN hee 


* Im nachsten Abschnitt befassen wir uns mit diesem Problem ftir ein Material von besonderen 
nichthomogenen Eigenschaften. 
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Die Bestimmung der Spannungsverteilung und der entsprechenden kritischen Mo- 
mente M, und M, kann dann in Anlehnung an ein Verfahren vorgenommen werden, 
welches dem im Abschnitt VII dargestellten analog ist. Es fallt nicht schwer, sich zu 
iiberzeugen, da} das Grenzmoment M,, = M, (welches das Grenztragvermogen des 
Bogens bestimmt) in beiden Fallen dasselbe ist. 

Selbstverstaindlich entspricht die Annahme R (r) = const dem Fall einer gleich - 
zeitigen Plastifizierung des ganzen Querschnitts; somit ist dann M Pee 

Ks sei noch folgendes bemerkt: Ist der Verlauf der Funktionen G, (r) + G, (vr) und 
K (r) festgelegt, so geniigt es — bei der Analyse der Entstehung und Fortpflanzungs- 
richtungen der plastischen Gebiete —, anstatt die Funktion R (r) zu untersuchen, fest- 
zustellen, ob der Wert des nach der Gleichung (7.1) berechneten Moments MV, erdBer 
oder kleiner ist als der Wert des Moments, welches aus der Beziehung 


K (b) 4b? 


M, = “G0 +6,0) BAP BY 


(8.6) 


fiir r = b ermittelt wird. Der kleinere Wert von M, deutet die Stelle des Auftretens der 
ersten plastischen Zone an. 


Diese Uberlegungen fiihren zu den folgenden Schliissen: 
a) die Entstehung des plastischen Gebietes beginnt an der ,,Unterseite’ des 
Bogens, d. h. fiir r = a, wenn 


K (b) (2 A* b2 — B*) a? K (a) 
G, (0)+4, (6) ~ (B¥—2A4*a)0  [G,(@) +G, (a) 


(8.7) 


b) die Entstehung des plastischen Gebietes beginnt an der ,,Oberseite des 
Bogens, d. h. fiir r = 6, wenn 


K (b) (2. A* b? — B*) a? K (a) 

G, (b) + G, (b) S (B* — 2 A* a?) b? [G, (a) + G,(a)] ’ 

c) das plastische Gebiet tritt gleichzeitig an der ,,Unterseite‘‘ und an der ,,Ober- 
seite‘‘ des Bogens auf, wenn die vorhergehende Ungleichung zur Gleichung wird. 


(8.8) 


IX. 


Die in den Abschnitten VII und VIII erhaltenen Ergebnisse besprechen wir an Hand 
eines Beispiels eines Materials besonderer Art, und zwar eines elastisch homogenen 
Materials, fiir welches somit 


G, =const, G,.— const, G, + G, = G, = const, (9.1) 


das jedoch gleichzeitig ein plastisch nichthomogenes Material ist, wobei die Funktion 
K (r), welche die plastischen Higenschaften des Materials charakterisiert (vgl. Ab- 
schnitt VI), vorerst noch nicht festgelegt wird. Ein solches Material ware somit in 
Klasse (II) unserer grundlegenden Klassifizierung ° einzureihen. 

Eine derartige Annahme kann zum Beispiel von Bedeutung sein, wenn ein thermisch 
behandeltes Material oder ein wahrend des Belastungsvorganges erwarmtes Material 
betrachtet wird (unter der Voraussetzung, da die veranderliche Temperaturverteilung 
keine merklichen zusitzlichen Warmespannungen hervorruft). In solchen Fallen werden 
die elastischen Eigenschaften kaum beeinfluBt, wohingegen die plastischen Higen- 
schaften in der Regel merkliche Ortsabhangigkeit aufweisen werden. 

Bei der Berechnung der Komponenten des Spannungszustandes gemaf (5.6) erhalten 
wir die uns aus der Literatur bekannten Ausdriicke (vgl. 1”) 
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4M /a®b? b r a 
o, = — In—+ 6?In=-+ a? In—], 
N ( r2 a b | (9.2) 
/ 2 h2 
a= — = aS In + b2In> + a?In< +0? —a!), 


worin 
2 


b 
N = (0? — a)? — 4a? b*(In-~) 
die Spannungen hangen also nicht von den elastischen Higenschaften des Materials ab, 
aihnlich wie etwa im Fall eines isotropen Materials. 
Bestimmen wir weiterhin die Stelle, an welcher die plastischen Deformationen 
zuerst auftreten. Die Materialkonstanten A* und B* nehmen die Werte an 


4 (b2 — a?) 


A* = b\a7 
Go (b? — a?)? — 4a? b? (in 3) | 
b (9.3) 
UGTA 2102 a 
B= “ 
b\2 
Go (b? — a*)? — 4a*b? (m i) | 
wobei immer b >a ist. 
Beim Vergleich des ersten kritischen Moments M,, gema8 Gleichung (7.1) 
x(a [2 a9 4028 (n 2) | 
{ (a — a’)? 4a Nass 
ME 2 (9.4) 


ae 
— 2|2 b2 ae (b= — a’)| 


und gemaé8 Gleichung (8.6) (vgl. Abschnitt VIII) 


K (b) iG — a?)? — 4a? b? (m =) | 
; [My] _, See ; co te 2 atin] ’ (9.5) 


wenn dabei K (b) = K (a) (die Funktion wachst mit dem Radius r), gelangen wir zu 
dem Schlu8B, da8 die Plastifizierung immer zuerst von r = a ausgeht, d. h. also auf der 
Seite der gréBeren Kriimmung beginnt. 

Wenn hingegen K (b) < K (a) (die Funktion nimmt mit dem Radius r ab), so be- 
ginnt die Plastifizierung bei r = b, d. h. also auf der Seite der kleineren Kriimmung, 
wobei jedoch 

(b? — at) — 2a*In~ 
K (6) < K (a) ; (9.6) 
2b*In > — (b? — a?) 


gelten mu; der Fall des Gleichheitszeichens ist gleichbedeutend mit dem gleichzeitigen 
Auftreten der plastischen Gebiete auf beiden Seiten des Bogens. (Dieses Ergebnis darf 
jedoch nicht mit der gleichzeitigen Plastifizierung des ganzen Bogensegments ver- 
wechselt werden.) 

Wir gehen gegenwartig dazu iiber, den Charakter der Veranderlichkeit der Funk- 
tion K (r) zu bestimmen, welche die Forderung der gleichzeitigen Plastifizierung des 
ganzen Querschnittes erfiillen soll. Die Funktion y (r), welche durch die Beziehung (8.2) 
bestimmt ist, nimmt die Form an 


b\2 
Go lo — a’)? —- 4a? b? (im <) | v2 


Yo (rT) = (9.7) 


2M, 2a? b?In — (b? — a?)r? 


Elastisch-plastische Biegung des nichthomogenen orthotropen Bogenstreifens 119 


Wenn bei gegebenen Werten a und b der Wert 


es et 
cara 2 a? b? In — 
B a 
r= lhe b? — a2 (9.8) 
nicht innerhalb des Intervalls [a, b] liegt, dann ist y (r) eine ausgesprochen zunehmende 
oder abnehmende Funktion in diesem Intervall. 
Aus der Beziehung (8.5) erhalten wir die nachstehende Bedingung, welche zur 
Erfillung der Forderung nach gleichzeitiger Plastifizierung notwendig ist 


K (r) wo (r) = const. (9.9) 


Liegt eine wohl bestimmte Funktion K (r) vor, so sind wir, unter Verwendung der 
in den Abschnitten VI und VII abgeleiteten Beziehungen, imstande, die Kompo- 
nenten des Spannungszustandes o, und og; fiir die elastischen und plastischen Zonen 
zu bestimmen. 


Anhang 
X. 


Fiir den ersten elasto-plastischen Zustand nehmen die Integrationskonstanten A,, 
B,, C,, welche aus den Bedingungen (7.3) ermittelt wurden, die Gestalt an 


Ay = 7-4 2K (01) [f (6) — f (ond) + ga [Gs (0s) + Ge (0) [ (au) — (@)I}, 
B, = Z- 42K (on) [g (0) — 9 (au)] + (Gr (ou) + Ge (ard) (on) — & (a)]}, 
C 


1 = 712K (on) Uf Og (01) — f (os) 9 O)] + gx (Or (0x) + Galo) [h (0n) — 
h (a)) 9 (6) + [Gy (@1) + Ge (01)] ih (0 1) —A (a) ft ( (b)f, 


(10.1) 


LI, = =F [Gy (01) + Ge (01)] [gy (6) — 9 (01) ] — [G1 (@1) + Ge (@1)] fAGy == 740)" 


Fiir den zweiten elasto-plastischen Zustand werden die Integrationskonstanten 
A,, B,, C, gemaB den Bedingungen (7.7) in der nachstehenden Form erhalten 


—— 


Ay = 7 { [h 0) +h (a) —h (92) — h (on) ~ 2K (02)] — 
Sate er aeN fi (Q2) — f (@1) +536 (es)||. 
siege yl) Pha) = B(ee) = Be) 2 (eal ost — 


— 201° Groh [9 (ee) —9 (01) +4 (02)]}; 


| 
Cy = qo {pala [h (01) — h(@)] + [h (6) + (a) — (02) — (Qu) — 


(10.2) 


wn 


— 2K (02)] [9 (01) — @:° f (e1)] + 2 ex” Ff (ex) as [9 (@2) — 
es Fen 100+ =o (03)|t, 
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worin 


Ly ={ or[f (es) — f(x) + ge@ (02)| — Ig (02) — 9 (0s) + 4 (of, 
G(r) = @, (") + Ge (") 


dabei besteht die zusitzliche Beziehung zwischen @, und @, 


4G (0) [1— 25] (@) + (6) — b (02) — B (0)] = 08° If (os) — F (0s)][K (es) + 


02” 
2 G(e 


G (09) 2 Glos) (10.3) 
+ Ge) K(o,)| — [9 (ex) — 9 (e1)]| K (02) + G2 Gen K (o,)|. 
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Die Parallelogrammplatte mit Einzellast 
Von E. Tungl, Wien 
Mit 2 Textabbildungen 


Zusammenfassung: Eine parallelogrammférmige Platte mit beliebigem Seitenverhaltnis und 
Winkel ist allseitig frei drehbar gestiitzt und tragt an beliebiger Stelle eine zur Mittelebene normal 
gerichtete Einzellast. Es wird die im Rahmen der Kirchhoffschen Plattentheorie mathematisch 
strenge Lésung fir die Durchbiegung, d. h. also die Greensche Funktion des Problems, ermittelt. 
Sie ist in Form einfach unendlicher Reihen als Funktion der Ortskoordinaten angegeben, wobei 
gewisse Reihenkoeffizienten aus Gleichungssystemen zu berechnen sind. 


Einleitung 


In der Fachliteratur sind zahlreiche Abhandlungen zu finden, die sich mit 
parallelogrammférmigen Platten befassen. In letzter Zeit werden besonders Modell- 
versuche bevorzugt. Homberg und Marx? legen als Ergebnis von systematisch 
durchgefiihrten Modellmessungen ein Buch mit Diagrammen und Tabellen fiir Haupt- 
momente und MomenteneinfluBwerte vor. Die angegebenen Groen wurden zum Teil 
von Vogt tiberpriift, der selbst auch an Modellversuchen? arbeitet. Fiir die Bemessung 
schiefer Platten werden auch Ergebnisse der Bruchlinientheorie nach Johansen® oder 
Bemessungsregeln und Schatzwerte+* beniitzt, die an Hand von verfiigbaren Zahlen- 
ergebnissen zusammengestellt wurden. 

An elastizitatstheoretischen Untersuchungen sind meist Naherungslésungen im 
Rahmen der Kirchhoffschen Plattentheorie vorhanden. Ehasz® vergleicht drei 
Naherungsmethoden, von welchen sich die Differenzenmethode als die zweckmaBigste 
erweist. Allgemeine Beziehungen zur Differenzenrechnung sowie Zahlenergebnisse fiir 
parallelogrammférmige Platten geben Vogt’, Jensen*®, Favre® und Nielsen?®. 
Es werden zwei- oder vierseitig frei drehbar gelagerte oder eingespannte Platten ver- 
schiedener Seiten- und Neigungsverhaltnisse untersucht und auBer Durchbiegungs- 
und Momentenwerten auch Hauptmomententrajektorien und EinflufBfelder* gebracht. 
Das Galerkinsche Verfahren wird von Guzman und Luisoni" beniitzt. Rongved” 


1H. Homberg und W. R. Marx: Schiefe Stabe und Platten. Diisseldorf: Werner Verlag 
1958. — Hauptmomente und deren Richtungen einer 40° schiefen Platte im Vergleich mit denen 
einer quadratischen Platte. Beton- u. Stahlbetonbau 58, 34 (1958). 

2 H. Vogt: Das statische Verhalten von zweiseitig aufgelagerten schiefwinkeligen Einfeld- 
platten. — Die statische Behandlung von schiefen Plattenbriicken tiber mehrere Felder. Beton- 
u. Stahlbetonbau 50, 291 (1955), 58, 209 (1958). 

3K. W. Johansen: Pladeformler. Kobenhavn: Polyteknisk Foreining, 2. Ausg. 1949. 

4H. Risch: Fahrbahnplatten von StraBenbriicken. Deutscher Ausschuf fir Stahlbeton, 
H. 106, Berlin: Wilh. Ernst u. Sohn, 4. Aufl. 1957. 

5 H. Vogt: Die Auflagerkrafte von zweiseitig aufgelagerten schiefwinkeligen Einfeldplatten. 
Beton- u. Stahlbetonbau 51, 184 (1956). — Beitrag zur Bemessung schiefwinkeliger Plattenbriicken. 
Bauplanung u. Bautechnik 11, 83 (1957). 

6 FL. Ehasz: Structural Skew Plates. Proc. Amer. Soc. Civ. Engrs. 71, 805 u. 1436 (1945), 
72, 405 (1946). 

7H. Vogt: Beitrag zur Berechnung schiefwinkeliger Platten. Diss. T. H. Hannover 1939 
(veroff.: Druck v. Charles Coleman, Liibeck). 

8 V. P. Jensen: Analyses of Skew Slabs. Bull Univ. Illinois No. 332 (1941). 

® H. Favre: Contribution 4 l’étude des plaques obliques. Schweiz. Bauzeitg. 120, 35, 51, 60 
(1942). — Le calcul des plaques obliques par la méthode des équations aux différences. Abh. 
Int. Ver. Briickenbau u. Hochbau 7, 91 (1943/44). 

10 N. J. Nielsen: Skaevvinklede Plader. Ingeniorvidenskabelige Skrifter Nr. 3, Kobenhavn: 
G. E. C. Gad, 1944. , ; 

1 A. M. Guzman und C. J. Luisoni: Elastizidad y Plastizidad de Placas Oblicuas. Univ. 


Nac. La Plata, Publ. Fac. Ci. Fiscicomat., ser. 2, 4, 452 (1953). 
122 P. Rongved: Uber das Randwertproblem schiefwinkeliger Platten. Diss. T. U. Berlin 1945. 
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entwickelte eine eigene Methode mit genahert erfiillten Randbedingungen. Grau- 
denz!3 und Heinen" haben nach dem Verfahren von Schultz-Grunow EinfluB- 
flachen ermittelt, welche den Randbedingungen punkteweise geniigen. Anisotrope 
schiefe Platten haben Naruoka und Yonezawa” mit Differenzenrechnung und 
auch auf dem Versuchswege behandelt. 

Mathematisch strenge Untersuchungen wurden fiir die Parallelogrammplatte mit 
gleichmaBig verteilter Normalbelastung 1 17 18 und ebensolcher Schragbelastung’® 
begonnen, jedoch vorzeitig abgebrochen. Lediglich Anzelius’* fiigt unvermittelt ein 
numerisches Ergebnis fiir den Torsionsmomentenverlauf an, welches aber nicht mit dem 
Resultat von Jensen” iibereinstimmt. SchlieBlich sei noch eine Dissertation erwahnt, 
aus welcher jedoch nur einige prinzipielle Erkenntnisse*! veroffentlicht sind. 

Im folgenden wird fiir die allseits frei drehbar gestiitzte Platte nach Abb. 1 die im 
Rahmen der Kirchhoffschen Plattentheorie mathematisch strenge Lésung fiir die 
Durchbiegung ermittelt. Sie ist in Form einfach unendlicher Reihen als Funktion der 
schiefwinkeligen kartesischen Koordinaten «, y angegeben. Gewisse Reihenkoeffizien- 
ten des homogenen Losungsanteiles sind jeweils aus Gleichungssystemen zu berechnen. 
Fiir zahlenmaBige Auswertungen in der Umgebung des Lastangriffspunktes empfiehlt 
es sich, von der partikularen Lésung bekannte Singularitaiten abzuspalten. 


I. Allgemeine Beziehungen in schiefwinkeligen Koordinaten 


Betrachten wir neben dem schiefwinkeligen Koordinatensystem x, y mit dem 
Neigungswinkel w zwei rechtwinkelige Systeme x,, y, und 2, y,, von welchen jeweils 
eine Achse mit einer des x,y-Systems zusammenfallt* (Abb. 1), so gilt 


Hp) 
Jz A7 t= a+ cy, Y= BY, (1a) 
Loum 62s Yo = CX + Y; 
Pijwead @ditd (2 fie 
Ope OES oy, s\oay at 
(1b) 
peop fe mals ’ 
Oyo oy’ d%, 8 \dau ay)’ 
wobei 
S=sinw, ¢€ = cosa. (2) 


* H. Graudenz: Beitrag zur Berechnung der MomenteneinfluBfelder schiefwinkeliger Platten. 
Diss. T. H. Hannover 1948. 

“4 R. Heinen: Beitrag zur Berechnung von Einflu&flachen schiefwinkeliger Platten. Diss. 
T. H. Hannover, 1948, veroff.: Ing.-Archiv 26, 268 (1958). 

Joe Naruoka und H. Yonezawa: Uber die Anwendung der Biegungstheorie orthotroper 
Platten auf die Led ates schiefer Balkenbriicken. Bauingenieur 32, 391 (1957). 

#6 A.Anzelius: Uber die elastische Deformation parallelogrammférmiger Platten. Bauinge- 
nieur 20, 478 (1939). ‘ Sse ce ee 
17 P. Lardy: Die strenge Losung der schiefen Platte. Schweiz. Bauztg. 67, 207 (1949). 

18 J. Krettner: Beitrag zur Berechnung schiefwinkeliger Platten. Ing.-Archiv 22, 47 (1954). 

HON Mr Mc Kenzie and M. Rothman: Note on the Bending of a Finite Parallelogrammic 
Plate Under Continuous Non-Normal Loading. Journ. Roy. Aero. Soc. 60, 134 (1956) 

20S. FuBnote ® auf Seite 121, 

21 G. Rumpel: Uber das Verhalten dinner Platten in den Eck i i 
coe: ckpunkten. Bauingenieur 33, 
‘ : Die positive z-Achse ist in die Bildebene hinein gerichtet, es handelt sich also um Links- 

ysteme. 
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Somit lautet der Operator A 


und die Plattengleichung 


aw aw 


1 [ o4w 
AAw= st Ee 40(sog ae ari) a 
Hierin bedeuten 


Eh 
12 (1 — %) 
die Biegesteifigkeit der Platte, » (x, y) ihre 
Normalbelastung pro Flacheneinheit, w die 
Durchbiegung senkrecht zur Mittelflache, 
h die Plattendicke, # und yw den Elastizitats- 
modul und die Querdehnungszahl des Werk- 
stoffes. 

Die SchnitigréBen (Biegemomente, Dril- 
lungsmomente, Querkrafte, Ersatzscher- 
krafte) je Langeneinheit kénnen sowohl 
fir das 2, y,-System als auch fiir das 
Xo, Yo-System aus einer bereits vorhandenen 
Losung w (x, y) durch Differentiationen nach 
x bzw. y erhalten werden. Wir brauchen 
nur in die bekannten Ausdriicke ftir recht- 
winkelige Koordinaten — wie z. B.: 


k= 


0*w 
m,,= —K |e ca P+ use| 


— die Gln. (1b) und (3) einfiihren und be- 
kommen 


im x,.y,-System (Abb. 2a) 
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ot ; 
2 (1+ 20) sare + Se] =P. 


I; 


m,, = —K |p dw + (1— w) SI, 
K 
Ms, y; (1 ae) LL) =| 
x ale) K[a(Aw)  _a(Aw) (5a) 
dary se p0gn” el ay at | 
Q Ai 1 
Y, Ty, = ae ? 
im X,., Y.-System (Abb. 2b) 
ew o2 
m,=—K|4w—(l—w) Fe], m= —K|wdw + l— a) SI, 
KI ew e2w 
tis = — (= | eae age 5b) 
a(Aw) a _ _ (dw) (5 
x. = | Ox Cc ay ’ Ey SS oy ? 
re OM xz 2 
= Git 3, 


g* 
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II. Partikulire Lésung 


Die zu untersuchende Platte ist in Abb. 1 dargestellt. Im Punkt mit den Koordi- 
naten a = u, y = vgreife eine Einzellast P an. Diese kann formal durch die Fourier- 
sche Doppelreihe”? dargestellt werden 


p (ey) = — SD sin cn sin B,v SiINnGy xX SMPrY, Cm=a= ictal 9 b,=b = “=. 

(6) 
Ein partikulares Integral w der Gl. (4) mit p (a, y) nach (6) findet man in einfacher 
Weise mit Hilfe des Ansatzes*’ 


m= lon 


w= >) > (Cn, Sinax sin By + Dn, cosax cos By]. (7) 
Die sich ergebenden Konstanten Onn und Dmn KOnnen mit 
b Sie esa A 2.62) Ter vem Tye 4cv(1 + v2 m2) m (8) 
Y= an? = TT + 20? (1 — 202) m? + 8 m4? = TI 4 202 (1 — 202) m2 + 4 m4 
geschrieben werden 
oe 4Ps4 Y : oo 4Ps* : : 
Cnn = Kapp C Sinau sin Bv, Dn n= — Happ P sin au sin pv; (9) 


Indizes bei », C und D wurden trotz Abhangigkeit von m und n weggelassen. Die Parti- 
kularlo6sung (7) lautet damit 


Wo =e 2 sin p vsin p yd, Osinau sinux— X'sin B vcospy S'Dsinau cos ax]. 

(10) 
Sie ist natiirlich in m, a, x, u einerseits und n, B, y, v anderseits symmetrisch. Davon 
kann man sich sofort durch Multiplikation des herausgehobenen Faktors 1/84 in die 
Ausdriicke fiir C und D und Einsetzen von y iiberzeugen. Die Darstellungsweise nach 
Gl. (10) erscheint hier zweckmaBig, weil die Doppelreihen im folgenden nach der (will- 
kirlich gewahlten) Laufzahl m summiert werden sollen. Es handelt sich also um die 
Summen von Reihen der Form 


> C cos at, > Dsinet. 


m=1 m=1 


Zu ihrer Ermittlung wird die Hilfsreihe 


Be i cos mT (11) 
m=1 


xt + 22 m2 + v4 m4 


beniitzt: 
=’ C cosmr = @ — 8 ct ; 
(12) 
: afl e@ a@ 
=D sin mt = 4cyr Flare 3 G2) =f | 
Nach Ausfithrung der Operationen (12) ist [vgl. Gl. (8)] 
y= bilan, 2? = 2 y2 (1 — 2c2), x= 1 (13) 


2 Vel. K. Girkmann: Flachentragwerke. 4. Aufl. Wien: Springer-Verlag 1956, Gl. (506), 
(514). 
*3 Vgl. FuBnote!’ auf Seite 122. 
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zu setzen. Die Summation der Reihe (11) mit x = 1 wurde bereits friiher?4 durchge- 
fiihrt. Auf dem gleichen Wege ergibt sich hier, giiltig fiir 


Ost 2a, Barve, (14) 
mit 
Biss oe oy 
ova py SS weNt 
1 (15) 
ee Sa Cc ane 
die Beziehung 
pi gens ze , 
eer gay er ep? Sof yt sin dt — 9 Sin yt cos #rt + 
1 3 
+ Goplyp a —cos2oa L(y sin 20a + 8 Gin 2y2) Cojyt cost + 
+ (8 sin 29x —y Sin2 yx) Sin yr sin dr]}, (16a) 
oder 
1 n 1 
= Bat Davy 0 Gj2pa codon 


-{[8 Sin yx cosda + y Cojya sin dz] Coj y(t — x) cos d (rx — 2) + 
+ [8 Goj yx sin da — y Sinyxcosdx] Giny (tr — 2) sind (tx — m)}. (16b) 


Die beiden Ausdriicke (16a, b) sind identisch, die zweite Anschreibungsart erweist 
sich jedoch als zweckmafiger fiir die Weiterrechnung. Aus (12) erhalt man mit (16), (15) 
und (13), nach Aufwendung einiger Rechenarbeit, geschlossene Ausdriicke fiir 
>) C cos mr und >) D sin mt. Ersetzen wir in diesen mt durch a ¢, dann ist 
a 

— = pc 
qo? yu=Bsa, da = poa, (17) 
y(t —a)=Bs(t—a), b(r—xz)=fc(t—a) 


und die gewonnenen Beziehungen lauten 


»' Ccosat = —5+ 


(oeszal 


a) {os (a) la = 363), (a) (a + As) ¢ (a) 88 Ayo. (a)| 4 


+ $2 (t — a) l(a — 20?) 9, (a) + (3 “fo As) (a) + 8? Ay % (a)| sts 


+ (1-2) [2 ¢- © 2 @ + a) pe @}, (18a) 


a 


3 Dsin at = — 5% Ay {o5 (¢-- a) [ps (@) — (peg + 40) (@)| — 


m=1 


— 9 ¢—a)|g2 (a) — (gpg + A2)o0 (@] + 
+(1—2) | ¢-9@—gt—ae(a|}, (18b) 


24K, Girkmann und E. Tungl: Naherungsweise Berechnung eingespannter Rechteck- 
platten. Osterr. Bauzeitschr. 8, 47 (1953). 
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mit dem Giiltigkeitsbereich 


Oat ata, (18¢) 
Hierbei wurden folgende Abkiirzungen beniitzt 
1 1 
Ay= oj 2B sa — cos2 fea = 2 (Gof? B sa — cos? Bea) ’ (19a) 
2sin2fca 2 Gin2 Bsa 


A 


Ai Coj2Bsa—cos2fhca’~?2 Gof 2Psa—cos2pea’ 


y, (d) = Cof Bsdcos Bed, ws (4) = Coj Bsd sin Bed, 

gv, (d) = Gin Bsdsin Bed, gy, (d) = Sin Bsd cos Bed. 

Die Richtigkeit der Gln. (18) kontrolliert man am einfachsten durch Entwicklung der 
Funktionen 9; (¢ — a) und (t — a) q; (¢ — a) in Fourierreihen mit der Periode 2a. 
Mit Hilfe der Beziehungen (18) kénnen nun die geschlossenen Ausdriicke fiir die 


m-Reihen in (10) angegeben werden. Wegen (18c) miissen sie fiir die beiden Platten- 
bereiche x < uund x = w gesondert ermittelt werden: 


td C cosa (te — 2) ae ee 


(19b) 


DEY Csi ; =— »' Ceose (ut x is 
~ NE EA OTL - OS iat 4S Cos CUE — th) oa - ee 


20 
+)’ Dsina(u— x)... ¢ 54, eM 

2>'Dsinaucosax = 5) Dsina (u+ 2) . 

= ry AND ano (a ye. Pe, 


Um die erhaltenen Ergebnisse einfach anschreiben zu kénnen, sollen die nachstehenden 
Abkiirzungen eingefiihrt werden. Hierbei bedeutet der Gro8buchstabe U, dai die 
Ausdriicke von der Koordinate u des Lastortes abhangen, die zweigestrichenen GroBen 
gehoren immer zum Bereich x < wu und die eingestrichenen zu x = uw. 
iby = —2 |, (u — a) pg (2) + Gs (u — a) v3 (2)] = Cj Bs (2a — u) cos Bou — 
— cos fc (2a — u) Cof B su, 
Tz = 2 [ps (u — a) 94 (a) — pa (w — a) 95 (a)] = Sin B s (2a — w) sin Bou — 
— sin Bc (2a — u) GinB su, 
U,=—-2 ile 20) [ps (u — a) Pr (2) + ps (u — a) po (a)] + 
+ (1 — ua) [pz (& — a) 3 (4) +p (u — @) 4 (2) ]} = 
= (u/a— 2c?) SinB s (2a—u) cos Bcu + (uja — 28?) cosBc(2a—u) SinBsu, 
= 2{(1 — 2c?) [p(w — a) yy (a) — oy (u — @) Y (a)] — 
— (1 — ua) [py (u — 2) ps (a) — G2 (u — @) Y (a)]} = 
ae (uja — 2c?) oj Bs (2a — wu) sin B cu + (w/a — 2 s*) sin Bc (2a — u) Cof Psu, 
U,=—2 {Ps (w — a) py (4) — py (uw — @) py (a) + (1 — u/a) [py (w — a) ps (a) — 
— P2(u — 4) G4 (a) ]} = 
—. (uja — 2) Goj B s(2a — u) sin B cu + (u/a) sin Bc (2a — u) Gof B su, 
Tg = 2 {4 (w — 2) gy (a) + 5 (w — 2) (a) — (1 — ula) [2 (w — @) 3 (a) + 
+ Yi (U— a) % (a)]} = 
= (u/a — 2) Sin Bs (2a — u) cos Bou + (ula) cos Bc (2a — u) SinBsu. 


U, 
(21a) 
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Die U; gehen aus den U. i hervor, wenn w durch u + a ersetzt wird, also etwa 


Pe=(uja—t) SinBs (a—u)cosBc (a+ u) + (ula + 1)cosBc (a—u) SinBs(a+u). (21b) 


Die gesuchten Ausdriicke lauten nunmehr 


fiir 2 <u: > Csinausnas = 
™m 


| 

= Fe 4o{|(seq + o* 4s) U,—s* A, U,+0,|o. (~) + | 
| 

) 


+ (Fea +? 41) Bs + st 41 B+ Be] 04 (2) — (ula) [Trpr(e) + Froa(ay, 


> Dsinaucosa x= 
m (22a 


= Fa Aol |(geq + 42) U,+ Us| 0: (2) — | 
a (ga + A,) U; e Te | 2 (x) + (x/a) [Oy P3 (x) — O20 (x)]} 


firzy > u: >’ Csinausin «ax — 


+ 
oe (peg + 6? As) Oa + 9% Ay ete 
oe 


“ (0, o1 (@ — a) + Uy, (4 —a)] 


> Dsinau cosax = (22b 


a 


= Te 4o{[(a0 + A.) Cs si Ts| ns (=r a) — 


a (sat A,) Cs =i T.| 2 (x — a) += [Ui9s (x — a) — Us (x—a)]}\. 


Fihren wir diese in die Partikulaérlésung (10) ein, dann entsteht 
fir x <u 


w= 9 5 SAE? 4 fx [U, SinB sa cos B (cx + y) — U, oj Bsxsin B (cx + y)| + 


ne °\a 


: (23a) 
+ 0, Gin psxsinp (cx + y) — Co Bsxeosp ex + yl}, 
firz > u 
yO AP 4,17" 10, Gin Bs (« — a) cos B (¢ (z — a) + y) — 
Wi — Ty Gof Bs (x — a) sing (c(w —a) + y)J +. re 


as U, Gin B s(x — a) sin B (c (w —a) + y) — 
ane Hie Co} B s (% — a) cos B (c (ea) + 9)}, 


wobei 
Ie b 
0 =e oe 
und 
= ] = = == 1 rr rr 
U,=(3,5+ Ay) 0, + Us, O.=(geg+ 4s) U,+ U; (25) 


Gl. (25) gilt auch fiir die U;. 


Die Ausdriicke fiir w,) nach Gl. (23a, b) gelten in y-Richtung iiber die gesamte 
Plattenbreite b, in z-Richtung hingegen sind die beiden, lings « = u zusammenstoBen- 
den Bereiche zu unterscheiden. Dieses Ergebnis wurde erhalten, weil in Gl. (10) will- 
kiirlich nach der Laufzahl m summiert wurde. Hatten wir die Summation nach der 
Laufzahl » durchgefiihrt, dann waren zwei analoge Ausdriicke fiir w, entstanden mit 
den Giiltigkeitsbereichen y < v bzw. y = v. Man erhiilt sie aus den Gln. (23a, b) durch 
Vertauschen von 2, u, a, a, m mit y, v, b, B, n: 


firy <v 


Oe) Ss Bly (V. Sin a sy cos a (x + cy) — V, Coj asysina (x +cy)|] + 


+ V, Cina sysina (x + cy) — P, oj asy 00s a (x + ey} - 
(23c) 


m2 


tty = QSOS" BL" [Pe Gin es ly — B)c08 a (2 + o(y —8))— 


ia Gof as (y — b) sina (w+ ¢(y — b))] + 
-{- V, Sin as (y — b) sina (x + ¢ (y — d)) — 
— Vz Coj as (y — b) cos « (# + ely —8)). (23d) 


Hierbei sind die B;, V;, V; analog den A; [GI. (19a)], U; und U; [Gl. (21a, b) und (25)] 
zu bilden, also z. B.: 


hae il tee 1 Es = 
By = of 2asb — cos2ach ” Vs =(3 at B,) Vago Ve, 


V, = Gof as (b — v) cosac (b + v) — cosac (b — v) Cojas (b +2). 


(26) 


III. Randbedingungen, Randwerte der partikuliéren Lésung 


Die Platte nach Abb. 1 soll an allen vier Randern frei drehbar gelagert sein. 
Unter Beachtung der Gln. (5) konnen die zugehérigen Randbedingungen (Durch- 
biegung und Biegungsmoment gleich Null) auch im schiefwinkeligen Koordinaten- 
system in der Navierschen Form 

w=0, Aw=0 (27) 


geschrieben werden, 
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Zur Berechnung der Ausdriicke fiir A w ist nachfolgende Zusammenstellung dienlich 


Tabelle 1 


a es 


w Aw AAw 


————— eee 


« Go} B sa cos B (cx + y) (2 B/s) Sin B sx cos B (ca + y) 0 
w Gin BsxcosB(cx+ y) | (2B/s) Coj Bsa cosP (ca + y) 0 
x Goj BsxsinB (cx + y) | (2B/s) Sin B sx sin B (cx + y) 0 
a Gin BsusinB (cx + y) | (2 B/s) Cof BsasinB (ca + y) 0 


sie gilt auch fiir x —a anstelle von x. Bei Vertauschen von « und y ist natiirlich B durch a 
zu ersetzen. 


Aus (23a, b) ergibt sich damit, wenn 


2% 2 Ps 
oe ee i (28) 


bedeutet, 


fir 
Aw) = Qs 2 as [U, Coj B sx cos B(ca + y) — U, SinBsasin B(cx + y)], (29a) 
Pur 2 


A w= a [U, Coj Bs (x — a) cos B (c (w — a) + y) — 
— U, Ginps(a—a ) sin B (c (wv — a) + y)]. (29b) 


Die fiir die Bereiche y © v giiltigen analogen Ausdriicke folgen aus den Gln. (23c, d); 
ihre explizite Anschreibung ist wohl nicht ndotig. 


Wir wollen nun die von der Partikularlosung verbleibenden Randwerte w, und A wy, 
bestimmen. Hierbei ziehen wir fiir den Rand x = 0 die Gl. (23a) bzw. (29a) und fiir 
den Rand 2 = a die Gl. (23b) bzw. (29b) heran. Fiir die Rander y = konst. sind dem- 
gemaiB die analogen Beziehungen zu bentitzen. Wir erhalten zunachst 


(W)x<=-0 =—Q ee oUscosBpy, (Wo)y»-0 =— Se V,, cos a x, | ; 
e 0 

(Wo)x=2 =—Q spe A, U,cosBy, (wWo)y=2=— ey, V,cosaa | as 

und 

(A wW9)x=-0 =Qa ahs U,cosBy, (Aw»)y=0 == se = Bo V, cosa, a 

(Atv))eme =Qa SA” Ay T, cos By, (wy)ys =O S82 B, V_cos ac. 


Diese Randwerte besitzen eine sehr einfache Form und erwecken den Anschein, daB 
man bei Uberlagerung der homogenen Lésung w, und Erfillung der Randbedingungen 
(27) mit w =w, + w, kaum nennenswerten Schwierigkeiten begegnen wird. Bei der 
Weiterrechnung stellt sich jedoch heraus, daB das erhaltene Gleichungssystem fir die 
Integrationskonstanten keine allgemeine Losung ermédglicht und zur numerischen 
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Behandlung wenig geeignet erscheint. Wesentlich zweckmaBiger ist der im folgenden 
beschrittene Umweg iiber eine Koordinatentransformation. 

Damit die Randbedingungen in gerade und ungerade Funktionen aufgeteilt werden 
kénnen, miissen wir das Koordinatensystem in die Plattenmitte verschieben. Die zu x 
bzw. y parallelen Achsen seien mit £ bzw. 1 bezeichnet (Abb. 1). Betrachten wir nun 
beispielsweise die beiden ersten Gln. (31) und zerlegen wir 


ASU a yy RAT SU Pere (32) 


dann kénnen wir uns die Randwerte (4w,y),=o und (4 w»),—. auch so entstanden 
denken, daB eine gerade Funktion von in € = 4+ a/2 


(A wo) =n 3 U2, cos B y (33a) 


und eine ungerade Funktion von & in = + a/2 


(A wo)u = —V Se U2, cos B y (33b) 


geliefert hat. Aus (32) erhalt man mit (19a) und (21) 


Gin Bs (a — uv) sinbcu — sinfc(a — u) Ginbsu 
2 (Coj B sa — cos B ca) : 

Gin Bs (a — u)sinBcu + sin Be (a — u) em pen, 
2 (Cof B sa + cos B ca) 


Uzg=2A, (0. gO 


(34a) 


Un=tA; (0, —U,) = 
Driicken wir weiters cos f y gemaB 
cos B y = cos B( + 5) = cos 5 * cos By — sin“¥ sin B 


durch 7 aus und beachten, daB fiir gerade n das zweite und fiir ungerade n das erste 
Glied verschwindet, dann lauten die Gln. (33) fiir die Randwerte der 4 w, in & =a/2 


(A Wo)e =Qa = a U2, COs cos B 7 — = — gs sin “— sin B n) 


eis (35a) 
(A wou =Qa |- 5 BE Us, cos cos Bn +3 ou. Us, sin“; sin B a 
In gleicher Weise entstehen fiir die Rander 7 = konst. mit 
Vog=4 By (Ve te Va), Vau=4 By (Vy ai V;) (34b) 
die Randwerte der Aw, in 7 = b/2 
(A wy)p =Qa > ana V2, cos cos aé = mma Voz sin“, sin a ‘| 
(35b) 


yy sina ah 
(Arr) 0, [— 33 ae Vou cos Sy cos wé + vi = een 5} sin a |. 


13,5 


An Hand der Gln. (35) kénnen die Randwerte der Durchbiegungen w, sofort angeschrie- 
ben werden. Man braucht bloB die Analogien zwischen den Gln. (30) und (31) zu 
beachten und erhalt 
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in € =a/2 
(e)); = |- % oo ~ Us, cos cos hn + 33 ss Us, sin sin Bn), 
4,6 ee 


(36) 


\7 sin ee 

(OP REE, B ba Usy cos" cos Bn — ay Uae sin "sin B 
74,6 45 Be 5 

mit 


aA (Os+ Us) =(5-2 Ped, ) Ug + SIPs to wl Fou + sin Be a—u )Go} Bs ) 
( 


__ 2Goj Bs (2a—vw) sin Beu + Gof Bs (a —u) sin Be (a+ u) — Coj Bs (a+u) sae (a—u) 
2 (Wof 2 Bs a — cos 2 Bea) 


2 (Coj Bsa — cos Bea) 


Usu=t A,(Us—Us) = (seat ) Oe u Woj fs ( ea ee Cof Psu 


2 (Coj Bsa + cos Bea) 
__ 2Cof Bs (2a—u) sin Bew — Coj Bs (a—u) sin Be (a+u) + Cof Bs ‘Bs (a+) sin Be (a—u) 
2 (Col 2 Bsa — cos 2 Bea) 


t 


IV. Homogene Lésung 


Wir brauchen nun eine Lésung w, der homogenen Plattengleichung (4), welche 
gemeinsam mit dem partikuliren Integral w, die Randbedingungen (27) erfiillt. 
Hiefiir ermittelt Lardy* die Funktionenklassen 


kexp{a[+i(a+cy) +sy]}, | 
kexp{B[+i(cca+y) tsa}, f 


wobei a und £ beliebig sind, wahrend fiir k wahlweise 1 oder x oder y zu setzen ist. 
Bei Beachtung der Gln. (1a) ist sofort ersichtlich, da8 es sich hierbei um die bekannten 
Lésungen der Scheibengleichung in rechtwinkeligen Koordinaten x,, y, und 22, y2 
handelt, die in bestimmten Linearkombinationen zusammengefaBt sind. Da wir ja die 
Randbedingungen im é,7-System erfiillen wollen, schreiben wir die Funktionen (38) 
mit € und 7 anstelle von x und y. Fiir « und f wahlen wir zweckmaBig die bereits mit 
denselben Buchstaben bezeichneten Groen [vgl. Gl. (6)]. Nun beabsichtigen wir die 
beiden Bedingungen (27) fiir gerade und ungerade Funktionen von & bzw. 7 getrennt, 
daher nur in € = +a/2 und 7 = +6/2 anzuschreiben. Fiir diese Gleichungen bendtigen 
wir also 16 frei wahlbare Integrationskonstanten und bilden daher den Ansatz 


(38) 


w, = S'|(, += = N,) cos B (c € +1) Goj Bs é+(N,+-_ Ne) cos (c£ +7) Gin sé 


n=1 


+ (WV, +=N,) sin 6 (c€ +7) Goj Bs é+(N +; Ne) sinf (cf +1) SinBas| 


i (a+ + M,) cosa (&+¢7n) Goj asy + (m,+ +. My) cosa (+ cn) Sinasy 
m=1 

+ (M,-+ 7 M,) sina (E + en) Cofasy + (Mit Me) sina (€ +07) Ginesy]. 

(39) 

Um Doppelindizes zu vermeiden, wurde fiir alle Konstanten, die den Index n fiihren 


sollten, der GroBbuchstabe J, fiir jene mit dem weggelassenen Index m der GroBbuch- 
stabe M gewahlt. 


* Siehe FuBnote?’ auf S. 122. 
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Mit Hilfe der Tabelle 1 ergibt sich aus (39) 
Aw, =22| Simi, cos (c & + 1) Sina & + Ne cosh (cE +7) Gof Bs + 
+ N,sin B (c é + 7) GinBs & + Ngsin B (c & + 7) Goj Bs ] + 


+ 2» Ym [M, cos a (é +c) Ginasy + M,cosa(é + cm) Cojasy + | 


m=1 


+ M,sina (& + cn) Sines + My sin e(E +1) Gof esr) (40a) 


In A w, sind also nur halb so viele Konstanten enthalten wie in w,. Wir werden daher 
zunachst die zweiten Bedingungen (27) erfiillen. 

Mit Beniitzung des Additionstheorems der Winkelfunktionen kénnen wir A w, in 
gerade und ungerade Funktionen der Koordinaten ¢ und 7 aufspalten: 


ZA{ Sin [We Goi fs € cos fof +N, Gin s Ein fe £) cos Bn + 
n=1 


Aw, —~abs 


+ (N, oj Bs EcosBcé — N, Sin Bs Esin Be &) sin By] + 
+ »'m [M, Cojasyncosacn + M, Ginasysinacyn + 
m=1 


+ My Gin 97 c08 eo + Mg Gof es sin ee] eos « &| + 


2 . 
+24 > n{[(N; Gin Bs Ecos Be &é + Ng Coj Bs EsinBc &) cosBy + 
+ (N, GinBs EécosBcé — N, oj Bs Esin fs &) sin By] + 
+ 2») m [M, €oj asncosacy — M; Ginasnsinacn+ 
m=1 

+ My Gin as cos eon — Me oj es sin wen) sine (40b) 
Hierbei ist die erste geschlungene Klammer in é gerade, die zweite in  ungerade und 
auBerdem sind sowohl in den n-Reihen als auch in den m-Reihen die in 7 geraden und 
ungeraden Funktionen zusammengefaBt. Fiir € = a/2 ist a € = m 2/2 und es verschwin- 


det fiir jedes m entweder cos a é oder sin a &. Wir erhalten daher mit der gleichen Be- 
zeichnungsweise wie in den Gln. (35) fiir die Randwerte der A w, in = a/2 


9 co 
(Bene = 34] & (Ve Co1®$2c0 fee. Sin PS sin ?S") nm cos Bn + 


+ 2” (M,Wojasn cos acy + M, Sinasy sin acy) m cos +- 


2,4,6 


+| (Ws Coj 5" cos "$4 — N, Sin Ps sin “S) m sin B + 


+ (M, Sin asy cosacyn+ M, Go} sy sin on) m cos "| ey 
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2 \ 
(Aw). = Fj] S (me Gin ton ees Gof "5 sin 5") n cos By + (41) 


+ wo” (M, oj asyn cos acy — M, Sinasy sin ae) sin"="| + 


oS »' (Ny, Sin eS N, Gof M5 sin *S*) n sin B 7 + 


a 
3 


(M1, Sin asy cos «ey — M, Goj asy sin wen) msin "™|t 


In den beiden Gln. (41) enthalt jeweils die erste eckige Klammer die geraden Funktionen 
von 7. Die Randwerte der A w, in 7 =b/2 sind vollkommen analog und brauchen daher 
nicht angeschrieben zu werden. Mit den Ausdriicken (41) und (35) kénnen wir nun die 
Randbedingung 

Aw, + Aw, =0;-in é = a/2 (42) 


anschreiben. Sie zerfallt in 
(A wo)e + (A w,), =0, (A wo)u + (4 w,)u = 0 


von welchen jede Gleichung selbst wieder in einen geraden und einen ungeraden Anteil 
von 7 aufzuspalten ist. Die Randbedingung (42) liefert somit die vier Beziehungen 


gue Us, cos 7008 n+ S'(N, Coj Pee cos oN) Sin 4 sin OS) cos By | 


2,4,6 


+ 3» (M,Cojasy cosacn+ M, Ginasy sin ac) m cos~5- =0, 
25,6 


aad Gee US sin“ 5 sin Bn +2", oj Coca Sin PS sin 5") nsingy 
1,3,5 
“- py M, Ginasyn cosacy + M, Cojasy sin a ¢ 7) m cos > = 0e (a3) 
2,4,6 
a — Ua: cos ~~ cos By i » (M. Sines cos $4 + Ny Gof 5" sins) n cospy 
ANG 


a5 = (M,ojasy cosacn—M, Sinasy snacyn)m sin =0, 


g srinee Up, sin "sin py +S (Wy Gin #2 cos 2&4 _. NV, Coj 7S" sin 5) nm sin B 7 


AG 


\ 


+ 3” (M, Sinasn cosacyn — M, Gojasy sin acy) msin ~~ =0. 
1, =, 

Analoge vier Gleichungen erhalten wir bei Anschreibung der Bedingung (42) fiir 7 = 
= 6/2. Diese kénnen auch direkt aus den Gln. (43) abgelesen werden, wenn 7, b, », 
n, B, N, U mit é, a, u, m, a, M, V vertauscht wird. Insgesamt besitzen wir also acht 
Gleichungen fiir die acht Unbekannten M,... Mz, N;... Ng. 

Bevor wir ihre Auflésung behandeln, mollen wir die weiteren acht Gleichungen fiir 
die Unbekannten M,... M,, N,...N, anschreiben. Sie entstehen aus der ersten 


Randbedingung (27), also aus 
Wo + w, =0mé =a/2. (44) 
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Die hier bendtigten Randwerte von w, sind in Gl. (36) angegeben, jene von w, sind aus 
(39) in gleicher Weise, wie eben die Randwerte von A w, gewonnen wurden, zu ermitteln. 


Wir erhalten dann die vier Gleichungen 


—Q pe Ua O08 cos By + 
2,4,6 
+ S'|(w, oj 4 aN, Gin 282 *) cos 5° + (MN, Gin 2 + 4N, Coj 84 Pet) sin "| cos By + 
n=1 
+ 3* (a Goiesn+ £m, inasn) cos wey + (My Sinasy + 5M, Gol an) sin acn] cosy —0, 
Q SP? y,, sin sin Bn + 
1,3,5 
+ »'|(a, Gof S* +40, Sin 5") cos 5" — (Ne Sin 984 4 4, Cof S*) sin | sin By + 
n=1 
+ 3 [(Ma Gof asm +} M, Sinwsn) sin wer + (My Sines y+ Ms Coj a7) cosacn|cos"y = 0, 
Q ane ange U 3.008 “cos Bn + rf 
are » |(¥s 60) Eero, Sin 5°) sin 5" + (We Gin 82 44N, Coj P34 *) cos 5 *| cos Bn + 
+ Sslleccotcon +44; incor concen ~ (afuGineon +] 3 Go eo) nae] sin =o 
57 sin Bo nn, 
— ee 7a Ueu Sin sin B y — 


Lys thye § 


[27 Gof 5% + 47, Gin 54) sin 5* — (av, Sin 5" + 4, Gof 5") cos “54 | sin By + 


(a0, Gofacn + 5 Me Sinasy) sine oy — (M, Cinasy+7M, Coj asm) cosacn| sin" =0. | 


= 

. 
& 
oO 


Die vier analogen Beziehungen fiir 7 = b/2 kénnen wieder, wie bei den Gln. (48), direkt 
abgelesen werden. 

Die Auflésung der Gleichungssysteme (43) und (45) erfolgt mittels Keattisatee 
vergleich. Hiezu miissen nachstehende Funktionen von 7 in Fourierreihen mit der 
Periode 2b entwickelt werden 


Coj asyn cosacy =» Cn’ cos B , + Gof asn sin acy = Ss gn cos Bn, 
£8: 13,5 
Cin asy sinacy - 3 d,’ cos B 1, + Gin a8n COS ACH Sel h,’ cos Bn, 
a a (46) 
Coj asyn sin acy = 2 é,’ sin Bn, F Gof a8 COS aCH = “ k,’ sin B n, 
8, 1, 8,5 
Sin aS COS aCH -2, f,’ sin B n, + Sin asy sin acy =» L,’ sin B n. | 
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a ae von geraden Funktionen (yj) mit verschwindendem konstantem 
ie 


= >’ a, cos B 7 
ist nur in einer Form méglich, wobei die Koeffizienten aus 
b2 ; 
4 ° 
a oh | ” (n) cos Bnd ny — ((6 —») cos B 4 dy 
9 b/2 


zu berechnen sind. Mit 
cos 8 7 = cos B (yn — b) cos B b — sin B (yn — b) sin Bb 
und unter Beachtung von f 6b = nz, kann fiir a, geschrieben werden 
6/2 
an =, | 9 (n) cos Bn dnd Gn’ = 1ROEB de (47) 
0 
Bei den ungeraden Funktionen y (7) wurde von den beiden geeigneten Fortsetzungs- 
moglichkeiten mit 


= >’ b, sin Bn 
jene gewahlt, wo 


6/2 


b, =o5 Jo (7) sin fran Jv (6 — yn) sin Byundy 


ist. In gleicher Weise wie oben entsteht hieraus 
b/2 
== | v(m) sinBndn, Waa 3) Bae (48) 


6 


Dn 


Die Fourierkoeffizienten ¢n’, dy’, gn’, hn’ sind somit als a, aus (47) zu berechnen, 
wahrend sich e,’, fn’, kn’, In’ als b, aus (48) ergeben. Die Ausrechnung liefert 


asb ae ac Ran pee 
=f Cit, sim, - ° + fy CoS 83> 
b aby acb 
dy’ =f, Gin sin “6° — fySoj=5- cos—- , 
asb ae acb OU yee 
Gi= fs Sin 5 — f, Gof cos 
5 3b acb 

i=], Sin 22" sin 2° 5 seas Sof - cos ; (49a 


asb ae acb asb acb als asb acb 


gn’ =4(f. Gof — f, Gin cos-4- —g, Sn - ® sin gb cos}, 


) 

. 
: b 8b b 

ha’ =4(f, Co 23" ein bij = eps — 9,00 —— 5 See — 9, Cof=- cos a) 
cb ; db. b sb cb 

kn’ = 4(f4 Gof SS si i ee snot 5 ee es) 


2 
asb acb asb asb as Ae 
2 


ot /pSo\-5 sit ff, Cin, cos “S— — 9,Sin 2S? sin 2S" + g 3 Coj =- 
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wobei 


nn 
SS RCE TR ee de = (fe fs — fifa) sin 9 


fo =[1 4+ 2% m2(1—20%)]F, go = (hfs + fafa) sin a 


soe (49b) 
fs =svm (v2 m*? + 1) Ff, Y3 = (fatac= [i Je) SUN soe 
PEND [a (ae tine WO oe Gu tii" — f.? + fs? — f,?) sin > 
und 
4sin b 
ie na(l 2m (l—2e)+y%mi] > ” an’ (49¢) 


Nach Einfiihrung der Reihen (46) in das Gleichungssystem (43) spalten wir zweck- 
mafig in gerade und ungerade Laufzahlen » und m auf. An der Bezeichnungsweise 
der Konstanten NV; bzw. M; mu8 somit erkenntlich sein, ob sie zu geradzahliger oder 
ungeradzahliger Laufzahl gehéren. Analog Gl. (46) beniitzen wir einen bzw. zwei Striche 
fiir ungerade bzw. gerade Laufzahlen. Ks ist also 

N, ie turn = '1,-3,.03. 5; fa’ fir m = 1, 3,.5..., 


M: 


DIN far nd 6 a8, i |My,’ fiir m = 2, 4, 6.... (80} 


Aus der ersten Gl. (43), die nunmehr nachstehende Form besitzt 


co 


Q S22" V7, cos 008 Bn + 4 DS (ey Mo” + dy’ Mr’) m cos ™S* cos Bn + 


1, 3,5 2,4, 6 


2,4,6 
a aE ( Ne Gof bee cose" + Ny” Gin ==" sin 5") n cos B+ 


»4,6 
oe Pe Coj Eat cos N,’ Sin Pe sin 7S") m, cos Bn == 0 
entstehen somit die erste und fiinfte der folgenden Beziehungen 


Q 


sin oe 


Us, cos ~ 5 tn (Ne CofE" cos 84 + Wy” Sin 284 si 2 ns") =n 


n ee Coj ?S" co ee N35” Sin 84 sin 2&4 l=; 


sin ae 


—@ Us, cos Sar te n(Ns” Sin $4 cos 4 + Ne” Coj = sin “*) = I 


n (1, % See Neo’ Gof "5" sin 5) =0, 


nm (Ne Gof 5" cos "5" + N71’ in PY sin") + OA ( (cn’ Me’’ +d, ‘M,'’)m cos "5 = 0, 


sin Bv 


cee ne eee = HOE Dre eine Ps sin 5) + 


2 
ae a (en Ms” + fx’ Ms’) m cos“>- =0, (51) 


2,4,6 be! | 
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B co 
m (No! Gin *$ cos 924 +. Ne’ Gof 88 sin Pst) +. 3h (oy! Me’ — dy’ Me’) m sin” = 0, 


Z 1, 3,5 


| 
Ur, sin“y 5} ~ + n (Nr ‘ Sin a cos #4 No’ Coj = sin “54 


Q 2 


= . “ NVIE 
+ 3 (6n’ My’ — tf,’ Me’) m sin —>- = 0; 
ra b 


die weiteren sechs angeschriebenen Beziehungen sind aus den restlichen drei Gln. (43) 
hervorgegangen. Wie bereits erwaihnt, haben wir noch ein analoges, zum Rand ee 


gehoriges Gleichungssystem zur Verfiigung. Die zweite und vierte Gl. (51a) und die 
dazu analogen Gleichungen liefern 


N,” = N,” Ty — gee M,’ = M,” X 28) tg ae 
(52) 
aes acb 
2 


a 


Nee erdrag a Te es 


Bei Beniitzung der ersten und dritten Gl. (51a) und der beiden analogen Beziehungen 


k6énnen wir damit die Konstanten V,”.. NV,’ und M,”.. M,’’ explizit ausdriicken 
ae Bea 
Nn.” sin Bv N Me sin fv nn 8 9 
Bae ViGien tae isa au, Ng = @ Pee Bears 782 
om Co} Ta 
ca 5 ca 
sin ie ee sy sin pv N It 2 
cats Q U U> N — Q Cos 
2 isa 2 ? n2 2 Bees ( 
Co) = Cin 
ae 
sinaw ma © 
M,” =@ ms COS > aed au 
Sin-—- | 
Die Konstanten mit ungeraden Laufzahlen (N;’... N’,, M;’... M,’) miissen aller- 


dings aus einem gekoppelten Gleichungssystem numerisch berechnet werden. Dieses 
Gleichungssystem besteht aus den Gln. (51b) und den dazu analogen Beziehungen. Die 
dort aufscheinenden Sumnien tiber gerade Laufzahlen sind hierbei gemaB 


co co 
WW , Hy wil sin ach . 
oF (er Me +d,’ M,/’) m cos = iia at oe, — fesin ew, 


Qu : 1+ ojasb cosach 
2|fasin ach la em ack sinau 


(54) 


co le) 
—y Mm It at 
oF (en M a abtT M “jm e os = Q > 
i 2,4,6 


durch bereits bekannte GréBen ausdriickbar. 


Bei der Auflésung des Gleichungssystems (45) gehen wir in gleicher Weise vor wie 
eben beim Gleichungssystem (43). Die Konstanten mit geraden Laufzahlen kénnen 
bei Beniitzung der Beziehungen (53) wieder explizit angeschrieben werden 
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in 62 N Tt ee oa psa 
sin pv NT 2 1 
Nee =Q ne - Cos om a Use | 2 es Tq ) 
Co] = 


pea 


P cos —— 
sin fv NI 2 


Bsa 
Na = Oe 008 = paa (Ua + # Ventas oe). 


cor 


eo 
sin pv na 2 1 psa 
Ne — Q an COS aa (U, Loe U, u Tq 5) i (55) 
a OG) 
sin fv na 5 side 
N, Q n2 89 joa( Use + U2 ,6tg 45 ar) 
ae 
ach 
sm au tite ae asb 
Mev =O cos So 5 (V 82 ee, 2 Vee 30-5 +) ovale, VE aysiye 
oj ar 


| 
n2 2 Cot psa\, 22 


Zur Ermittlung der Konstanten mit ungeraden Laufzahlen (N,’... N,’, M,’... M,’) 
verbleibt ein Gleichungssystem, bestehend aus den nachfolgend ergescenehanen und 
den dazu analogen Beziehungen 


ae een 


ee Coj 284 = So ose! +N,’ Sin PS sin ee (2, 


co 


aie >” (Cn My’+ Gig We ob lO DENS a On ee) soo BE I 0, 
2,4,6 “l 


Ny’ 6oj 724 sin ae Sin 8 co 55" 4 4(My Sin “34 sin 204 — N;' Cof"5" cos 25) 


+ Se! My — fe MYL? Me BEM) ain oe EROS Gein REP 


2 
1,3,5 a n 


Ny’ Gin 3% sin 22*_ Cope “co so 44 (v Ns/ Gof 3 sin — 


alee 

% 

[= 
bo 


— SH fy! Mal +e! Ma’ + by! My! 41, ‘M,'’) cos ae radu 


2.4.6 n2 


N, Gin “cos "S" + W,’ Gof 284 sin 224 4 4 L(V, she en Sin "sin 25") 


Fel (Cn My! — dy’ My’ + hy’ My’ — gy’ M,') sin — = 0. 
(56) 


Hierbei sind die Konstanten N,’...N,’ und M;'...Ms' aus der Auflosung des Glei- 
chungssystems (43) bereits bekanne Mit Hilfe von ( 53) und (55) kénnen ferner die Reihen 
in der ersten und dritten Gl. (56) durch bekannte GroBen ausgedriickt werden 
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Ny , wr a , r , " m 
GM eM hig’ Mg! 1g,’ Me") cos = 


2.4.6 2 | 

a vi, sin aw f sin acb 

—Q = m2 \fo Veettn Vee —Gaasb lh VesetGooth Stg asb) Vast, | 
ey. (57) 

= (fn Ms oe en’ MM," Ken’ MM.” ue By M,’’) cos 5" 53) 

ssi Vy, Sina w { sin acb Co) asb + cos acb 

=e m2 (fo + Siren ts) Vou (fy Veu t+ $94 Vou) Gin asb ie | 


Das vorliegende Problem der Parallelogrammplatte unter einer Einzellast gemaB 
Abb. 1 ist somit gelést. Die Biegeflache ist zu berechnen aus 


w(x, y)= Wot w,, (58) 


die zugehérigen SchnittgréBen ergeben sich dann aus den Gln. (5). Fiir die partikulare 
Lésung w, ist eine der Gln. (23a—d) mit den HilfsgréBen nach Gln. (2), (6), (24), 
(19a), (21), (25), (26) heranzuziehen. Die homogene Lésung w, ist in Gl. (39) durch die 
Koordinaten £, 7 ausgedriickt. Eine gesonderte Anschreibung von w, (x, y) ist nicht 
notig, da die beiden Koordinatensysteme durch bloBe Transformation ineinander iiber- 
gehen. Die Integrationskonstanten N; und M; sind unter Beachtung von (50) aus 
(53), (55) bzw. aus den Gleichungssystemen (51b), (56) zu ermitteln. Hierbei werden 
die weiteren Hilfsgr6Ben nach Gl. (34), (37), die Fourierkoeffizienten nach Gl. (49a—c) 
und die Reihen Gl. (54), (57) bendtigt. . 

(Hingegangen am 23. April 1959) 


Buchbesprechungen 


Experimentelle Untersuchungen an verschieden stark konvergenten, schlanken Rotationskérpern 
bei maSig hohen Uberschallgeschwindigkeiten. Von H. R. Voellmy. (Mitteilungen aus dem 
Institut fiir Aerodynamik an der Eidgendssischen Technischen Hochschule in Ziirich: Nr. 24.) 
Mit 37 Abbildungen, 828. Brosch. sfr. 17.60, DM 17.—. 


Der Verfasser der vorliegenden interessanten Untersuchungen hat bei Windkanalmessungen 
an schlanken Rotationskérpern im Uberschallkanal die Beobachtung gemacht, da8 die gemessenen 
Werte von den gerechneten abweichen. Er fibrte nun eine Versuchsreihe durch, um die Ab- 
weichungen von der Potentialtheorie zu erfassen. Untersucht wurden fiinf schlanke Rotations- 
kérper, die die gleiche Spitzenform, aber verschieden stark konvergente Heckteile hatten. 

In der Theorie weist der Verfasser auf die bekannten Linearisierungen hin und zeigt in sehr 
ausfiihrlichen Kapiteln die sich ergebenden Resultate bei verschiedenen Vernachlassigungen. In 
einer Reihe von Versuchen bei Mach-Zahlen 1,38 und 1,67 werden die zu untersuchenden Groen 
bestimmt und mit den berechneten verglichen. Als Ergebnis findet er ein Abweichen der Theorie 
vom Experiment bei einer Anstellung, die groRer als 3 Grad ist, wie auch bei starkerer Krimmung 
des Heckteiles. Besonders bemerkt wird noch, daf die Grenzschicht die Umstromung stark 
beeinfluBbt. 

Untersuchungen dieser Art sind deshalb immer interessant und aufschluBreich, weil sie recht 
gut zeigen, wie weit potentialtheoretische Berechnungen mit moglichst genauen Messungen tiber- 
einstimmen und anwendbar sind. R. Bruniak, Wien 


Hawranek/Steinhardt: Theorie und Berechnung der Stahlbriicken. Nach einem hinterlassenen 
Manuskript von A. Hawranek: Vollstindig neu bearbeitet von O. Steinhardt, Mit 269 Abb.., 
XII, 426 S. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1958. Geb. DM 66.—. 

Professor Hawranek hatte wihrend des zweiten Weltkrieges ein Werk gleichen Titels im Wege 
des Deutschen Stahlbauverbandes nach einem vereinfachten Druckverfahren herausgegeben. Die- 
sem Werk waren die neueren Entwicklungen des Stahlbriickenbaues, wie sie sich in den Dreibiger- 
jahren anbahnten, zugrunde gelegt. Seither hat der Stahlbriickenbau abermals gewaltige Fort- 
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sehritte aufzuweisen. Es ist daher sehr dankenswert, da8 Prof. Steinhardt dieses Werk durch eine 
miihevolle vollstiindige Neubearbeitung zu einem dem gegenwirtigen Stand gerecht werdenden 
Lehrbuch gemacht hat. Demnach gliedert sich das Werk in das Grundlegende aus Festigkeitslehre 
und Elastizititstheorie, Theorie der Platten und der Tragerroste, Stabilitatsprobleme, Fahrbahn- 
tafeln aus Stahl und Beton, Vollwand- und Fachwerktrigerbriicken, Bogen- und Hangebriicken 
sowie Verbundbriicken. 

Hervorgehoben seien die praktischen Methoden zur Berechnung von Flichentragwerken und 
yon Trigerrosten, die eingehende Behandlung der Stabilitatsprobleme sowie die Dynamik der 
Hingebriicken. Besonders wertvoll sind ferner die zahlreichen Zahlenbeispiele und die vielen 
Literaturangaben. 

Die Fachwelt wird es ohne Zweifel sehr begriif%en, nunmehr ein modernes, sich durch Klarheit 
auszeichnendes Lehrbuch der Theorie und Berechnung der Stahlbriicken zur Verfigung zu haben, 
das die vielen neuen Erkenntnisse enthalt, die bisher nur in Zeitschriften verstreut zu finden waren. 


FB. Ozitary, Wien 


Kolloquium iiber Ermiidungsfestigkeit. Stockholm 25.—27. Mai 1955. Verhandlungen. Editors 
W. Weibull, F. K. G. Odqvist. With 194 Figures. XI, 339 8S. Berlin-Géttingen-Heidelberg: 
Springer-Verlag. 1956. Geb. DM 46.50. 


Den Inhalt des Buches bilden Forschungsbeitrage bekannter Wissenschaftler aus zwolf Landern, 
die an diesem Kolloquium iiber neuere theoretische und experimentelle Fortschritte auf dem Gebiet 
der Dauerfestigkeit berichteten. In 35 Beitrigen werden in der Hauptsache Probleme aus dem Ge- 
biet der Werkstoffestigkeit gegeniiber Schwingungsbeanspruchung und die physikalischen Grund- 
lagen der ,,Ermiidung‘‘ behandelt, die derzeit noch nicht vollig geklart sind. Die einzelnen For- 
schungsbeitrage vermitteln ein Bild tiber den Stand der Grundlagenforschung hinsichtlich der 
Physik des Dauerbruches, des Ablaufes der Werkstoffschadigung in schwingungsbeanspruchten 
Metallen, der Ausbreitung des Daueranrisses, der Anwendung mathematisch-statistischer Methoden 
bei Dauerversuchen, dem Prot-Verfahren zur abgekiirzten Ermittlung der Dauerfestigkeit, der 
Dauerfestigkeit bei veranderlichen Spannungsamplituden (Betriebsfestigkeit), bei Kontaktermiidung, 
bei mehrachsigen Spannungszustaénden sowie bei hohen Temperaturen. Die Referate sind teils in 
deutscher (9), teils in englischer (19), teils in franzésischer (7) Sprache verfaBt und durch zahlreiche 
Schrifttumsangaben ergénzt. 

Ihr Inhalt spricht vor allem den Materialpriifer, Metallphysiker, Festigkeitstheoretiker. und 
Werkstoffachmann an. 

Das Buch mu diesem Leserkreis wirmstens empfohlen werden, da es einen Markstein in der 
Entwicklung der Festigkeitsforschung darstellt. K. Létsch, Wien 


Einfiihrung in die Theorie der StrOmungsmaschinen. Von A. Betz. Mit 168 Textabb., XVI, 2728. 
Karlsruhe: Verlag G. Braun. 1959. Ganzleinen DM 36.—. 


Die verschiedenen Arten von Strémungsmaschinen haben sich konstruktiv fast véllig unabhangig 
voneinander entwickelt. Dies hatte aber leider auch eine getrennte Entwicklung der zugehdérigen 
Stromungstheorien zur Folge. Der Verfasser stellt sich nun die Aufgabe, diese Theorien, von ihrer 
gemeinsamen Grundlage ausgehend, zu einer Einheit zu verschmelzen. Eine wesentliche Rolle 
spielt dabei der Begriff des Stromungsgitters, und dessen Theorie ist denn auch ein betrachtlicher 
Teil des vorliegenden Buches gewidmet. Voran geht eine Ubersicht tiber die Grundbegriffe der 
Stromungslehre inkompressibler und kompressibler Medien und ein Abschnitt tiber die Stromung 
in Rohrleitungen und Kandalen. Der restliche Teil des Buches (etwa die Hilfte) behandelt das 
Grundsatzliche der Strémungsmaschinen: Pumpen, Turbinen, Propeller, Windrider, Raketen und 
Strahltriebwerke. Ein Anhang bringt Tafeln und Diagramme. 

Mit Ausnahme des Abschnittes iiber die Gitterstromung ist das Buch durchaus elementar 
gehalten und ohne besondere Vorkenntnisse leicht verstindlich. Wie nicht anders zu erwarten, ist 
es klar und sauber geschrieben — Flichtigkeiten, wie ,,... miiBte Warmeenergie in mechanische 
umgesetzt werden, was dem zweiten Hauptsatz der Warmelehre widerspricht‘ (Seite 57), werden 
bei einer Neuauflage leicht auszumerzen sein. Vor allem unseren angehenden Ingenieuren mu8 


man den Band warmstens empfehlen; er wird ihnen einen wesentlichen Teil ihrer wissenschaft- 


lichen Grundausbildung vermitteln. 
H. Parkus, Wien 


¢ 
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Peltemnensekals: Dies iit eine tiefgreifende Um-. ae e 
zur Folge gehabt. Neue Probleme hat die Ub réragung wesentlic 
gestellt. Erstmalig sind auch Gleichstromanlagen f solche Zwecke gebi ut 
ist vieles zur rationelleren Ausniitzung der Baustoffe, zur Erhohung der B triebss 
zur Verkleinerung der Abmessungen geschaffen worden. Neue Isoli rer 
und haben zum Teil gro8e praktische Bedeutung erlangt. Die Forschun i 
-Gebiet der Koronaerscheinungen an Freileitungen, der Physik. des Sto8 
der Blitzvorginge, des Lichtbogens und vieler anderer Vorgiinge erzielt, die 
technik ihre praktische Anwendung gefunden haben. Die Leistungen ' 
bei verhiiltnismabig kleineren Abmessungen um ein Mehrfaches vergr 


Jahren. Auch bei der Erzeugung hoher Gleichspannungen, i in der Réntg rohr ch 
Bau von Hochspannungsapparaten fir die Atomphysik sind wesent be Fortschritte Z 
-zeichnen. Alle diese Neuerungen sind in der vorliegenden vierten J sat 
sichtigt worden, so da die Rothsche ., Hochspannungstechnik“ als § dardw er. den 
lange geworden ist, den modernsten Aueprapien ae wird. 2 A aS 
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_ Das elektrische Feld. — Die festen Isolierstoffe der Hdphapantumetogas _ Fa 
stoffe der Hochspannungstechnik. — Die Luft und andere Gase im elektrischen 
Zusammenbau mit festen Baustoffen. — Luft im Zusammenbau mit festen Stott 
magnetische Schwingungs- und StoBvorgiinge. — Wirkliche Anforderungen an di 
Festigkeit im praktischen Betrieb und Prifvorschriften. — Einrichtungen fir d 
prifraum. — Der Lichtbogen in Luft und unter Ol. — Hochstromprobleme in Ho 
anlagen. — Hochspannungsanlagen fiir Wechselstrom. — pte a ees ] 
— Literatur- und Sachverzeichnis. 
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